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ABSTRACT 

A  thermal  design  procedure  for  air  cooled  electronic 
circuit  boards  has  been  developed  for  the  Hewlett-Packard 
ilodel  9845  desktop  computer.   The  system  of  interactive 
programs,  called  THERMELEX,  performs  thermal  analysis  of 
printed  circuit  boards  to  predict  either  the  junction  tem- 
peratures for  given  power  dissipation  levels  or  the  maximum 
power  levels  for  given  junction  temperature  limits.   The 
system  includes  the  following  features:   totally  interactive 
with  all  input  in  question  and  answer  format,  simple  data 
verification  and  correction  capabilities,  ability  to  store 
and  retrieve  circuit  board  descriptive  data  totally  under 
program  control,  wide  variety  of  output  formats  including 
tabular  and  graphical.   By  using  internal  selection  of  heat 
transfer  correlations,  the  THERMELEX  system  depends  only  on 
input  of  physical  parameters  for  thermal  predictions. 
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of  circuit  board  picture  on  screen 
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Xi     Distance  from  the  entrance  of  the     [m] 
cooling  air 

Y      Direction  perpendicular  to  cooling 
air  flow  also  referred  to  as  the 
vertical  direction  due  to  orienta- 
tion of  circuit  board  picture  on 
screen  of  computer 

Zb     Distance  between  circuit  boards       [m 
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I.   INTRODUCTION 

A.   BACKGROUND 

Electronic  components  generally  convert  a  significant 
fraction  of  the  input  power  into  internal  joulian  heating. 
When  the  components  are  large,  widely  separated,  and  air  is 
permitted  to  circulate  freely  among  them,  this  heat  generated 
within  the  component  is  transferred  to  the  environment  through 
natural  convection.   Indeed,  millions  of  pieces  of  home  elec- 
tronic equipment  have  operated  reliably  for  years  depending 
only  on  natural  circulation  for  cooling.   The  major  emphasis 
over  the  last  15  years  has  been  a  continuing  effort  toward 
a  reduction  in  physical  size  of  components  and  increased 
component  density  within  electronic  equipment.   The  military 
is  particularly  in  need  of  smaller,  more  complex,  yet  relia- 
ble equipment  that  must  often  be  sealed  from  an  extremely 
hostile  environment.   This  results  in  the  need  for  more 
efficient  cooling  methods. 

The  age  of  large  scale  integration  (LSI)  is  here. 
Although  the  power  dissipated  in  each  active  junction  within 
a  component  has  greatly  decreased,  the  large  number  of  heat 
sources  on  each  chip  has  created  serious  heat  removal 
problems.   There  is  a  genuine  need  for  electronic  designers 
to  not  only  be  aware  of  the  problems  associated  with  higher 
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temperatures,  but  they  must  also  be  prepared  to  solve  those 
problems.   References  [1]  through  [10]  are  a  sampling  of 
the  many  sources  that  indicate  the  increased  emphasis  on 
cooling  problems  throughout  the  electronics  community. 

The  need  to  operate  electronic  equipment  with  maximum 
junction  temperatures  below  those  levels  that  result  in 
failure  is  well  recognized,  but  even  moderately  high  temper- 
atures in  electronic  components  result  in  progressive 
deterioration  and  reduced  reliability.   A  generally  accepted 
thumb  rule  is:   for  every  10  deg  C  increase  in  junction 
temperature,  the  lifetime  of  a  component  will  decrease  by 
one-half.  [3,  7]   The  increasing  costs  associated  with  each 
failure  make  it  imperative  to  address  reliability  during 
the  equipment  design  phase  and  to  provide  sufficient  cooling 
to  maintain  temperatures  as  low  as  practical. 

There  are  many  methods  of  removing  the  heat  dissipated 
within  electronic  components.   These  methods  include  the 
natural  air  convection  previously  mentioned,  thermo-electric 
devices,  heat  pipes,  cold  plates  and  even  complex  refrigera- 
tion systems  which  use  pumped  liquid  coolants  [8].   The 
complex   systems  required  for  these  solutions  can  lead 
to  lower  reliability  through  increased  probability  of 
failures  in  the  cooling  system.   While  many  of  these  methods 
can  allow  extremely  high  power  densities  for  specific  appli- 
cations, the  most  widely  used  method  for  cooling  of  components 
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mounted  on  printed  circuit  boards  is  forced  air  cooling. 
Air  is  readily  available,  abundant,  non-corrosive,  non-toxic, 
non-flammable,  dielectric,  and  is  easily  pumped  using  readily 
available  fans  and  blowers.   For  those  situations  where  there 
is  a  need  to  seal  the  equipment,  interior  cooling  air  is  often 
circulated  through  heat  exchangers. 

All  too  often  the  problem  of  cooling  the  electronic  com- 
ponents is  attacked  after  the  circuit  boards  have  been  designed 
or  even  produced  and  assembled  [9].   In  some  cases,  the  need 
for  increased  cooling  is  recognized  only  after  repeated 
failures  in  service  have  resulted  in  unhappy  users  and  ex- 
cessive repair  costs.   At  best  both  the  electronic  design  and 
the  thermal  design  progress  concurrently  but  are  worked  on  by 
separate  design  groups.   These  groups  may  have  conflicting 
range  goals  that  can  interfere  with  the  important  long  range 
goal  of  reliability. 

Many  of  the  efforts  to  directly  involve  the  electrical 
designer  in  the  thermal  problems  are  directed  towards  overall 
system  cooling.   Reference  [11 ]  details  a  thermal  design 
program  called  VENTBOX.   VENTBOX  treats  a  cabinet  enclosure 
with  entire  circuit  boards  modeled  as  distributed  heat 
sources.   This  program  fails  to  address  the  problems  of  in- 
dividual components  and,  therefore,  is  unable  to  analyze 
the  circuit  board. 

Electronic  circuit  analysis  programs  are  often  utilized 
in  the  thermal  analysis  of  circuit  boards  [12  and  13] .   This 
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technique  requires  the  development  of  an  equivalent  thermal 
circuit,  usually  by  a  packaging  engineer,  and  thus  the  thermal 
design  is  removed  from  the  hands  of  the  electronic  designer. 
Reference  [14]  details  a  thermal  analysis  program  for  circuit 
boards  that  is  much  more  accessible  to  an  electronic  designer 
since  the  inputs  are  physical  dimensions  and  types  of  compo- 
nents rather  than  equivalent  thermal  network  parameters. 

All  these  programs  share  a  major  weakness  in  that  they 
depend  on  a  large  general  purpose  computer  facility.   This 
not  only  can  result  in  excessively  long  turn-around  times 
in  batch  processing,  but  also  the  expense  of  computer  time 
may  become  a  factor.   In  addition,  the  input  data  typically 
consists  of  long  lists  of  numbers  that  must  be  in  the  proper 
format  with  the  correct  option  selection  codes.   Likewise 
the  output  also  consists  of  even  longer  lists  of  numbers 
with  the  key  information  hidden  in  their  midst. 

One  solution  to  automated  thermal  design  and  analysis 
of  avionics  systems  is  being  developed  as  a  joint  project 
by  the  Air  Force  Flight  Dynamics  Laboratory  at  Wright 
Patterson  Air  Force  Base  and  Boeing  Aerospace  Company.   The 
ITAD  (Integrated  Thermal  Avionics  Design)  system  is  expected 
to  include  many  existing  analysis  programs  and  a  large  ever 
growing  data  base  containing  reliability  parameters.   It 
is  anticipated  that  ITAD  will  be  accessed  by  user  through 
both  batch  methods  and  interactively  through  color  graphics 
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terminals.   The  scope  of  this  project  is  enormous  and  the 
benefits  to  the  system  designer  will  be  many,  provided  he 
has  access  to  the  large  computer  at  Wright  Patterson  Air 
Force  Base  [15,  16] . 

Another  solution  on  a  much  smaller  scale  is  to  create 
a  program  for  a  compact,  stand  alone  desktop  computer  that 
is  easily  accessible  to  the  designer  of  electronic  circuit 
boards.   Such  a  system,  if  interactive  and  "friendly"  to 
the  casual  user,  would  allow  rapid  evaluation  of  various 
circuit  board  designs  at  the  conceptual  level.   Rejection 
of  those  designs  with  poor  thermal  characteristics  could 
then  occur  before  they  leave  the  drafting  table  saving  both 
future  efforts  and  dollars. 

B.   OBJECTIVES 

The  main  objective  of  this  thesis  was  to  develop  an 
interactive  thermal  analysis  program  utilizing  the  Hewlett- 
Packard  9845  desktop  computer.   It  was  considered  important 
to  include  the  following  features: 

1.  Simplicity  of  operation:  All  data  input  and  control 
of  program  flow  are  done  in  a  question  and  answer  format  in- 
cluding specific  instructions  with  each  question. 

2.  Graphical  data  checking:   Display  circuit  board 
replicas  on  the  screen  both  for  verification  of  input  and 
to  provide  a  more  useful  form  of  output. 
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3.  Data  correction  capabilities:   Use  interactive  methods 
to  allow  correction  of  portions  of  the  data  without  the  need 
to  repeat  all  the  input. 

4.  Automatic  storage  and  retrival  of  data:   Circuit  board 
descriptions  should  be  written  to  and  read  from  mass  storage 
devices  using  simple  questions  and  answers  rather  than  requir- 
ing specific  knowledge  of  operations  of  the  devices. 

5.  Analysis  methods  hidden:   Various  empirical  heat 
transfer  correlations  should  be  used  based  on  the  physical 
descriptions  and  the  user  should  be  relieved  of  the  need  to 
make  decisions  concerning  the  details  of  the  heat  transfer 
analysis  methods. 

6.  Sensitivity  analysis:   Allow  automatic  parameter 
changes  to  investigate  the  effects  on  the  thermal  performance 
of  these  changes  with  plotted  data. 


18 


II.   DESCRIPTION  OF  THE  THERMELEX  SYSTEM 

A.   GENERAL 

1.   The  Computer 

The  name  THERMELEX  will  be  used  to  refer  to  the 
system  of  interactive  programs  created  to  perform  thermal 
analysis  of  air  cooled  electronic  circuit  boards.   The 
programs  are  written  for  the  Hewlett-Packard  Model  9845 
desktop  computer  using  the  Hewlett-Packard  extended  version 
of  the  basic  language.   A  sophisticated  operating  system 
hides  the  complexity  of  the  9845  from  the  user  and  provides 
protection  from  his  mistakes  yet  provides  the  power  and 
flexibility  needed  for  application  programs  such  as  the 
THERMELEX  system. 

The  computer  used  to  create  and  debug  THERMELEX  is 
the  9845A  with  64  k  bytes  of  read/write  memory  (option  203) , 
internal  printer  (option  500)  and  graphics  package.   In 
addition,  dual  tape  drives  and  dual  floppy  disks  were  utilized 
in  the  development  of  the  system.   These  additional  mass 
storage  devices  provided  considerable  increase  in  the  con- 
venience and  speed  of  storage  and  retrieval  operations  but 
THERMELEX  is  designed  to  operate  with  only  the  standard  tape 
drive.   Recent  advances  in  the  internal  electronics  have 
resulted  in  this  model  being  superseded  by  the  9345B  model 
with  larger  memory  capabilities  and  several  other  features. 

19 


THERMELEX  is  available  in  either  an  A-version  or  B-version 
for  use  in  the  respective  model  of  the  9845. 

Figure  1  shows  the  9845A  but  the  9845B  is  identical 
in  physical  appearance  except  for  the  name  tag.   The  screen 
at  the  top  will  display  the  messages  from  the  system  to  the 
user,  what  is  typed  by  the  user  and  the  output  from  the 
THERMELEX  programs.   This  output  may  be  printed  on  the  screen 
in  what  is  known  as  the  alpha  mode  or  it  may  be  presented  as 
pictures  and  words  from  the  graphics  mode.   Output  of  either 
mode  may  also  be  produced  on  the  thermal  printer  above  the 
keyboard  in  the  inclined  area.    Directly  underneath  the 
screen  there  are  four  pull-out  reference  cards  that  explain 
error  messages  and  other  operating  conditions  that  may  occur. 
At  the  extreme  upper  right  corner  of  the  sloping  section  is 
the  standard  tape  transport  (:T15)  for  mass  storage  of 
programs  and  data. 

The  keyboard  contains  not  only  a  set  of  standard 
typewriter  keys  but  also  several  other  groups  of  keys  that 
are  important  to  THERMELEX.   The  numeric  keys  at  the  lower 
right  allow  easy  entry  of  numbers  as  an  alternative  to 
those  in  the  alphanumeric  group  of  keys  and  allows  numeric 
calculations  to  be  performed  even  during  execution.   Between 
these  groups  are  the  gold  colored  program  control  keys.   All 
responses  to  questions  are  followed  by  pressing  the  CONT  key 
at  the  bottom  of  this  column  of  keys.   The  group  of  special 
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function  keys  in  the  upper  right  are  defineable  by  the  user 
or  from  program  control.   They  are  also  used  in  THERMELEX 
as  priority  interrupts  to  simplify  the  input  of  data.   This 
use  of  these  keys  is  further  explained  later.   The  final  key 
of  interest  to  the  THERMELEX  system  user  is  the  AUTOST  key 
in  the  lower  right  corner  of  the  EDIT/SYSTEM  FUNCTIONS  group. 
This  key  allows  automatic  loading  and  execution  of  the  first 
program  in  the  THERMELEX  system.   Explanation  for  its  use 
is  found  in  Appendix  1,  the  User  Instructions.   For  further 
explanations  of  the  features  of  the  Hewlett-Packard  model 
9845  computer,  see  Ref  [17]  supplied  with  the  computer. 
2.   THERMELEX  System 

Due  to  the  memory  limitations  of  the  9845A,  THERMELEX 
is  divided  into  three  major  programs  which  are  generally 
brought  into  case  from  mass  storage  under  program  control  as 
they  are  needed.   However,  each  program  is  designed  to  stand 
alone  which  can  allow  the  experienced  user  to  bypass  some  of 
the  questions  and  answers  required  to  help  the  inexperienced 
users.   AUTOST  is  the  first  program  of  the  three;  here  several 
pages  of  user  instructions  (see  Appendix  1)  may  be  printed, 
the  special  function  keys  are  defined  and  the  user  is  directed 
along  a  path  to  follow  towards  the  other  programs.   In  BOARDS , 
the  circuit  board  descriptive  data  is  entered,  verified  and 
stored  on  any  one  of  the  mass  storage  devices.   The  user  may 
also  elect  to  produce  a  printed  copy  of  the  input  data  for 
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his  records.   The  THERML  program  performs  the  thermal  analy- 
sis and  provides  output  in  various  tabular  or  graphical  forms 
These  three  major  programs  will  be  described  in  more  detail 
in  the  pages  that  follow. 

There  are  also  several  smaller  files  in  the  system. 
BDSKEY  and  STDKEY  alter  the  definitions  of  the  special 
function  keys.   TNAMES  contains  a  listing  of  all  variable 
names  and  program  section  names  with  explanations  to  aid  in 
any  future  modifications  of  the  THERMELEX  system.   DEMO-P 
and  DEMO-T  are  data  files  containing  example  circuit  board 
descriptive  data  to  allow  demonstrations  of  the  system  and 
to  help  the  new  user  become  familiar  with  the  capabilities. 

B.   AUTOST 

1.   General 

AUTOST  is  the  driver  program  for  the  system.   This 
name  allows  the  program  to  be  automatically  loaded  and  exe- 
cuted if  the  AUTOST   key  is  latched  down  and  if  the  tape 
containing  the  THERMELEX  system  is  in  the  standard  right- 
hand  tape  drive  (:T15)  when  the  main  power  switch  is  turned 
on.   Regardless  of  how  the  program  is  loaded,  the  first 
question  will  be  concerning  the  location  of  THERMELEX.   The 
program  will  ask  which  mass  storage  device  contains  the 
system  since  this  is  necessary  for  correct  program  control. 
At  this  point,  the  program  defines  the  default  mass  storage 
device  using  the  "MASS  STORAGE  IS  " "  command,  and  the 
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default  mass  storage  device  should  not  be  altered  while  using 
THERMELEX. 

2.  Instructions 

AUTOST  will  provide  a  printed  set  of  user  instruc- 
tions either  on  the  screen  in  short  segments  or  on  the  thermal 
printer  in  8.5  by  11  inch  pages.   These  instructions  present 
an  overall  system  view  and  are  intended  to  familiarize  the  in- 
experienced user  and  supplement  the  extensive  instructions 
and  prompts  that  are  presented  on  the  screen  in  all  sections 
of  THERMELEX. 

3.  Special  Function  Keys 

AUTOST  redefines  the  special  function  keys  as  required 
for  the  system  and  can  provide  a  paper  key-code  overlay  as  a 
substitute  for  the  plastic  model  (Hp  Part  #7120-6164) .   Figure 
3  shows  both  styles.   The  paper  style  may  be  used  as  a  guide 
for  filling  in  the  appropriate  blanks  on  the  plastic  version. 

4 .  Option  List 

Finally,  AUTOST  provides  a  list  of  options  for  the 
user  to  select  from.   He  may  load  either  of  the  other  two 
programs  in  THERMELEX  or  produce  a  complete  copy  of  the 
programs  in  the  system  through  selection  of  the  appropriate 
number  from  the  menu  of  options.   This  technique  is  employed 
throughout  the  system  whenever  possible  because  of  the  sim- 
plicity involved  in  entering  a  single  number  over  other 
methods  of  option  selection. 
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5.   The  Copy  Option 

The  copy  option  is  included  as  a  built-in  feature 
since  there  are  many  separate  files  in  THERMELEX  and  a 
separate  command  is  required  for  each  file  to  be  copied. 
Any  mass  storage  device  may  serve  as  the  destination  for 
the  system;  however,  when  copying  to  tape,  it  is  suggested 
that  it  be  blank  due  to  the  number  of  open  blocks  required. 
Following  completion  of  this  option,  program  control 
returns  to  option  list.   It  is  expected  that  the  first 
option  (keyboard  entry  of  circuit  board  description)  is  the 
most  likely  to  be  used  and  the  description  that  follows 
assumes  this  path. 

C .   BOARDS 

1 .   General 

The  major  purpose  of  the  BOARDS  program  in  the 
THERMELEX  system  is  to  provide  a  simple  method  for  input 
of  the  descriptive  data  needed  in  THERML  for  the  creation 
of  the  thermal  model  and  the  analysis.   This  data  includes 
the  geometric  and  material  properties  of  the  circuit  board 
as  well  as  what  components  are  on  the  circuit  board.   These 
components  are  limited  to  DIP'S  (Dual  Inline  Packages)  and 
chip  carriers  (FLAT  Packs) .   The  thermal  model  used  for  the 
analysis  of  the  board  is  a  finite  difference  scheme  and  the 
physical  location  of  the  components  is  approximated  by  the 
locations  of  the  finite  regions  created  by  the  user. 
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The  descriptive  data  may  be  entered  from  the  keyboard 
in  response  to  simple  questions;  or,  it  may  be  retrieved  from 
mass  storage  data  files  previously  recorded  using  BOARDS . 
The  instructions  included  as  Appendix  1  provide  specific 
explanations  of  how  to  use  this  program. 

2 .  Checking  of  Input  Data 

All  input  data  is  displayed  using  either  the  graphics 
capabilities  of  the  9845  or  with  printed  lists  produced  on 
the  screen  or  paper.   The  user  may,  therefore,  easily  verify 
the  validity  of  the  data  he  has  entered.   This  feature  is 
included  whenever  possible  throughout  THERMELEX. 

3 .  Correction  of  Input  Data 

There  are  many  opportunities  to  change  the  descriptive 
data  using  the  same  question  and  answer  format  of  interactive 
programming.   At  no  time  is  it  necessary  for  the  user  to  have 
extensive  knowledge  of  the  program's  internal  workings  or  the 
machine  level  commands  performed  within  the  program.   As  a 
further  means  of  correcting  past  errors,  the  backup  feature 
built  into  the  THERMELEX  system  makes  it  extremely  easy  to 
return  to  previous  questions.   This  feature  is  accessed  by 
depressing  the  special  function  key  (Ko)  in  place  of  a 
response  to  a  request  for  input.   Program  control  will  jump 
back  to  the  previous  question  to  allow  re-entry  of  data. 
Repeated  use  of  Ko  allows  backing  up  to  any  desired  point 
and  resumption  of  program  flow  from  that  point. 
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4.   The  Circuit  Board 

To  describe  the  geometry  of  the  circuit  board  under 
investigation,  the  user  must  enter  the  length,  width,  thickness 
and  the  thermal  conductivity  of  the  circuit  board.   With  no 
input,  the  thermal  conductivity  will  be  assigned  a  default 
value  of  0.2  9  Watts/m-degC ,  a  representative  value  for  the 
bonded  glass  laminates  generally  used  for  construction  of 
circuit  boards.   The  length  parameter  is  defined  to  be  in  the 
direction  of  air  flow  regardless  of  which  dimension  is  the 
largest. 

In  order  to  create  the  finite  difference  thermal 
model,  the  board  must  be  divided  into  logical  regions  by 
placing  any  number  of  equally  spaced  vertical  lines  and  any 
number  of  equally  spaced  horizontal  lines  on  the  board  up  to 
a  total  of  50  regions.   More  regions  would  be  possible  in  the 
9845B  due  to  the  larger  memory  capacity  but  this  would 
require  some  program  changes.   Each  region  will  either  contain 
a  component  assumed  to  be  centered  in  the  region  or  will  be 
empty.   Since  there  may  be  regions  with  no  components,  there 
is  no  unique  set  of  descriptive  data  for  a  given  circuit 
board  design.   The  user  is  encouraged  to  try  different  com- 
binations of  horizontal  or  x  regions  and  vertical  or  y 
regions  to  describe  a  circuit  board.   Figure  4  shows  two 
such  possibilities  for  a  board  with  six  14  pin  DIP's.   Both 
the  six  region  model  and  the  42  region  model  describe  the 
same  board.   Experience  is  helpful  in  making  decisions  as 
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to  how  to  divide  the  board  and  some  boards  may  not  fit  into 
the  constraints  imposed  by  THERMELEX,  but  many  will. 

5.  The  Components 

Each  of  the  defined  regions  will  be  empty  or  contain 
one  of  the  following  components:   14,  16,  24,  40  pin  DIP's 
with  either  vertical  or  horizontal  orientation  and  16,  24, 
40,  64  pin  chip  carriers.   These  components  will  be  assured 
to  be  centered  in  the  region  (with  the  exception  of  the 
40  pin  DIP  which  will  be  assumed  to  occupy  two  regions) . 
The  input  of  component  type  for  each  region  is  via  the 
special  function  keys.   A  special  paper  key-code  overlay  for 
use  in  this  component  input  section  may  be  provided  if  a 
plastic  version  is  not  available,  see  Figure  5.   As  a 
flashing  cross  appears  in  each  of  the  regions,  the  appro- 
priate special  function  key  is  depressed.   This  defines  the 
type  of  component,  length,  width,  orientation  and  draws  the 
component  on  the  screen  for  visual  verification  of  correc- 
tions.  After  all  regions  have  been  defined,  corrections 
are  possible  through  removal  and  replacement  of  components 
using  the  same  special  function  keys. 

6 .  Thermal  Conduction  Paths 

Each  of  the  defined  regions  on  the  board  may  exchange 
energy  with  the  adjacent  regions  via  conduction  through  the 
board  itself  and  any  thermally  conductive  material  added  to 
the  board  such  as  electrical  leads.   In  addition,  many  circuit 
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board  designs  involve  metal  conduction  rails  that  provide 
both  mechanical  support  and  a  method  to  transport  the  heat 
from  the  componets.   There  are  provisions  in  BOARDS  to  model 
either  the  electrical  leads  or  the  conduction  rails.   For 
circuit  boards  without  conduction  rails,  the  user  specifies 
the  average  lead  width  (mm) ,  thickness  and  thermal  conductivity 
of  the  lead  material.   These  geometric  parameters  are  used 
to  determine  the  thermal  resistance  of  a  single  lead  connect- 
ing the  region  centers  or  nodes. 

When  conduction  rails  are  included  on  the  circuit 
board,  the  effects  of  the  electrical  leads  are  ignored.   To 
allow  modeling  of  rails  that  may  have  different  widths,  the 
concept  of  a  conduction  path  unit  (CPU)  is  introduced.   A 
CPU  is  defined  to  be  .1  mm  in  width  but  the  user  specifies 
the  thickness  and  conductivity  of  the  material.   These  para- 
meters are  used  to  determine  the  thermal  resistance  of  a 
single  CPU  connecting  the  region  centers. 

With  the  thermal  resistance  of  single  CPU  or  elec- 
trical lead  determined,  the  user  need  only  specify  the  number 
of  such  resistances  that  connect  the  circuit  board  regions. 
For  example,  if  a  conduction  rail  is  5.8  mm  in  width,  it  may 
be  modeled  as  containing  58  CPUs.    While  the  concept  of 
CPUs  is  totally  artificial,  this  concept  does  allow  modeling 
of  the  Navy  Standard  Electronic  Module  (SEM)  and  Improved 
Standard  Electronic  Module  (ISEM)  for  those  situations  when 
these  modules  are  used  in  forced  air  cooled  systems. 
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With  the  physical  description  of  the  circuit  board 
now  complete,  the  user  may  elect  to  produce  a  picture  of 
the  circuit  board  from  a  dump  of  the  graphics  to  the  internal 
printer.   Pressing  special  function  Key  3  will  produce  a 
picture  such  as  Figure  6.   This  picture  may  be  used  as  a 
final  verification  of  the  validity  of  the  circuit  board 
description  that  has  been  entered.   The  remaining  descrip- 
tive data  to  be  entered  concerns  the  component  parameters. 

7.  Temperature  or  Power  Input 

The  user  must  specify  either  the  average  power  to  be 
dissipated  in  each  component  or  the  maximum  junction  tempera- 
ture.  When  power  levels  are  specified,  the  steady  state 
junction  temperatures  will  be  calculated  in  THERML .   When 
junction  temperatures  are  specified,  the  maximum  allowable 
component  power  levels  will  be  calculated.   In  either  case, 
the  user  must  also  specify  a  case  to  junction  thermal  resis- 
tance for  each  component.   This  is  usually  supplied  by  the 
manufacturer  and  provides  the  means  in  the  thermal  model  to 
link  the  component  power  level  and  junction  temperature  to 
the  case  surface  temperature.   Figure  7  shows  the  data 
summary  sheet  provided  after  the  data  has  been  entered  and 
verified. 

8.  Storage  of  Data 

Although  it  is  possible  to  load  THERML  and  perform 
the  thermal  analysis  directly  since  all  data  is  passed 
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through  a  common  block,  it  is  strongly  urged  that  the  circuit 
board  description  be  stored  on  mass  storage.   A  few  simple 
answers  allows  the  storage  to  tape  or  disk  for  retrieval  at 
a  later  time.   This  data  file  may  be  accessed,  verified, 
altered  and  recorded  back  to  mass  storage  using  BOARDS  to 
investigate  the  effects  that  changes  in  design  have  on  the 
temperatures  or  power  levels  obtained  in  THERML.   Both 
Appendix  1  and  the  program  provide  easy  to  follow  instructions 
for  retrieval  of  previously  recorded  descriptive  data  files. 

Totally  under  program  control,  the  third  major  program 
in  the  THERMELEX  system,  THERML  will  be  loaded  at  the  direction 
of  the  user  and  execution  started  at  the  proper  entry  point. 

D.   THERML 

1.  ■  General 

Recall  from  the  previous  discussion  that  BOARDS  es- 
tablishes all  the  parameters  that  describe  the  circuit  board. 
While  the  major  purpose  of  THERML  is  to  analyze  the  thermal 
performance,  it  must  first  establish  the  parameters  that 
describe  how  the  circuit  board  is  installed  and  cooled,  such 
as  the  spacing  between  boards,  volumetric  air  flow  per  board 
and  inlet  temperature  of  the  cooling  air.   Using  the  total 
set  of  board  and  installation  parameters,  the  thermal  model 
is  set  up,  solved  and  the  results  presented. 

The  approach  taken  in  the  thermal  analysis  is  to 
construct  an  approximate  thermal  network  [7]  for  every  path 
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of  heat  transfer  from  the  components  and  the  circuit  board. 
The  thermal  resistance  (R)  for  each  path  is  calculated  from 
the  set  of  parameters  and  if  a  temperature  difference  (AT) 
for  that  path  is  known,  the  expression: 

will  determine  the  rate  of  heat  transfer  (Q)  for  that  path. 
When  the  rate  of  heat  for  a  given  path  is  known,  the  expression 

AT  =  (Q)  x  (R) 

will  determine  the  temperature  difference. 

An  overall  heat  balance  for  the  components  and  the 
circuit  board,  with  all  the  heat  that  is  generated  going  into 
the  cooling  air  stream,  is  used  to  determine  either  the  steady 
state  junction  temperatures  (when  component  power  levels  are 
specified)  or  the  component  power  levels  (when  the  junction 
temperatures  are  specified) . 

2 .   Air  Flow  -  Thermal  Model 

The  cooling  air  flow  is  assumed  to  come  from  an  in- 
finite heat  sink  which  supplies  air  at  a  constant  inlet 
temperature.   All  the  heat  dissipated  within  the  components 
is  assumed  to  enter  the  air  stream  with  no  heat  conducted 
into  the  card  guides  or  electrical  connectors.   These  guides 
and  the  mechanical  support  sections  of  the  connectors  are 
typically  made  of  plastic  with  high  resistance  to  heat 
transfer  and  this  assumption  leads  to  conservative  estimates. 
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As  the  air  travels  from  inlet  to  outlet,  it  is  assumed 
to  remove  heat  from  both  the  component  surfaces  and  the  surface 
of  the  circuit  board.   It  is  further  assumed  to  travel  in  air 
lanes  defined  by  the  height  of  each  region  perpendicular  to 
air  flow  and  not  mix  until  the  outlet.   As  the  air  removes 
heat  from  each  region,  the  local  air  temperature  (Tair)  will 
increase.   The  process  is  described  by  the  general  equation: 

Tair  at  X  +  AX  =  Tair  at  X  +   Q  added  in  AX 


(Mair)  x  (Cpair) 

where  Mair  =  Mass  flow  rate  of  air  [kg/sec] 

Cpair  =  Heat  capacity  of  air  [J/kg-degk] 
AX  =  length  of  a  region  in  the  air  flow  direction 
This  process  results  in  the  temperature  of  the  cooling  air 
stream  being  modeled  as  a  series  of  steps  as  the  air  travels 
from  inlet  to  outlet  within  an  air  flow  lane.   The  local  air 
temperature  and  the  convective  resistance  determine  the  local 
convective  heat  transfer. 

The  air  velocity  is  determined  from  the  physical 
description  entered  by  the  user  and  this  air  velocity  is 
used  to  determine  the  convective  resistances.   Recall  that 
circuit  board  spacing  (Zb) ,  board  height  (Hb)  and  volumetric 
air  flow  rate  (Fair)  are  part  of  the  descriptive  data  set 
previously  entered.   With  no  components  present,  the  air 
velocity  (Vair)  would  be  defined  by: 

tt  •        Fair 
Vair 


(Zb)  x  (Hb) 
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However,  the  components  tend  to  reduce  the  air  flow  area 
(Aair)  by  blockage.   The  average  blockage  area  is  determined 
by  calculating  the  average  number  of  components  in  a  vertical 
column  of  regions  perpendicular  to  the  air  flow.   In  addition 
the  average  width  (Wavg)  and  height  (Havg)  of  the  components 
is  determined  and  the  average  air  flow  area  calculated  using: 

Aair  =  (Zb)  x  (Hb)  -  (Navg)  x  (Havg)  x  (Wavg) 

This  area  is  used  to  determine  the  air  velocity  from: 
Fair 


Vair  = 


Aair 


In  addition  the  average  area  is  used  to  determine  the  wetted 
perimenter  (Perim)  of  the  air  duct  formed  by  the  boards 
according  to: 

Perim  =  2Hb  +  2Zb  +  2  (Navg)  x  (Havg) 

The  hydraulic  diameter  (Dh)  is  therefore: 
4  Aair 


Dh 


Perim 


which  is  also  used  in  calculations  of  some  of  the  heat, 
transfer  coefficients. 

3.   Components  -  Thermal  Model 

All  heat  generated  within  the  components  must  be 
transferred  away  through  conduction  to  the  board,  radiation 
to  the  board  and  to  some  radiation  sink,  and  through  con- 
vection to  the  cooling  air.   Figure  3  shows  each  of  the 

thermal  networks  for  heat  transfer. 
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DIP  components  are  usually  mounted  by  soldering  the 
electrical  leads  to  the  circuit  board  electrical  conductors 
(Figure  9) .   These  connector  pins  also  act  as  low  resistance 
thermal  conductors  between  the  component  and  the  circuit 
board.   The  total  thermal  resistance  from  the  component  to 
the  circuit  board  may  be  calculated  from  the  cross  sectional 
area  of  each  pin  (Axpin) ,  the  length  of  the  pin  (Lpin) ,  the 
thermal  conductivity  of  the  pin  material  (Kpin)  and  the 
number  of  pins  (Npin)  using  the  expression: 

Rpin  =  LPHL 


(Kpin)  x  (Apin)  x  (Npin) 

There  is  also  a  small  gap  between  the  bottom  of  the 
component  case  and  the  circuit  board.   Figure  9  shows  the 
mounting  for  DIP  cases  where  this  gap  is  on  the  order  of 
one  mm;  for  chip  carriers,  this  distance  is  smaller  yet. 
Due  to  the  small  distance  (Dgap)  involved,  it  is  assumed 
that  no  convection  occurs  in  this  gap  and  that  the  conduc- 
tive resistance  of  the  air  gap  (Rgap-cond)  may  be  calculated 
from  the  expression: 

Rgap-cond  =  (Kair^P  (Agap) 

where  Agap  represents  the  surface  area  of  the  bottom  of  the 

component  and  Kair  represents  the  thermal  conductivity  of 

the  air. 

The  component  case  will  also  radiate  energy  to  the 

circuit  board.   Assuming  the  gap  to  act  as  two  parallel 
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plates  of  equal  areas  with  emissivities  of  Epsb  and  Epse 
the  heat  transfer  rate  may  be  calculated  (Ref  [18]  )  from: 

n  _.  (Sig)  x  (Agap)  x  (Te4  -  Tb4) 
+  -*r=z 1 


where : 


Epsb    Epse 


Sig  =  Stefan-Boltzman  constant 
•8   Watts 


(5.67  x  10 


2        4 
m   -  degK 


Te  =  Component  surface  Temperature  (degK) 

Tb  =  Circuit  board  Temperature  (degK) 

4      4 
This  radiation  term  may  be  simplified  by  expanding  (Te   -  Tb  ) 

in  a  Taylor  series  about  Te  and  retaining  only  the  linear 

portion  of  the  series.   When  this  is  done,  an  effective  gap 

radiation  resistance  (Rgap-rad)  may  be  calculated  from: 

Rgap-rad  =    Epse  +  Epsb  -  (Epse) (Epsb) 

(r)  x  (Sig)  x  (Epse)  x  (Epsb)  x  (Te  ) 

Since  the  component  case  temperature  (Te)  is  an  unknown,  this 

resistance  will  be  recalculated  as  the  solution  progresses. 

These  three  resistances  (Rpin,  Rgap-cond,  Rgap-rad) 
may  be  combined  into  a  total  resistance  between  the  component 
and  the  circuit  board  (Rtote-b) .   When  the  component  is  a 
chip  carrier,  the  air  gap  and  lead  length  are  assumed  to  be 
one-tenth  that  of  the  DIP  case  [13]. 

Although  the  radiation  heat  loss  from  the  component 
is  extremely  small  and  normally  neglected  in  thermal  analysis 
of  electronic  circuit  boards,  it  is  included  here  for  completeness. 
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The  component  will  radiate  to  the  other  components  and  the 
back  of  the  adjacent  circuit  board.   For  the  purposes  of  the 
thermal  model  under  discussion,  it  is  assumed  that  radiation 
is  from  the  top  of  the  component  to  the  adjacent  board  and 
that  this  adjacent  board  temperature  is  the  same  as  the  board 
temperature  below  the  component.   It  is  further  assumed  that 
the  area  of  the  circuit  board  is  much  greater  than  that  of 
the  component.   Using  the  same  linearization  technique  pre- 
viously discussed  results  in  the  expression: 

Rtop-r  =  1    ~    EPSe 


(4)  x  (Sig)  x  (Epse)  x  (Ae)  x  (Te3) 


for  the  effective  thermal  resistance  for  radiation  from  the 
top  of  the  component  to  the  adjacent  circuit  board. 

Convection  heat  transfer  to  the  air,  the  final  heat 
loss  path  from  the  component,  depends  on  the  heat  transfer 
coefficient  and  the  area  for  that  mode  of  heat  transfer. 
The  component  is  assumed  to  experience  two  separate  modes  of 
convection.   Those  portions  of  the  component  that  are  perpen- 
dicular to  the  air  stream  are  assumed  to  experience  a  stagna- 
tion form  of  air  flow  while  the  top  and  sides  see  parallel 
air  flow.   Recall  that  air  flow  is  defined  as  being  from  left 
to  right  as  required  for  data  entry  in  the  BOARDS  program. 
When  component  types  are  entered  as  being  horizontal  or 
vertical  using  the  special  function  keys,  the  component 
length  (Le)  and  width  (We)  parameters  are  set  such  that  Le 
is  along  the  air  flow  direction  and  We  is  perpendicular  to 
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the  air  flow.   The  standard  correlation  for  plane  stagnation 
flow  may  be  written  as:   Ref  [19] 

Hstag  =  (.57)  x  (Kair)  x  (Pr*4)  x  (■ 


(We)  x  (Gnu) ' 

where  Hstag  =  Heat  transfer  coefficient  for  areas  receiving 
stagnation  air  flow 

Pr  =  Prandtl  number  for  air 

Gnu  =  Kinematic  viscosity  for  air 
Using  the  component  height  (He) ,  the  area  for  stagnation  heat 
flow  (Astag)  may  be  calculated  from: 

Astag  =  (2)  x  (We)  x  (He) 

The  portions  of  the  component  that  are  parallel  to 
the  air  flow  (top  and  sides)  are  assumed  to  experience  the 
same  heat  transfer  coefficient  that  would  occur  in  a  smooth 
duct  modified  by  a  roughness  factor  (Fr)  to  account  for  the 
presence  of  the  components  acting  to  increase  this  heat 
transfer  coefficient  (Kpara) .   Hpara  depends  on  the  Reynolds 
number  (Re)  determined  by: 

=  (Vair)  x  (Dh) 
Gnu 

In  the  entrance  region  of  the  duct  formed  by  the  circuit 
boards,  up  to  a  distance  of  approximately  ten  hydraulic 
diameters,  the  flow  is  assumed  to  be  laminar.   Reference  [20] 
contains  the  following  correlation  for  the  heat  transfer 
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coefficient  in  the  entrance  region  of  a  smooth  duct: 

r  <** 

Mllf,t  -  (.664)  x  (Kair)  f  (Gz)  x  (l+(7.3)  x  lGz 
Ct  "  (1.1)  x  (Dh)     I   Pr" 


where  Gz  is  the  Graetz  number  defined  as: 

Gz  _  (Re)  x  (Pr)  x  (Dh) 
Xi 

and  Xi  is  the  distance  from  the  entrance  of  the  cooling  air 
flow. 

The  value  of  the  heat  transfer  coefficient  predicted 
from  the  above  equation  has  been  found  to  be  low  from  compari- 
sons to  experimental  data  of  Ref  [12]  and  Ref  [14] .   The 
expression: 

Fr  =  1    +  5  x  (Perim  -  (2)  x  (Bdh)  -  (2)  x  (Zb)  ) 

Perim 

has  been  created  to  adjust  the  predicted  heat  transfer  co- 
efficient for  the  rough  duct.   The  expression  reduces  to  one 
when  no  components  are  present  and  is  less  than  two  for  all 
reasonable  board  constructions.   The  resulting  heat  transfer 
coefficient  that  is  used  for  the  parallel  sides  (Hpara)  is 
therefore: 

Hpara  =  (Hduct)  x  (Fr) 

For  positions  beyond  the  entrance  length  the  flow 
may  be  laminar  or  turbulent  and  the  appropriate  heat  transfer 
correlation  must  be  chosen.   The  transition  from  laminar  to 
turbulent  is  assumed  to  occur  at  a  Reynolds  number  of  1000 
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since  the  components  act  as  turbulence  promoters  for  the 
air  flow.   For  laminar  flow,  the  expression: 

Hduct  =  (5-4)  jj  'Kair> 
Dh 

is  used  and  the  result  modified  with  the  roughness  factor 
(Fr)  to  obtain  Hpara.   For  turbulent  flow  the  Dittus-Boelter 
relationship  is  used: 

Hnar3  _  (.023)  x  (Re,S)  x  (Pr*4)  x  (Kair) 
Hpara ^ 

Heat  is  convected  not  only  off  the  top  and  sides  of 
the  components  but  also  from  the  surface  of  the  leads  or 
pins.   Since  the  thermal  conductivity  of  the  pin  material 
is  so  high,  their  surface  is  assumed  to  have  the  same  temper- 
atures as  the  surface  of  the  component.   The  effective  com- 
ponent parallel  side  area  for  heat  transfer  is,  therefore, 
determined  from: 

Apara  =  (Le)  x  (We)  +  2 (Le)  x  (He)  +  (Aspin)  x  (Npin) 

where  Aspin  is  the  surface  area  of  the  pins  that  experience 
parallel  air  flow.   The  thermal  resistances  for  convection 
from  the  component  may  then  be  calculated  from: 

RPara  =  (Hpara)  x  (Apara)        and' 

1 

Rstag    (Hstag)  x  (Astag) 
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These  two  resistances  may  then  be  combined  as  parallel  re- 
sistances to  give  a  total  convective  thermal  resistance  from 
the  component  to  the  air  (Re-conv)  of 

Re-conv  =  (Rpara)  x  (Rstag) 
(Rpara)  +  (Rstag) 

Returning  now  to  the  basic  premise  that  all  the  heat 
produced  within  the  component  must  be  transferred  through  one 
of  the  heat  paths  illustrated  in  Figure  8,  the  heat  balance 
equation  for  the  Ith  component  is: 

Pow  (I)  =  Qto  board  +  Qto  air  +  Qrad  off  top 

utilizing  the  general  expression  for  heat  flow  as  a  function 
of  temperature  difference.   This  may  be  expressed  as: 

Pnw    m    •=   Te(I)    -   Tb(I).       Te(I)    -   Tair(I) 
UJ    '        RToTe-b(I)  Re-conv (I) 

Te(I)    -   Tb(I) 


Rtop    -    R(I) 


This  expression  may  be  rearranged  and  solved  for  the  compo- 
nent temperature  (Te) .    The  resulting  expression  is  used  to 
eliminate  Te  from  the  final  set  of  equations  that  are  solved 
for  the  set  of  board  temperatures. 

4 .   The  Circuit  Board  -  Thermal  Model 

The  circuit  board  has  been  divided  into  regions  as 
determined  by  the  user  while  entering  the  circuit  board 
description  in  program  BOARDS.   Each  of  the  regions  will 
lose  or  gain  heat  by  the  normal  processes  of  conduction, 
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convection  and  radiation.   In  this  thermal  model,  radiation 
from  the  circuit  board  regions  is  neglected  due  to  the  lower 
temperatures  of  the  circuit  boards.   The  board  temperature 
(Tb)  is  assumed  to  be  uniform  within  a  given  region  for 
purposes  of  convective  heat  transfer  to  the  cooling  air 
stream.   For  purposes  of  calculating  the  heat  conducted 
between  regions,  the  temperature  differences  are  assumed  to 
exist  over  the  distance  between  the  centers  of  the  regions. 
These  region  centers  are  the  nodes  in  the  finite  difference 
model  used  to  analyze  the  circuit  board. 

The  surface  of  the  board  forms  part  of  the  air  duct 
previously  mentioned  in  the  discussion  concerning  the  calcu- 
lation of  the  convective  heat  transfer  coefficient  from  the 
parallel  sides  of  the  components  (Hpara) .   The  heat  transfer 
coefficient  of  the  board  varies  with  distance  from  the  air 
centered  in  the  regions,  the  convective  heat  transfer  co- 
efficient for  the  board  (Hb)  is  assumed  to  be  equal  to  Hpara. 
The  area  for  heat  transfer  from  each  region  (Areg)  includes 
both  sides  but  excludes  any  area  under  the  component  (Ae) 
since  the  air  is  assumed  not  to  flow  in  the  small  gap  between 
the  component  and  the  circuit  board.   The  area  (Areg)  is  cal- 
culated from  the  equation: 

=  (2)  x  (Bdl)  x  (Bdh)  .  Ae 
y  Nreg 

where  Bdl  =  Length  of  circuit  board 

Bdh  =  Height  of  circuit  board 
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The  convective  resistance  from  the  circuit  board  region 
(Rb-conv)  is  therefore: 


Rb-conv 


(Hb)  x  (Areg) 

The  conduction  of  heat  within  the  circuit  board 
material  and  any  added  conductive  material  is  assumed  to 
occur  between  the  nodes  located  at  the  centers  of  each 
region.   For  a  given  region,  this  heat  flow  is  assumed  to 
occur  only  out  the  four  sides  of  the  region.   The  distance 
between  nodes  in  the  horizontal  or  air  flow  direction  (Lr) 
is  calculated  from  the  circuit  board  length  (Bdl)  and  the 
number  of  regions  in  the  horizontal  direction  (Nxr) .   The 
expression  used  is: 

Lr  -  Bdl 
Lr    Nxr 

Likewise,  the  vertical  spacing  between  nodes  (Hr)  is  calcu- 
lated as: 

Bdh 
Hr  =  NyF 

These  distances  are  used  with  the  circuit  board  thickness 
(Thb)  and  the  circuit  board  material  thermal  conductivity 
(kb)  to  determine  the  base  board  conductive  resistances  in 
both  the  vertical  (Rb-ver)  and  horizontal  direction  (Rb-hor) . 
The  expressions  used  are: 

Hr 
Rb-ver  = 


(kb)  x  (Lr)  x  (Thb) 
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and 


Rb-hor  - 


(kb)  x  (Hr)  x  (Thb) 


The  conductive  resistances  in  the  material  added  to 
the  circuit  board  is  treated  in  the  same  manner  since  these 
materials  are  also  assumed  to  connect  the  nodes.   The  con- 
ductivity of  the  material  (kl)  and  physical  area  (al)  have 
been  specified  for  either  a  conduction  path  unit  (CPU)  or 
electrical  lead.   These  are  used  in  the  expressions: 

Rl-ver  =     Hr 


Rl-hor  = 


(kl)  x  (Al) 
and 
Lr 


(kl)  x  (Al) 

to  determine  the  thermal  resistance  for  conduction  in  a 
single  conduction  path.   This  resistance  is  divided  by  the 
number  of  such  paths  to  determine  the  total  conductive  re- 
sistance of  the  material  added  to  the  circuit  board. 

This  resistance  is  assumed  to  be  in  parallel  with 
the  conductive  thermal  resistance  of  the  bare  circuit  board 
and  a  net  conductive  thermal  resistance  is  calculated  for 
each  of  the  four  directions  using  a  normal  product-over-sum 
formula  for  parallel  resistances.   Figure  10  shows  how  a 
typical  region  would  thus  be  connected  to  the  four  adjacent 
regions.   The  values  of  these  conductive  thermal  resistances 
are  stored  in  a  two  dimensional  array,  Rl  (I, J),  where  the 
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first  index  represents  the  region  number  and  the  second  index 
is  the  direction  number  (1,  2,  3,  4).   These  numbers  represent 
bottom,  right  side,  top,  and  left  side,  respectively.    For 
example,  Figure  10  shows  how  Rl  (1,2)  connects  region  number  I 
and  region  number  I  +  1.   For  those  regions  on  the  edges  of 
the  circuit  board,  the  conductive  thermal  resistance  in  the 
direction  off  the  board  are  set  to  very  high  value  due  to  the 
assumed  adiabatic  boundary  conditions  for  all  edges.   The  set 
of  conductive  thermal  resistances,  Rl  (I, J) .  is  used  in  the 
heat  balance  for  the  circuit  board. 

With  the  component  treated  as  the  source  of  heat  for 
the  circuit  board,  one  may  again  apply  a  simple  heat  balance 
for  a  region  as  follows: 

Qe-b  =  Qconv  +  Q,  +  Q2  +  Q   +  Q 

where  Qconv  represents  the  heat  transfer  to  the  cooling  air 
and  Q,  through  Q,  represent  the  heat  conducted  to  the  adja- 
cent regions  as  shown  in  Figure  10.  This  expression  may  be 
written  for  Region  I  as: 

Te(I)  -  Tb(I)  _  Tb(I)  -  Tair(I)    Tb(I)  -  Tb(I  +  Nxr) 
Rtota-b         Rb-conv  Rl(I,i) 

L   Tb(I)  -  Tb(I  +  1)  ^  Tb(I)  -  Tb(I-Nxr) 


(RKI,2)  RKI, 3) 


,  Tb(I)  -  Tb(I-j) 
R1(I, 4) 


As  previously  stated,  the  heat  balance  equation  for  the  component 
derived  in  Section  3  may  be  solved  for  Ted)  and  that  expression 
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used  to  eliminate  Te(I)  from  the  equation  above.   The  only 
unknowns  will  then  be  the  board  temperatures.   A  heat  balance 
for  every  region  may  be  performed  resulting  in  Nreg  simultan- 
eous algebraic  equations.   The  coefficients  of  these  equations 
may  then  be  collected  into  a  matrix  and  solved  using  an  LU 
decomposition  method  [21]. 
5.   The  Analysis 

In  performing  the  thermal  analysis  of  the  circuit 
board,  there  are  two  separate  situations  to  be  considered: 
1.   Specified  component  power  -  determine  steady  state 
junction  temperature;  2.   Specified  junction  temperature  - 
determine  the  maximum  power  allowable.   For  the  first  situa- 
tion a  case  temperature  (Te)  is  assumed  for  each  component 
since  this  is  necessary  to  determine  the  effective  radiation 
resistances.   All  thermal  resistances  and  the  local  air 
temperatures  are  then  calculated  using  the  given  variables. 
The  resulting  system  of  Nreg  simultaneous  equations  is  solved 
using  LU  decomposition  for  the  temperature  of  each  region  of 
the  circuit  board  Tb. 

This  set  of  board  region  temperatures  is  used  to 
determine  the  set  of  component  case  temperatures  using  the 
expression  for  Te  determined  from  the  heat  balance  of  the 
component.   These  component  temperatures  are  compared  to 
those  from  the  previous  iteration.   If  the  largest  difference 
between  any  of  the  temperatures  is  less  than  the  convergence 
criterion  established  by  the  user  (default  .1  degC),  the 
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results  are  presented  in  the  graphics  mode  on  a  fascimile 
of  the  circuit  board  as  shown  in  Figure  11. 

In  the  display  of  the  circuit  board,  each  region 
contains  the  component  type,  junction  temperature,  power 
level  and  case  temperature  for  the  component  in  that  region. 
Empty  regions  contain  only  the  temperature  of  the  circuit 
board.   Those  junction  temperatures  that  are  within  5%  of 
the  maximum  are  starred  (**)  for  easy  reference.   This  output 
is  dumped  to  the  internal  thermal  printer  to  provide  a  hard 
copy. 

If  convergence  has  not  been  reached,  the  new  case 
temperature  is  used  to  calculate  new  radiation  thermal  re- 
sistances and  the  new  system  of  equations  solved.   Figure  12 
shows  the  intermediate  display  of  all  temperatures  that  is 
presented  on  the  screen  while  the  next  iteration  is  in 
progress.   If  longer  than  20  lines,  the  maximum  display  area 
for  output  on  the  screen,  this  display  may  be  moved  up  or 
down  with  the  display  control  keys  at  the  center  top  of  the 
keyboard  (see  Figure  2) .   Convergence  typically  occurs  in 
less  than  four  iterations  when  solving  for  component  junction 
temperatures . 

For  the  second  situation  with  specified  junction 
temperatures  (T j ) ,  a  power  level  of  .25  watts  is  assumed  for 
each  component.   The  component  case  temperature  may  then  be 
determined  from  the  expression: 

Pow 
Te  =  Tj    - 


Rj-c 
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where  Rj-c  represents  the  case  to  junction  thermal  resis- 
tance set  by  the  user.   As  in  the  first  situation,  all  the 
thermal  resistances  and  the  local  air  temperatures  (Tair) 
are  calculated.   It  is  important  to  note  that  Tair  depends 
on  the  assumed  power  levels  since  these  values  will  change 
with  each  iteration.   The  resulting  system  of  simultaneous 
equations  is  again  solved  for  the  board  region  temperatures 
and  the  component  case  temperatures.   The  resulting  compo- 
nent power  levels  for  each  region  are  determined  from  the 
expression 

Te  -  Tb 
Pow  - 


Rtote-b 


This  component  power  level  for  each  component  is  compared  to 
the  results  of  the  previous  iteration  or  the  assumed  values 
for  the  first  iteration.   Convergence  is  assumed  when  the 
largest  percentage  difference  from  these  comparisons  is 
less  than  the  convergence  criterion  established  by  the  user 
(Default  1%) . 

The  output  is  presented  in  the  same  manner  as  before 
except  those  power  levels  within  5%  of  the  minimum  are 
starred  (**)  in  this  situation  (Fig.  13) .   Since  both  the 
radiation  thermal  resistances  and  local  air  temperatures 
depend  on  the  power  levels,  more  iterations  are  required 
before  convergence  is  reached.   Typically  less  than  six  are 
sufficient. 
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6.  What  Now  Option  List 

Following  a  complete  cycle  through  THERML,  there  are 
several  options  available  to  the  user.   A  different  circuit 
board  may  be  analyzed  either  by  entering  the  descriptive  data 
set  from  the  keyboard  or  by  retrieval  from  mass  storage.   The 
previously  recorded  data  set  may  be  read  in  using  THERML  with 
no  data  checking  capabilities  or  BOARDS  may  be  loaded  under 
program  control  to  allow  data  checking  or  changes  to  the 
circuit  board. 

In  addition  the  same  circuit  board  may  be  reanalyzed 
with  a  new  set  of  installation  parameters,  i.e.,  board  spacing, 
inlet  air  temperature  and  volumetric  air  flow  rate  per  board. 

7.  Sensitivity  Analysis 

Another  option  available  allows  the  user  to  investi- 
gate the  effects  of  various  air  flow  rates.   The  user  specifies 
a  maximum  air  flow  rate  and  five  separate  analyses  are  performed 
for  each  of  five  air  flow  rates  up  to  the  maximum  specified. 
Figure  14  shows  how  the  results  of  this  sensitivity  analysis 
are  presented  for  the  situation  of  specified  power  levels 
while  Figure  15  is  an  output  for  the  second  situation  of 
specified  junction  temperatures.   The  outputs  from  each 
analysis  such  as  Figures  11  and  13  are  not  presented  during 
the  sensitivity  analysis  since  the  graphics  mode  is  used  but 
an  output  such  as  Figure  12  may  be  presented  on  the  thermal 
printer  if  records  of  individual  components  are  needed  for 
each  of  the  separate  air  flow  rates. 
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A  final  option  available  to  the  user  is  termination. 
It  is  important  to  terminate  the  THERMELEX  system  under 
program  control  to  insure  normal  key  definitions  are  returned 
and  graphics  parameters  correctly  assigned. 
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HI.   RESULTS  AND  CONCLUSIONS 

The  THERMELEX  system  offers  the  designer  of  electronic 
circuit  boards  the  means  to  predict  the  thermal  performance 
of  air  cooled  circuit  boards  and  avoid  the  problems  that 
often  surface  only  after  the  equipment  is  in  use.   The 
system  is  easy  to  use  yet  provides  valuable  data  in  a 
variety  of  formats  that  can  help  the  designer  to  make  im- 
portant design  decisions  regarding  circuit  board  layout  and/or 
cooling  air  flow  parameters. 

The  THERMELEX  system  has  been  tested  with  a  wide  variety 
of  circuit  boards  to  insure  that  various  combinations  of  the 
possible  components  and  empty  regions  will  create  valid 
descriptive  data  sets  and  reasonable  results.   In  all  test 
cases  the  results  have  been  satisfactory  when  compared  to 
expected  results.   For  example,  higher  component  power 
levels  result  in  higher  junction  temperatures  and  decreasing 
air  flow  produces  higher  temperatures  with  all  else  the  same. 
Several  tests  were  made  that  used  the  results  of  a  power 
level  prediction  as  input  to  the  same  circuit  board  to 
insure  that  predicted  junction  temperatures  were  the  same 
as  those  specified  for  the  original  test. 

In  addition  to  the  above  tests  for  general  validity, 
direct  comparisons  to  experimentally  determined  component 
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case  temperatures  reported  in  reference  (14)  were  made. 
The  circuit  board  consisted  of  25  equally  spaced  14  pin  DIPs 
and  is  depicted  in  Figure  6.   Three  different  air  flow  rates 
and  four  different  component  power  levels  were  used.   The 
results  of  these  experiments  are  shown  in  Figure  16.   These 
same  flow  rates  and  power  levels  were  used  as  inputs  for 
THERMELEX  and  the  predicted  maximum  case  temperatures  are 
also  included  in  Figure  16.   At  the  lower  power  levels,  the 
agreement  with  the  experimental  results  is  encouraging  but 
for  the  higher  power  levels  there  are  significant  differences, 
In  particular,  the  predicted  maximum  case  surface  temperature 
shows  a  much  stronger  dependence  on  the  air  flow  rate  than 
the  experimental  results  would  indicate.   It  is  believed 
that  several  effects  neglected  in  the  the  thermal  model 
become  significant  for  low  flow  rates  and  higher  component 
power  levels.   For  low  air  flow  rates,  the  effects  of  natural 
convection  heat  removal  become  more  important,  thus  decreas- 
ing the  actual  surface  temperature.   In  addition,  at  higher 
temperatures,  the  conduction  into  the  electrical  connectors 
and  mechanical  supports  will  also  tend  to  hold  the  surface 
temperatures  lower  for  the  experimental  results.   Further 
work  is  needed  to  resolve  the  differences  between  the 
experimental  results  and  those  predicted  by  THERMELEX. 
Particular  attention  towards  refinement  of  the  thermal  model 
is  required. 
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IV.   RECOMMENDATIONS  FOR  FUTURE  WORK 

While  the  present  version  of  THERMELEX  can  be  a  valuable 
tool  for  predicting  thermal  performance  of  electronic  circuit 
boards,  improvements  and  refinements  would  be  useful  in 
several  areas.   The  first  recommendation  would  be  for  more 
experimental  verification  with  particular  emphasis  towards 
developing  better  empirical  heat  transfer  relationships  for 
use  in  the  present  thermal  model.   The  thermal  model  should 
also  be  expanded  to  include  both  the  effects  of  natural  con- 
vection and  boundary  conditions  other  than  adiabatic.   In 
particular  the  user  should  be  able  to  specify  a  constant 
temperature  for  one  or  more  of  the  physical  circuit  board 
boundaries.   This  would  allow  the  modeling  of  installations 
that  include  metal  card  guides  or  cold  plates. 

The  present  model  does  not  include  the  interaction  that 
can  occur  between  circuit  boards  mounted  closely  together. 
These  effects  could  be  included  in  the  radiation  sink 
temperatures  "seen"  by  a  board  or  included  as  effects  on  the 
local  air  temperature  for  the  bottom  and  top  of  the  circuit 
board. 

Additional  improvements  could  also  be  made  in  order  to 
increase  the  number  of  circuit  board  designs  that  may  be 
analyzed  with  THERMELEX.   The  limited  set  of  components 
could  be  increased  to  include  DIPs  with  different  numbers 
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of  pins  and  other  flat  pack  case  styles.   Discrete  components 
such  as  resistors,  capacitors  and  switches  are  also  able  to 
transfer  heat  to  the  air  stream  and  could  also  be  included. 
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Figure  4.   Alternate  Region  Definition  for  Single  Circuit  Board 
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THE  ABOVE  DATA  13  FOR  NELONPGS  TE3T  BOARD 
BOARD  LENGTH  (defined  along  air-  flow>=  142  mm        HEIGHT  =  114  mm 
BOARD  THICKNESS*  1.443  mm  CONDUCTIVITY  =  .2342  Watts/M-K 
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Figure  7.   Data  List  of  Circuit  Board 
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Figure  9:  DIP  Mounting  on  Circuit  Board 
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APPENDIX  B 


10 
20 
30 
48 
50 
60 
?0 
30 
90 
100 

1  10 
120 
130 
140 
150 
160 
170 
130 
190 
200 
210 
220 

2  30 
240 
250 
260 
270 
230 


AUTOST  DRIVER  PROGRAM  FOP  

**************  THERMELEX         ****** ********* 

R  SYSTEM  OF  PROGRAMS  FOR  THE  HP  9*45 

PERFORM  THERMAL  ANALYSIS  OF  ELECTRONIC  CIRCUIT  BOARDS 

PREPARED  AT  THE  NAVAL  POSTGRADUATE  SCHOOL  MONTEREY 
R.  A.  FOLTZ  LCDR  USN 

I  I  I  I  I    RUTOST   |  |  |  |  | 


GOTO  240 

ABOUT  SPECIAL  FUNCTION  KE' 


>r-ke' 


PRINTER  IS  0 

PRINT  CHR*<27  'S  "8,1  34T" 

DIM  S*C60] 

COM  Map 

PRINTER  IS  16 

GOSUB  Set sto 

GOSUB  I  nit 

IF  finsf^BACKJJP"  THEN 

GOTO  Start 
290  St  in-::   GOSUB  Error 
300  Start:   3*="N0TE  TO  USER 
310   ON  ERROR  GOTO  3t  art  -: 
320    GOSUB  Pagehead 
330    PRINT  "     The  Pol  lowing 
d.  " 
340    PRINT  "Ir,  general  these  key  deffinitic 

ii  waiting" 
350    PRINT  "for  in  i  npu 
_up  opt  i  on" 

360    PRINT  "<k0>  return 
370    PRINT 
330    PRINT  "This  should  be  used  when  the 

i  nput  t  h at " 
390    PRINT  "can  be  corrs 
400    PRINT 

4  10    PRINT  "If  you  do  not  have  a  plastic 
6164)" 

420    PRINT  "a  paper  one  will  be  pro 
430    PRINT  LIN'  1  • 
440    PRINT  "      KEY  'k0' 
450    PRINT  "      KEY  'kl' 
460    PRINT  "      KEY  'k2' 
470    PRINT  "      KEY  'k3' 
480    PRINT  "      KEY  'k4' 
490    DISP  "Freis  CONT  WHEN  RE 
500    PAUSE 

510  PRINT  "  KEY  'k5'  will  enable  1 nt 
520  PRINT  "  KEY  'K6  ,,M  l  l  rewind  let 
530    PRINT  "      KEY  'K7'  wi 1 1 


1  b. 


1  ■=  i 


5  e  n  o  t  e 


id 


i e  Run  1 i ght  i n  1 ower 
he  pr  o q r  am  t  o  a  p 


ed  b 


Ght.   ,; 

e  y  i  o  u  s  q  u  e  i 
r  decides  h 

the  data." 
c  o  d  e  o  >■>  e  r  1 


r  am 
Ba< 


ded 


1  back  up  to  the  previou 
1  d  i  sp  1  ay  p  r  i  n  X  e  d  o  u  t  p  u  t 
1  d  i  s  p  1  -x\>  plotted  o  u  t  p  u  t 
1  print  copv  of  plotted 
I  d  i  sab 1 e  internal  p r i  n t  > 
ADY  FOR  MOPE  EXPLANATION: 


1    printei 

-,d    tape, 
d    right     hand    tape 


•i  1  ■.'    when    the    pro 
l  e  r    i  s    o  f  f  • .       T  h 
:-n    or    control     point 
■i  ■•  s    m  ad  e    i  o  m  e    e  r  r  o  r 

f  i  1  led  l n  (hp  tt  7: 


i  on.  " 

(EXIT  GRAPHICS 
(GRAPHICS)" 

(DUMP  GRAPHICS 
(PRINTER  IS  16 


(PRINTER  IS  0 
kPEWIND  IT14 
(REWIND   :T15 
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AUTOST  ..  DRIVER  FOR  THEPMSLEX 

540    PRINT  "      KEY  "K14'will  answer  Yes. " 

550    PRINT  "      KEY  '  K  1  5  '  u i  11  answer  No." 

560    PRINT  LIN<2> 

570    PRINT  "In  general  when  answering  questions  or  selecting  an  option  pressing 

CONT   " 
5S0    PRINT  "with  no  other  entry  wi  1  1  assign  either  the  default  or  the  first  opt- 
ion in" 

590    PRINT  "parentheses  (the  first  option  in  the  list?  as  the  desired  option." 
600    PRINT  LINO) 
610    PRINTER  IS  16 
620    GOSUB  Over! ay 
630    IF  Ans*="BACK  OP"  THEN  260 
640    DISP  "LOADING  KEY  DEFF INI T IONS" 
650    IF  Map  =  0  THEN  LORD  KEY  "BBSKEY" 
660    WAIT  1000 
670  GOTO  S90 

630  Options*:   GOSUB  Error 
690  Options:   S*="MflIN  PROGRAM  OPTIONS" 
700   ON  ERROR  GOTO  Optionse 
710    GOSUB  Pagehead 
720    PRINT  "     The  THERMELEX  program  package  consists  of  three  major  sections, 

one  of" 
730    PRINT  "which  is  in  core  now  providing  these  instructions.   The  other  two  a 
re  called" 

740    PRINT  "SCARPS  and  THE- ML.   Both  have  the  capability  to  read  circuit  bo 
ard  da'  a" 
750    PRINT  "from  a  Mass  Storage  Device  <eg  Tape  or  Disk)  but  only  BOARDS  can: 

760  PRINT  "(a)  Input  circuit  beard  description  from  keyboard..."; 

770  PRINT  "(b)  Allow  graphical  checkingot  descriptions.."; 

730  PRINT  "■,  c)  Hake  corrections  to  descriptions  as  required..."; 

790  PRINT  »<d>  Record   descriptions  to  mass  storage..." 

SCO  PRINT  LINO),  "THERML  performs  the  thermal  analysis  and  output,  of  results 

810    PRINT  LIN<1>, "YOUR  OPTIONS  APE:"; 

820    PRINT  "    1.  LOAD  BOARDS  TO  INPUT  DESCRIPTION  FROM  KEYBOARD." 

830    PRINT  LINU),SPfl<20>,"2.  LORD  BOARDS  TO  READ  DESCRIPTION  OFF  MASS  STORAG 

E" 

840    PRINT  LINU),SPfl<20),  "3.  LOAD  THERML  TO  READ  DESCRIPTIONS  OFF  MASS  STOF'A 

GE"!  WITH  IMMEDIATE  ANALYSIS." 

350    PRINT  LIN<1>,SPA<20>, "4.  COPY  ALL  PROGRAMS  " 

860    Ans*="l" 

870    INPUT  "YOUR  CHOICE  ?  < 1 , 2, 3, 4> " , fins* 

830    IF  Ans*="BACKJJP"  THEN  Start 

890    Ans=INT< VRL> flnsl) • 

900    IF  (flns>0)  AND  Cflns  3)  THEN  930 

910    GOSUB  Err  in 

920    GOTO  Options 

930    IF  Ans=l  THEN  Map  =  2 

940    IF  flns=2  THEN  Map*3 

950    ON  *nz    GOSUB  Gen, Gen, Therml , Copy 

960    GOTO  Opt i ons 

970  Gen:   DISP  "LOADING  GENERAL  CIRCUIT  BOARD  PROGRAM  " 

9S0   LOAD  "BOARDS" , 1 

990   STOP 

1000  Therml: DISP  "WORKING  LOADING  THERML  " 

1010  LOAD  "THERML", 10 

1020  STOP 
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HUTOST  ..  DRIVER  FOR  THEFNELEX 

1030  Copye: GOSUB  Error 

1040  Cop'/:  '  THIS  SECTION  FOR  PROGRAM  REPRODUCTION  ONTO  ANOTHER  MASS  STORAGE 

1050  ON  ERROR  GOTO  Copye 

10-50  S*="  COPYING  THERMELEX  PROGRAMS" 

1O70  GOSUB  Pagehead 

1O30  PRINT  "    This  section  allows  easv  reproduction  of  THERMELEX  from  one  mass 

storage" 

1090  PRINT  "device  to  another.   The  mass  storage  medium  <taps  or  disk)  must  be  i 

nitiali zed. " 

1100  PRINT  "It  the  medium  intended  as  the  destination  is  new  or  you  wish  to  eras 

ae  al 1  " 

1110  PRINT  "files  then  the  use  of  option  2  will  perform  this  task.   Since  the  TH 

ERMELEX" 

1120  PRINT  "='  stem  requires  only  approximatly  500  records, it  is  possible  to  use 

a  medium" 

1130    PRINT    "that    already    contains    files    to    be    retained;    however,    this    may    not    al 

low" 

1140  PRINT  "sufficient  space  for  the  data  files  containing  the  circuit  board  des 

c  r  i  pt  i  ons . " 

1150  PRINT  "Option  1  will  record  in  the  available  spaces  if  possible  but  it  is  r 

•commended" 

1160  PRINT  "that  an  entire  Mass  storage  media  be  devoted  to  THERMELEX." 

1170  PRINT  LINO),  "YOUR  OPTIONS  APE:"; 

1130  PRINT  "    1.  COPY  ALL  PROGRAMS  WITHOUT  INITIALIZATION." 

1190  PRINT  LIN<l),SPfi<20), "2.  COPY  ALL  PROGRAMS  AFTER  INITIALIZATION." 

1200  flns*="l" 

1210  INPUT  "YOUR  CHOICE  FROM  ABOVE  (1  or  2'1",An  =  * 

1220  IF  Rns$«"BflCK  UP"  Then  Options 

1230  Ans=INT(VAL'  Ans J  '  > 

1240  IF  Cflns=l)  OP  (flns«2)  THEN  1270 

125*  GOSUB  Err in 

1260  GOTO  Copy 

1270  Rns*=Msusl* 

12S0  EDIT  "WHAT  MASS  STORAGE  DEVICE  CONTAINS  THE  THERMELEX;  SYSTEM  •.  :  T  1  5  ,  :  T  1  4  ,  :  FS 

, ETC)", fins* 

1290  IF  flns**"BfiCK_UPM  THEN  Copy 

1300  IF  AnsxCl.M   ":"  THEN  ftns$=" :  "8<flns* 

1310  Msusl*=Ans* 

1320  Ansf=":T14" 

1330  EDIT  "WHfiT  MASS  STORAGE  DEVICE  CONTAINS  THE  DESTINATION  MEDIUM  CT14,:T15#: 

F8,ETC>", Ans* 

1340  IF  Ans*="EACK_UP"  THEN  Copy 

1350  IF  flns*Cl,  1K>":  "  THEN  Ans*  =  "  :  "iflns-T 

1360  IF  Ansf'  Msusl*  THEN  1450 

1370  BEEF 

1330  PRINT  PRGE,LIN<  10  >,  SPA'.  10)  .  "HOW  CAN  I  COP'/  FROM  "jMsusl*;"  TO  " ;  fins* 

1390  PRINT  SPA UO),  "TRY  AGAIN" 

14O0  WAIT  2000 

1410  BEEP 

1420  WAIT  1500 

1430  BEEP 

1440  GOTO  Copy 

1450  Msus2*=Ans* 

1460  IF  Ans  =  t  THEN  Copy2 

1470  BEEP 

1430  PRINT  PAGE.  LIN'  10>  ,  TAB'.  15),  " 

1490  PRINT 

1500  PRINT  TAB'  15  >  ,  "  +  **•*-*  *********  **Wf=lRN  I  NG***********  <-♦****" 


<3) 
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1510    PR  HIT    LIN<2> 
1520    PRINT    TflB<15; 


fiUTOST     ..     DRIVER    FOR    THEPMELEX 


THIS    PPOCRAM    HILL    ERASE    ALL    FILES    ON 


;:* 


1530  PRINT  LINCI 
1540  PRINT  LHK 


,  THE-  15 >,  " 

,SFfi  10>,"USE  KevO  TO   BACK  UP   IF  THIS  IS  NOT  YOUR  DESIRE" 

1550  PRINT  LIN<l>,SPft<20>, "PRESS  CONT  IF  THIS  IS  OK  " 

1560  WAIT  20© 

1570  EEEP 

1530  INPUT  fins* 

1590  IF  flns*="BflCK_UP"  THEN  Copy 

1600  DISP  "WORKING  INITIALIZING  " ; Msus2$ 

1610  INITIALIZE  Nsus2* 

1620  IF  Msus2*C2,2]<>"T"  THEN  Copy2 

1630  DISP  "WORKING  SETTING  UP  ";Msus2* 

1640  CREATE  " DUMMY "&Msus2*,  420 

1650  CREATE  "ENDSPfl"*cMsus2*,  1 

1660  PURGE  " DUMMY "&Msus2* 

1670  Copy2: !  IN  THIS  SECTION  THE  ACTUAL  COPYING  TAKES  PLACE 

1630  RESTORE  1750 

1690  FOR  1=1  TO  3 

1700  READ  Nai»4 

1710  COPY  Nam*&Msusl*  TO  Nam*S«Msus2* 

1720  DISP  "COPYING  "jName*;"  FROM  "jMsusl*;"  TO  " ; Msus2* 

1730  WAIT  1500 

1740  NEXT  I 

1750  DATA  BDS^EY.STDr  EY, DEMO-T, DEMC-P, flUTOST, BOARDS, THEPML, TNfiMES 

1760  RETURN 

1770  ' 

1780  Set stoe: GOSUB  Error 

1790  sets  to: !  THIS  SECTION  DETERMINES  WHERE  THE  THERMELX  SYSTEM  IS  AND  SETS 

1806         !  MASS  STORAGE  TO  THAT  PLACE 

1810  ON  ERROR  GOTO  Set st oe 

1S20  S*=" LOCATION  OF  THE  THERMELEX  SYSTEM" 

1830  GOSUB  Pagehead 

1840  PRINT  "     It  is  necessary  that  the  mass  storage  media  that  contains  the  TH 

ERMELEa" 

1S50  PRINT  "system  b«  located  in  the  default  mass  storage  device.   Therefore,  th 

e  def au 1 t " 

1860  PRINT  "mass  storag 

1 1  *  r  e  d " 

1870    PRINT    "during    the 


ge    devices' 

1830    PRINT 


be 


s    will    be    converted    by    the    program.       This    should    not    be 

.<se    of    THERMELEX.       Data    files    located    on    other    mass    stor. 

sed    by    appending    the    appropriate    device    code    to    the    file- 


name. • 

ie90    Ans*=":T15" 

1900   EDIT  "CHANGE  OR  ENTER  BELOW  THE  LOCATION  OF  THERMELEX  ?",Ans* 

1910   Ans$=UPC$<Ans*) 

1920   IF  fins*Cl,n<>":"  THEN  flns$»" :  "&Ans* 

1930   IF  <Ans*C2,2]  =  "F"  ■  OP  <Ans*C2, 23="T" )  OR  CPns*C2, 23=" Y" >  OR  <f 

)  THEN  1970 

1940   BEEP 

1950   GOSUB  Err  in 

1960   GOTO  Set sto 

1970   Hsusl*=fins$ 

1930   MASS  STORAGE  IS  Msusl* 

1990   RETURN 

2000  Over  lav:   !  THIS  SECTION  PREPARES  AN  OVERLAY  FOR  THE  USER 

2010  Ans*="Y" 

2020  INPUT  "DO  YOU  DESIRE  A  PAPER  KEY-CODE  OVERLAY  <Y  OR  N)?",Ans* 
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RUTQST  ..  DRIVER  FOP  THERHELEX 

2030  IF  CUPC*<flns*>=0N")  OP  Cflnsf="BfiCK  UP")  THEM  RETURN 

2040  PRINTER  IS  O 

2050  PRINT  SPfl<l5>, "TEAR  OFF  FOP  GENERA^  PURPOSE  KEY  CODE  OVERLAY  ",LIN<2) 

2060  PRINT  "|Back_up|  EXIT   |  ENTER-]   DUMP   [~PRT  I S  |  F  P  T  IS  |  REWND  |  REWND  | 

PLACE" 
2870  PRINT  "|  <ko)   | GRAPHIC | GRAPHIC | GRAPHIC  |    16   |    0    |   :  T 1 4  |   !T15  | 

ABOVE  KEYS" 

2030  PRINT  "TEAR  HERE  THEN  PRESS  CONT " 

2090  PRINT  LlNkl; 

2100  PAUSE 

2110  PRINT  »|  <kS>   |         |         |  |         |  EDIT   |   YES   |   NO    | 

PLACE" 
2120  PRINT  "|         |         |         |  |         |  LINE   |         |         | 

BELOW  KEYS" 
2130  PRINT  LIN'S) 
2140  PRINTER  IS  16 

2150  RETURN  !  END  OF  Overlay 

2160  ! 

2170  Pagehead:  '  THIS  SECTION  PLACES  DESIRED  HEADING  ON  Pi  BLANK  CRT 
2188  PRINT  PAGE,TAB<34-LEN<S*>/2),  "**  " ; CHR*< 132) ; S*; CHR*<123) \ "    **",LIN<2) 
2190  RETURN 

2200  Errin: !  this  SECTION  ALERTS  THE  USER  TO  AN  ATTEMPT  TO  INPUT  INVALID  DATA 
2210  BEEP 
2220  PRINT  PAGE 

2230  DISP  "INPUT  OUT  OF  RANGE TRY  AGAIN" 

2240  WAIT  2500 

2250  EEEP 

2260  WAIT  1500 

2270  PEEP 

22S0  RETURN  !  END  OF  Err  in 

2290  Error: ! 

2300  PRINTER  IS  16 

2310  PRINT  LIN<1«),TAB<20), "ERROR  NUMBER  " ; ERPN: "HRS  OCCUPED  IN  LINE  " : EPRL 

2320  DISP  "PRESS  CONTINUE  WHEN  READY  TO  RESUME  PROGRAM  FLOW" 

2330  PAUSE 

2340  ! 

2350  Pagecrt : !  THIS  SECTION  BREAKS  THE  INSTRUCTIONS  INTO  PACES  FOR  THE  CRT 

2360  DISP  "Press  CONT  when  ready  t"c-r    mors  instrctions" 

2370  PRINT  LIN<B1  ink  s  > 

2380  PAUSE 

2390  RETURN 

2400  ! 

2410  Pageprt:!  THIS  SECTION  BREAK'S  THE  PRINTED  INSTRUCTIONS  INTO  PAGES 

24  20  PRINT  LIN(Blanks) 

2430  PRINT  TflBOS),  "Pg.  ";Pagenum 

2440  PRINT  LIN<1),"_ 

2450  IF  P=Lgzr,ur>,>5~ THEN  RETURN 

24  60  PRINT  LIN<4>,  TAB  ■:  37-LEN  <  S*  :■  2  >  \  CHR*<  132)  ;  S*;  CHR*<  123) ;  TAB<76 ) ;  "Pg.  ";  Pagenum 

♦  1  ,  L  I  N  <  2  ) 

2470  RETURN 

24S0  ' 

24  90  Inste:GOSUB  Error 

2500  In*t:!  THIS  SECTION  PREPARES  A  SET  OF  WRITTEN  INSTRUCTIONS  FOR  THE  USER 

2510  ON  ERROR  GOTO  Inste 

2520  S*»" INSTRUCTIONS  FOR  USE  OF  THERMELEX" 

2530  GOSUB  Pagehead 

2540  PRINT  "     A  writ,  ten  set  of  instructions  car-,  bs  prepared  for  the  user  that 

will  hel p" 
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RUTOST  ..  DRIVER  FOR  THERflELEX 

2550  PRINT  "on?  to  become  familiar  with  the  THERMELEN  system  without  actually  us 

ing  the  " 

2560  PRINT  "the  computer.   These  are  intended  to  only  supplement  the  set  of  dire 

Ct  ions" 

2570  PRINT  "given  during  the  e  ecution  of  the  system  of  programs.   These  lnstruc 

t  ions  will" 

25S0  PRINT  "normally  be  presented  on  the  screen  but  if  you  desire  a  hard  copy  of 

the" 
2590  PRINT  "instructions,  ENTER  YPPT  rather  than  Y. " 
2600  flnsT="N" 

2610  INPUT  "DO  YOU  DESIRE  PRINTED  INSTRUCTIONS  (N  or  Y  OR  YPRT)?",Ans* 
2620  IF  Ans**"BACK_UP"  THEN  RETURN 

2630  IF  (UPC#<fins*Cl, n>*"YM)  OR  <UPC*<Ans*)="N" )  THEN  2660 
2640  GOSUB  Err  in 
2650  GOTO  Inst 

2660  IF  UPC*<finsS>="N"  THEN  RETURN 
2670  PRINTER  IS  16 
2630  PRINT  PAGE 

2690  IF  UPC*<Ans*Cl,2])="YP"  THEN  PRINTER  I?  0 
2700  P=16 

2710  IF  UPC*<flns*Cl,2])="YP"  THEN  P=0 

2720  PRINT  LIN<4>,TAB<37-LEN<S*)  2>,"**  " ; CHRSC 132) ; 3*: CHRSC 128) ; "  **",LIN(2) 
2730  PRINT  "     The  THERMELEX  system  is  designed  to  perform  thermal  analysis  of 
air  cooled" 

2740  PRINT  "electronic  circuit  boards.   The  system  consists  of  three  major  progr- 
ams :  " 

2750  PRINT 

2760  PRINT  LIN<1),THB<17), "1.  RUTOST    Driver  and  system  setup." 
2770  PRINT  LIN<1),TAB<17), "2.  BORFDS    Input  and  recording  of  data." 
2788  PRINT  LIN<1),TAB<17>, "3.  THEFML    Analysis  and  output  of  results." 
2790  PRINT  LIN.  1  .. 

2800  PRINT  "The  RuTOST  program  establishes  special  key  deffinitions  and  acts  a 
s  the  driver" 

2310  PRINT  "for  the  other  programs  within  the  THERMELEX  system.   The  key  deffini 
t i ons  are  " 

2320  PRINT  "valid  throughout  the  system  except  when  entering  the  specific  coupon 
era  type" 

2330  PRINT  "codes  during  execution  of  BOARDS  .  " 
2840  B1 ant  s=0 

2850  IF  P=16  THEN  GOSUB  Paqecrt 
2860  PRINT 

2870  PRINT  "     This  driver  program  is  entered  by  placing  the  tape  in  the  primar 
y  tape" 
2SS0  PPINT  "transport  CT15)  at  the  right  side  above  the  special  keys.   You  then 

: " ,  l  I  n  <  i ) 

2890  PRINT  TAB<18>, "TYPE:  LORD  "; CHRSC34) ; "RUTOST: T15" ; CHR*<34> J", 1      Then  Pre 
ss   EXECUTE 
2900  PPINT 

2910  PRINT  "If  the  9345A  is  off.  the  driver  may  be  loaded  by  inserting  the  tape- 
in  :  T 1 5 , " 
2920  PRINT  "latching  the  AUTO  ST  Key  in  the  down  position,  and  turning  the  power 

on.  " 
2930  PRINT  "In  either  case  the  first  question  will  be  concerning  the  location  of 

the  " 
2940  PPINT  "THERMELEX  system.   The  program  will  convert  the  default  mass  storage 

device" 
2950  PRINT  "to  what  ever  device  contains  the  system.   It  is  important  to  not  cha 
nge  this" 
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2960  PRINT  "while  using  THERMELEX.   These  instructions  may  be  provided  on  paper 

or  CRT  and" 

2970  PRINT  "a  paper  keycode  over  lav  may  be  provided  as  a  guide  to  the  key  deffin 

i  t  i  ons" 

29S0  PRINT  "that  is  loaded  by  RUTQST.   m  plastic  overlay  is  available  from  Hew 

I ett -Packard" 

2990  PRINT  "'part  #  7120-6164)  which  nay  be  filled  in  using  the  paper  overlay  as 

a  gui  de. " 
3000  PRINT  "The  or  her  programs  are  entered  from  r.he  option  selection  list  (Menu) 
in  flUTOST. " 
3010  PRINT  "Normally  BOARDS  is  next  loaded  to  allow  r.he  input  of  circuit  board 

d  e  s  c  r  i  p  t  j  y  e  " 
3020  PRINT  "data  from  the  keyboard.   The  instructions  that  follow  apply  to  this 
opt i on. " 
3030  PRINT 
3040  Blanks=4 

3050  IF  P=16  THEN  GOSUB  Pagecrt 

306O  PRINT  "     The  user  should  have  either  the  circuit  board  or  a  facsimile  av  a 
i  1  ab 1 e  " 

3070  PRINT  "while  he  is  using  the  BORPBS  program.   Orient  the  board  such  that 
the  cool  i  ng" 
30SO  PRINT  "air  flow  is  from  the  left  to  the  right..   Measure  the  length  '.'defined 

in  the  air " 
309O  PRINT  "flow  direction),  width  and  thickness  of  the  circuit  board  in  mil  lime 
ter s  Cmm; . " 

310O  PRINT  "It,  ,j  necessary  to  divide  the  board  into  evenly  spaced  imaginary  reg 
ions  in  both" 

3110  PRINT  "horizontal  and  vertical  directions.   The  centers  of  each  of  these  re 
gi  ons  will" 
3120  PRINT  "form  -he  nodes  of  a  finite  difference  grid  to  be  used  in  the  thermal 

analysis" 
3130  PRINT  "of  the  electronic  circuit  board.   The  imaginary  lines  that  are  creat 
ed  divide" 

3140  PRINT  "the  circuit  board  into  regions  that  will  either  be  empty  or  will  con 
tain  one" 

3150  PRINT  "component.   fin  exception  to  this  rule  is  the  case  of  40  pin  DIP'S  wh 
ich  occupy" 

3160  PRINT  "two  regions  due  to  their  large  size.   fin  attempt  to  approximate  the 
real  circuit" 
3170  PRINT  "board  with  a  model  requires  some  flexibility  and  imagination  at  this 

poi  nt .  " 
3160  PRINT 

3190  PRINT  "     The  thermal  model  used  tor  the  analysis  will  assume  that  each  co 
mponent  i  s" 

32O0  PRINT  "centered  in  the  region.   The  following  components  can  be  included:" 
3210  PRINT 

3220  PRINT  "  DIPS  (horizontal  or  vertical',  CHIP  ChPPIEPS  " 

3230  PRINT  "  14,16,24,40  pins  16,24,40,64  pins" 

3240  PRINT 
3250  Pagenum=l 
3260  Blanks-4 

3270  IF  P=0  THEN  GOSUB  Pageprt 
3230  IF  P=16  THEN  GOSUB  Pag^crr 

3290  PRINT  "The  limitations  imposed  by  restricting  the  components  to  the  list  ab 
owe  will  " 

3300  PRINT  "prevent  some  circuit  boards  from  being  precisely  modeled  using  THEFM 
ELEX  but  " 

3310  PRINT  "are  necessary  to  provide  for  the  simple  input  of  'h*  descriptive  dat 
a.  " 
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3320  PPINT 

3330  PPINT  "     When  BOARD'S  is  loaded  using  option  #1  in  RUTQST  a  control  pa 

r  am  e  t  e  r  is" 

3340  PRINT  "assigned  which  causes  the  program  to  begin  at  the  section  which  asks 

for  the" 
3350  PRINT  ".geometric  parameters  of  the  bane  circuit  board.   You  will  be  asked  t 
o  enter  the" 

3360  PRINT  "length  of  the  board  in  millimeters.   Enter  the  length  using  either  t 
he  number" 
3370  PRINT  "pad  at  the  far  right  of  the  keyboard  or  the  numbers  above  'he  letter 

keys" 
3330  PRINT  "followed  by  the  C0NT  Key.   The  next  question  will  appear  at  the  bo 
ttom  of  the" 

3390  PRINT  "scran    while  the  length  that  was  just  entered  will  be  printed  immedi 
at  1 y  belcu" 
3408  PRINT  "the  instructions.   Answer  each  question  remembering  to  use  the  C0NT 

Key  and" 
3410  PRINT  "check  the  values  on  the  screen.   If  a  mistake  is  made  press  K0  and 

the  program" 
3420  PRINT  "will  ask  the  question  again.   Rfter  the  length,  height,  thickness,  c 
onduc t i vi  ty, " 

3430  PRINT  "number  of    regions  in  the  'X'  direction  and  the  number  of  regions  in 
the  'V  " 

3440  PRINT  "direction  have  been  entered,  the  program  will  ask  for  a  short  (less 
than  50  " 

3450  PRINT  "characters:'  descriptive  title.   This  title  will  become  part  of  the  d 
at  a  set  for" 

3460  PPINT  "the  circuit  board  and  will  be  used  to  identify  the  output.   Two  exam 
pies  are: " 
3470  PRINT 

34S0  PRINT  SPPK5J,  "TEST  BOhPD  FOP  PROJECT  UMPTYFP8TS    .^r       DESCRIPTIVE  TITLE 
#1234. ?ft  .  " 
3490  Elanks=y 

3500  IF  P=l€  THEN  G0SUB  Pagecrt 
3510  IF  P=0  THEN  PRINT 

3520  PRINT  "     The  -.cn-n    will  now  shift  fi 
de  and  draw" 

3530  PRINT  "an  empty  board  with  the  imaginai 
es.   This" 
3540  PRINT  "picture  will  remain  on  the  screen  for 


ut  it  may" 

3550  PRINT  "be  returnee 


•  he  B'Zr-a  n  u  s  i  n  g  k  2 


he  alpha  mode  to  the  graphics  mo 
g  ions  m  ar  k  e  d  o  f  f  w  •  t  h  d  o  1 1  e  d  1  i  n 
-ily  3  sec  after  it  is  drawn  b 
1  1  o  w  f  u  r  t  h  e  r  t  i  m  e  f  o  r  c  h  e  c  k 


3560  PRINT  "hi  will  return  the  screen  to  the  alpha  mode  for  the  next  section." 

3570  IF  P=fl  THEN  PPINT 

3536  PRINT  "     The  special  function  keys  will  now  be  used  to  input  the  contents 

of  each  of" 
3590  PRINT  "the  regions  defined  on  the  blank  board".   This  is  the  last  chance  to 
use  the  " 

3600  PRINT  "Back -up  option  <k0>  if  the  board  is  not  correct.   fl  new  Key  code  ove 
r 1  ay  will" 

3610  PRINT  "be  provided  if  necessary  and  must  replace  the  general  purpose  overla 
y  on  the" 

3620  PRINT  "keyboard.   The  keys  defined  in  this  section  act  as  priority  interrup 
t s  and  will" 

3630  PRINT  "allow  the  screen  to  remain  in  the  graphics  mode  as  the  components  ar 
e  def i red. " 

3640  PRINT  "fi  flashing  cross  will  move  to  each  region  in  sucession  and  wait  for 
a  valid" 
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3650  PRINT  "key  to  be  pressed.   Every  region  must  be  defined  even  if  empty.   It 

I*  ■ 

3660  PRINT  "important  to  wait  for  the  cursor  to  move  to  the  next  empty  region  be 

fore  keys" 

3670  PRINT  "are  pressed.   The  back-up  option  will  NOT  work  while  defining  comp 

onents  but" 

3630  PRINT  "there  will  be  an  opportunity  to  make  corrections  in  the  next  section 

3690  PRINT  "When  all  regions  have  been  defined,  the  screen  Mill  return  to  the  al 

pha  mode  to" 

3700  PRINT  "ask  if  there  are  any  corrections.   fit  this  point  k2  will  return  the 

graph i cs" 

3710    PRINT    "to    the    screen    and    kl    will     return    the    alpha    mode    to    the    screen." 

3720    Bl anks=0 

3730  IF  P=16  THEN  gosuB  Pagecrt 

3740  PRINT 

3750  PRINT  "To  make  corrections,  the  numbers  of  ALL  incorrect  regions  are  ente 

red  and  the" 

3760  PRINT  "components  erased.   When  there  are  no  more  to  be  removed,  pre n  ing 

COMT  with" 

3770  PRINT  "no  entry  will  allow  all  empty  regions  to  be  redefined.   R  flashing  c 

ross  will" 

3780  PRINT  "visit  each  empty  region  followed  by  further  correction  opportunity  a 

s  needed. " 

3790  PRINT  "When  the  board  is  correct,  the  general  purpose  key  deffi nit  ions  are 

returned" 

3800  PRINT  "and  the  normal  key  code  overlay  should  be  replaced. 

3310  PRINT 

3320    Pagenum=2 

3830    Blanks=ll 

3S40    IF    P=0    ThEN    G0SUB    Pageprt 

3S50    PRINT    "  DIP    components    generally    are    produced    in    two    different    case    styl 

es,     ceramic" 

3860  PRINT  ".;<■     plastic.   These  are  different  sizes  and  this  affects  the  thermal 

analysis" 

3370  PRINT  "performed  in  THEPHL.   The  next  section  provides  several  methods  to 

spec  i  fy  the" 
3830  PRINT  "case  styles  for  the  DIP'S.   These  are:  1)  all  plastic... 2)  all  ceram 
i  c  .  .  .  " 

3890  PRINT  "3)  majority  ceramic  ''user  specify  which  are  p  1  ast  i  c  )  .  .  .  4  j  majority  c 
eramic  (user" 
39O0  PRINT  "specify  which  are  c e r am i c ) . . . 5)  user  specify  case  style  for  each  DIP 

component . " 
3910  PRINT  "In  general,  option  3  or  4  require  the  user  to  input  the  number  of  re- 
gions that" 

3920  PRINT  "are  either  plastic  or  ceramic  and  are  useful  ir.r    specifying  the  Styl 
e  when" 

3930  PRINT  "when  only  a  few  are  different.   When  there  is  a  nearly  equal  mix,  op 
t  i  on  5" 
3940  PRINT  "requires  a  1  for  plastic  style  and  a  0  for  ceramic  style.   After  all 

DIP'S  are" 
3950  PRINT  "specified,   a  list  of  region  numbers  is  displayed  on  the  screen  (in 
al pha  node) " 

3960  PRINT  "with  the  plastic  cases  in  inverse  video  (numbers  black  on  light  back 
round)  and" 

3970  PRINT  "ceramic  cases  in  normal  video.   Correct  as  needed  and  again  no  corre 
ct  i  ons  i  s" 
39S0  PRINT  "indicated  by  pressing  CONT  Key  with  no  entry." 
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3990  IF  P=e  THEN  4040 

4090  DISP  "PRESS  CONT  WHEN  READY  FOR  ANOTHER  PAGE  OF  INSTRUCTIONS" 

4010  PRINT  LINC2) 

4020  PAUSE 

4030  GOTO  4050 

4O40  PRINT 

4050  PRINT  "     The  next  section  concerns  the  thermally  conductive  paths  that  ar 

e  added  " 

40&0  PRINT  "to  the  circuit  board.   Both  the  electrical  conductors  and  any  therma 

1  conduct  i  on" 

4070  PRINT  "paths  or  rails  will  trznsfer    heat  between  regions.   When  any  conduct 

i  on  " 

40£0  PRINT  "rails  are  present,  the  effects  oi     the  electrical  leads  will  be  negl i 

gable  and" 

4090  PRINT  "ignored.   To  aid  in  the  modeling  of  the  conduction  rails  it  is  neces 

&  ar  y  to" 

4100  PRINT  "introduce  the  concept  of  a  Conduction  Path  Unit.  (CPU).   Each  of  the 

rails  is" 

4110  PRINT  "modeled  as  containing  an  integer  number  CPU's.   The  width  of  a  CPU  i 

s  defined  " 

4  120  PRINT  "to  be  .1  ram  while  the  user  must  define  the  thickness  (mm)  and  therma 

1  " 

4130    PRINT    "conductivity    'Uittifl-O    of    the    material     that     forms    the    conduction    r 

ai Is." 

4140  PRINT  "The  thermal  conductivities  of  several  aluminum  alloys  and  copper  are 

presented" 
4150  PRINT  "to  aid  the  user.   If  there  are  no  conduction  rails  pr-mnt  .  the  elec 
tncal  lead" 

4  1  i 0  PRINT  "description  is  entered  in  the  m u  c  h  t  h  e  s  am  e  w ay :  t  h  i  c  k n e  s  s  < m m  > ,  av  e 
rage  width" 

4170  PRINT  "..mm)  and  the  conductivity  ( Wat  t  s-  H-C  > .   This  will  default  to  334  Wat 
tS''M-C,  the" 

4  180  PRINT  "value  for  copper.  This  is  the  most  commom  material  for  traces  on  cir- 
cuit boards . " 
4190  IF  P=0  THEN  4240 
4200  PRINT  LIN<2) 

4210  DISP  "PRESS  CONT  WHEN  READY  FOR  ANOTHER  PACE  OF  INSTRUCTIONS" 
4220  PAUSE 
4230  GOTO  4250 
4240  PRINT 

4250  PRINT  "    In  the  next  section  the  numbers  of  leads  or  CPU's  that  cross  the 
l nt  er i  or" 

42?*  PRINT  "boundaries  between  regions  are  entered.   The  display  will  alternate 
between" 

4270  PRINT  "the  graphics  and  the  alpha  mode  with  the  flashing  cross  or  cursor  mo 
wing  to  the" 

4230  PRINT  "appropriate  location  on  the  board  to  be  defined.   The  graphics  will 
remain  on" 

4290  PRINT  "the  screen  for  two  seconds  and  then  return  to  the  alpha  mode  for  an 
i  r.put  1  i  ne.  " 

4300  PRINT  "It  ijill  ask  for  the  input  of  the  number  of  leads  or  CPU's  between  tw 
o  spec  i  f  i  c  " 

4310  PRINT  "regions;  ie  if  the  conduction  rail  is  6.7    rum  wide,  there  are  £?    CPU' 
s  in  that" 

4320  PRINT  "path.   Special  function  key  3  fk3>  may  be  used  to  dump  the  graphics 
to  the" 

4330  PP1NT  "printer  for  a  hard  copy  of  the  graphics  picture.  It  is  important  to 
account " 
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4340  PRINT  "for  the  leads  on  both  sides  of  r  he  circuit  board  when  entering  this 
number .  " 

4350  PRINT  "fis  before  follow  all  entries  with  the  CONT  key.   Eec ause  of  the  pr- 
ogram flow" 

4360  PRINT  "in  this  section,  the  Back_up  option  <k8)  will  not  function.   Do  not 
attempt  to" 

4370  PRINT  "make  corrections  or  return  to  a  previous  section  until  after  all  ent 
ries  have" 

43S0  PRINT  "been  made.   Hll  corrections  will  be  made  in  the  next  section." 
4390  PRINT 

4400  PRINT  "     In  the  thermal  model  used  by  the  THEPMELE:;  system,  all  heat  is  t 
ransf ered" 

4410  PRINT  "to  the  cooling  air  stream.   None  of  the  heat  is  transferred  out  the 
edges  of" 
4420  PRINT  "the  circuit  board.   For  this  reason,  the  sides  of    the  regions  at  the 

edges" 
4430  PRINT  "will  be  labeled  with  0  to  indicate  that  no  leads  or  CPU's  cross  this 

boundary. " 
4440  PRINT  "When  this  heat  loss  from  the  board  edges  and  connector  is  neglected, 

the  resul is" 
4450  PRINT  "will  be  conservative.   These  assumptions  may  be  unacceptable  for  cir- 
cuit boards " 
4460  PRINT  "used  in  installations  where  a  significant  fraction  of  the  total  heat 

di ssipated" 
4470  PRINT  "is  conducted  aw  a"  from  the  board  through  card  guides  or  connectors." 
4430  PRINT 
4  4  90  P  ag  e  n  u m  = : 
4500  Blanks =6 

4510  IF  P=0  ThEN  G0SUB  Pageprt 

4520  PRINT  "     When  all  regions  have  been  specified,  the  screen  will  return  to 
the  alpha" 

4530  PRINT  "mode  and  ask  if  an'.'  corrections  are  needed.   Use  the  special  functio 
ns  keys" 

4540  PRINT  "as  required  to  either  return  to  gr aph i cs  '  al oh a  or  produce  a  hard  cop 
y .   Enter" 
4550  PRINT  "the  numbers  of  the  connect i ng" regi ons  that  have  incorrect  values  for 

the  numbers" 
4560  PRINT  "of  leads  or  CPU's.   his  before,  CONT  with  no  entry  signals  the  lack 

of  further" 
4570  PRINT  "corrections.   The  physical  description  of  the  circuit  board  is  now  c 
ompl et  e .  " 

4580  PRINT  "All  that  remains  is  to  specify  either  the  component  power  levels  <Wa 
tts)  or  the" 

4590  PRINT  "junction  temperatures  for  each  component." 
4600  PRINT 

4610  PRINT  "     When  the  maximum  junction  temperatures  for  the  components  are  sp 
ec  i  f i  ed, " 

4620  PRINT  "the  maximum  steady  state  power  le"el  that  will  result  in  that  temper- 
ature is" 

4630  PRINT  "calculated  for  each  component.   When  the  stead','  state  power  of  each 
component " 
4640  PRINT  "is  specified,  the  stead','  state  junction  temperatures  are  calculated. 

4650  PRINT  "In  either  situation  it  is  also  necessary  to  specify  the  case  to  June 

t  ion" 

4  660  PRINT  "thermal  resistance  R.j_C  'Peg  C  -  W  at  t  :  .   The  thermal  model  uses  the  co 

mponent " 

4670  PRINT  "surface  temperatures  for  heat  transfer  calculations  and  Rj_C  provide 

s  the  1 i nk " 
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4680    PRINT    "between    the    component     power,     junction    temperature    and    the    surface    t,e 

mper  ature . " 

4690  PRINT  "If  Ri_.:  is  specified  as  zero  >>0>  then  the  case  surface  temperature  w 

ill  be" 

4700    PRINT    "equal    to    the    junction    temperature." 

4710    PRINT 

4720  PRINT  "     After  all  component  power  levels  or  junction  temperatures  have  b 

een  set , " 

4730  PRINT  "the  complete  list  is  displaced  in  the  alpha  mode  to  allow  checking  a 

r.d" 

4746  PRINT  "changes.   Corrections  are  inserted  by  entering  the  incorrect  region 

number  then" 

4750  PRINT  "entering  the  correct  values  in  responce  to  the  question  asked.   If  d 

es  i  red" 

4760  PRINT  "a  hard  copy  of  all  data  for  the  circuit  board  will  be  provided  from 

the  pr i nter . " 

4770  PRINT 

4780  PRINT  "     The  next  section  allows  recording  r,f    the  circuit  board  descripti 

ve  data  " 

4790  PRINT  "on  to  a  mass  storage  media.   This  should  be  done.   The  recorded  data 

file  may  be" 
4800  PRINT  "accessed  at  some  later  time  either  in  B 0 ft P D S  for  editing  and  chang 
es  or" 
4310  PRINT  "directly  in  THEPHL  to  perform  the  thermal  analysis.   B0ARB3  uill 

create  the  " 
4320  PRINT  "the  necessary  data  file  to  record  the  description  but  if  a  data  file 

of  " 
4830  PRINT  "sufficient  size  exists,  it  may  be  used.   Fill  old  data  in  that  file  w 
i  1  1  " 

4340  PRINT  "be  permanently  lost.   Information  concerning  the  minimum  data  file  s 
i  z  e  is" 

4850  PRINT  "presented  by  the  program  to  help  in  this  decision.  " 
4360  PRINT 
4370  PRINT  "     Built  in  error  traps  throughout  the  entire  THERMELEX  system  will 

s  av e  the  " 
4330  PRINT  ">Mir    from  loss  of  previous  inputs  in  case  of  errors.   Information  ab 
out  the" 

43?u  PRINT  "error  is  presented  and  using  the  error  message  guide  on  the  plastic 
slide  ou»s" 

43O0  PRINT  "below  the  =,cran,     the  user  may  be  able  to  m  ak  e  corrections.   When  re 
ady,  as" 
4910  PRINT  "signaled  by  CONT  liy,     program  flow  returns  to  a  previous  point  and 

oper  at l on" 
4920  PRINT  "resumes.  The  data  recording  section  where  file  names  are  entered  is 
p  ar  t  i  c  u  1  ar  1  y  " 

4930  PRINT  "prone  to  spelling  irrori    or  entry  of  incorrect  mass  storage  unit  spe 
c  i  f  i  ers. " 

4940  PRINT  "Be  careful  but  remember  th  system  will  attempt  to  catch  those  inevit 
able  errors" 

4950  PRINT  "and  act  on  them  before  they  become  fatal." 
4960  PRINT 
4970  PRINT  "    The  end  of  a  cycle  through  BOARD;  is  now  complete.   As  the  page 

he  adi  rig" 
4930  PRINT  "that  will  appear  states,  WHAT  NOW  ? .   More  work  with  BOARDS  is  pos 
sible.   The" 

4990  PRINT  "user  may  either  input  another  circuit  board  description  from  the  key 
board  and" 
50O0  PRINT  "record  that  to  mass  storage;  or,  he  may  retrieve  a  different  circuit 

board" 

(12) 

89 


RUTQST  ..  driver  fop  .  thefhele:: 

5010  PRINT  "descriptive  data  set  from  mass  storage  for  checking  and  corrections. 

5020  PRINT  "Another  option  is  to  perforin  the  thermal  analysis  of  the  circuit,  boa 

rd  using" 

5030  PRINT  "the  third  program  in  the  THERMELEN  system,  THERML. " 

5040  PRINT 

5050    Pagenum=4 

50&0    El anks=6 

5070  IF  P  =  0  THEN  GOSUB  Pageprt 

50S0  PRINT  "     Set  ore  an  explanation  of  the  program  flow  in  THERML.  some  expl 

anat  i  on  of" 

5090  PRINT  "the  use  of  B0RRDS  to  edit  previously  recorded  data  files  is  in  ord 

er.   This" 

5100  PRINT  "option  may  be  accesed  through  the  use  of  flu TOST  as  explained  on  th 

e  first  page" 

5110  PRINT  "of    these  instructions  or  the  BOARDS  program  may  be  directly  enters 

d  into  the" 

5120  PRINT  "computer  with  the  command:      LOAD  "; CHR*< 34  •;  "BOARDS: T 15 ,  1 "; CHR*  •  .3 

4) ; "         then  Press  EXECUTE  " 

5130  PRINT 

5140  PRINT  "Of  course  the  appropriate  mass  storage  specifier  should  be  used  both 

here  and" 
5150  PRINT  "Linen  entering  the  file  name  for  any  previously  recored  data  file.   S 
ee  the  " 

51c*0  PRINT  "Operating  and  Programming  Manual  for  the  9845  if  the  above  is  not  cl 
ear .  " 

5170  PRINT 

51S0  PRINT  "    Following  the  input  of  the  file  name  that  contains  the  descriptiv 
e  dat  a,  " 

5190  PRINT  "the  data  file  will  be  read  and  all  the  variables  that  are  normally  s 
et  from" 

5200  PRINT  "the  keyboard  will  be  defined.   The  circuit  board  will  be  drawn  on  th 
e  screen" 

5210  PRINT  "in  graphics  mod*  and  changes  to  the  specific  components  that  occupy 
the  regions" 

5220  PRINT  "may  be  made.   However;  it.  js  not  possible  to  change  the  number  of  re 
g i o  n s  in" 

5230  PRINT  "any  way.   This  major  revision  of  the  board  must  be  done  as  a  new  boa 
rd  i  nput . " 

5240  PRINT  "Correction  to  case  styles  is  next  and  then  the  numbers  of  leads  or  C 
PU's  " 

5250  PRINT  "between  regions  will  be  labeled  on  to  the  graphics  picture  of  the  ci 
rcuit  board" 
52&0  PRINT  "Changes  will  be  allowed  to  the  numbers  but  not  the  physical  sires  of 

the  leads" 
.5270  PRINT  "or  CPU's;   this  also  requires  a  new  board  input.   A  list  of  power  le 
vels," 

5280  PRINT  "junction  temperature  and  case  to  junction  thermal  resistances  is  pre 
sented" 

5290  PRINT  "for  checking  and  corrections  as  required.   This  edited  descriptive  d 
at  a  set  may" 
5300  PRINT  "then  be  recorded  either  back  to  the  same  data  file  or  to  a  new  file. 

5310  PRINT 

5320  PRINT  "    There  are  two  demonstration  data  files  included  in  the  THERHELEK 

system. " 

5330  PRINT  " D E M 0 - T  contains  a  typical  circuit  board  description  needed  t o  dete 

rm l ne  the" 
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5340  PRINT  "steady  state  junction  temperature  <  i  e  component  power  levels  are  spe 

cif  ied>  .  " 

5350  PRINT  "DEMO-P  contains  a  different  circuit  board  for  which  maximum  powers 

levels  are  " 
5360  PRINT  "determined  in  the  thermal  analysis.   The  us*  of  either  of  these  can 
be  helpful  " 

5370  PRINT  "to  familiarize  the  user  with  both  3 Oft EDS  and  'he  last  program  in  t 
h€  THERELEX" 

53S0  PRINT  "system,  THERML. " 
5390  PRINT 
5400  PRINT  "    THER,'1L  creates  the  mathematical  model,  performs  the  analysis  of 

the  c  i  rcui  t  " 
5410  PRINT  "board  and  produces  the  output.   The  variables  used  in  the  thermal  mo 
del  are  set" 

5420  PRINT  "using  either  a  data  file  from  mass  storage  or    they  are  passed  in  a  c 
or.imon  bl  ock  " 

5430  PRINT  "when  THEPML  is  loaded  from  BOARDS.   There  are  no  correction  oppo 
rtuni ties" 

5440  PRINT  "for  the  circuit  board  descriptive  variables  in  this  program." 
5450  PRINT 

5460  PRINT  "     Assuming  that  THERML  has  been  loaded  from  the  What  Now  option 
list  in  " 

5470  PRINT  "EOF!- PS,  the  first  input  is  the  inlet  t  emperat  ure<deg  C '■  of  the  coo 
1 i  ng  air." 

54S0  PRINT  "After  this  is  entered  and  printed  on  the  ir.reen,     the  volumetric  air- 
flow" 

5490  PRINT  "rate  <MA3/'sec)  is  requested.   Since  the  same  set  of  general  purpose- 
spec  i  a  1  " 

5500  PRINT  "function  keys  are  valid  in  this  program,  the  keycode  overlay  should 
remain  and" 
5510  PRINT  "k0  (Back-up)  may  be  jsed  to  return  to  a  previous  question  when  input 

errors    ar e " 
5520  PRINT  "made.   The  next  input  is  the  spacing  between  the  circuit  boards  for 
rack  mount  =  .  " 

5536  PR  IN"   This  spacing  and  the  board  height  determines  the  sire  of  the  i mm  agin 
ary  air  duct " 

5540  PRINT  "containing  the  circuit  board.   When  combined  -nth  the  amount  of  air- 
flow for" 

5550  PRINT  "each  circuit  board  the  air  velocity  and  heat  transfer  correlations  a 
re" 

5560  PRINT  "determined.   Since  all  the  power  dissipated  on  the  circuit  board  mus 
t  be" 

5570  PRINT  "removed  by  the  air  flow,  it  is  important  for  these  values  to  be  as  p 
r  e  c  i  s  e " 

55S0  PRINT  "as  possible.   After  the  thermal  analysis  is  complete,  there  will  be 
opport un i  tes" 

5590  PRINT  "to  alter  these  important  air  flow  paramters  and  examine  the  effects 
on  the" 

5600  PRINT  "output .  " 
5610  PRINT 
5620  Pagenum*5 
5630  Blanks=6 

5640  IF  P=0  THEN  GQSUB  Pageprt 

5650  PRINT  "     Before    the  analysis  begins,  the  convergence  criteria  must  be  set 
.   The" 

5660  PRINT  "analysis  consists  of  assuming  a  solution  then  defining  a  set  of  simu 
taneous" 
5670  PRINT  "equations  '-.one  for  each  region)  which  are  solved  for  either  the  powe 
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5689  PRINT  "or  the  junction  temperatures.   The  results  are  compared  to  the  previ 
ous  resul t s" 

5690  PRINT  "'the  assumed  values  for  the  first  itteration)  and  if  the  largest  dif 
f erence " 

5700  PRINT  "between  these  values  is  less  than  a  maximum  set  by  the  user,  the  res 

ults  are" 

5710  PRINT  "presented  in  graphics.   If  the  convergence  criteria  is  not  met,  the 

resul is" 

5720    PRINT     "are    printed    to    the    screen     in     alpha    mode     and     another     itteration     is     en 

tered.     u 

5730  PRINT  "The  closer  to  zero  one  sets  the  convergence  criteria,  the  longer  the 

proces  s" 
5740  PRINT  "takes.   Default  values  for  the  convergence  criteria  are  .5  Deg  C  or 
\\    change  in" 

5750  PRINT  "power.   These  typically  require  2  to  6  itterations  before  they  are  m 
et.   Each" 

5760  PRINT  "itteration  takes  5  to  103  sec  to  perform  depending  on  the  number  of 
regi  ons. " 

5770  PRINT  "fl=  the  solution  precedes  the  user  may  view  the  intermediate  results 
to  watch" 

5730  PRINT  "the  progress." 
5790  PRINT 

5S00  PRINT  "     The  final  results  are  written  onto  a  picture  of  the  circuit  boar 
d  with" 

5310  PRINT  "each  region  containing  component  type,  junction  temperature,  power  1 
evel  and" 
5320  PRINT  "case  temperature.   Empty  regions  contain  only  the  temperature  of  the 

circuit" 
5830  PRINT  "board.   Those  junction  temperatures  within  5";  of  the  maximum  and  tho 
se  power" 

5340  PRINT  "le"els  within  5":  -.'f    the  minimum  are  starred  H**)  for  easy  rifir^nc^ 
as  trouble" 

5250  PRINT  "spots..   This  graphics  output  is  automatically  dumped  to  the  printer 
to  insure" 

5860  PRINT  "that  a  hard  copy  of  the  results  exist." 
5370  PRINT 

5880  PRINT  "     This  completes  a  cycle  through  THEPML  and  again  the  question  i 
s  What  Now  '"' 

5890  PRINT  "  E  C  f>.  F  D  ■£,  may  be  returned  to  core,  to  make  changes  to  the  circuit,  boa 
r  d  or" 

5900  PRINT  "to  enter  a  new  circuit  board.   Another  analysis  of  the  same  circuit 
board" 

5910  PRINT  "may  be  done  with  a  new  set  of  air  flow  paramters  or,  a  new  circuit  b 
oard" 
5920  PRINT  "descriptive  data  set  may  be  read  in  from  mass  storage  with  THEPNL . 

5930  PRINT 

5940  PRINT  "     fin  additional  option  is  sensitivity  analysis.   Plots  of  the  maxi 
m  u  m  " 

5950  PRINT  "Tjunc  vs.  fiir  Flow  Pate  or  the  minimum  Power  vs.  fiir  Flow  Rate  may  b 
e  produced. " 

5960  PRINT  "A  maximum  air  flow  rate  is  specified  and  five  separate  analyses  are 
per  formed" 

5970  PRINT  "and  the  results  plotted  and  dumped  to  the  printer.   The  maximum  air- 
flow rate" 

5980  PRINT  "specified  should  be  an  integer  multiple  of  five  to  make  for  better  1 
ook  i  ng  a/.e-  .  " 

5990  PRINT  "It  is  also  possible  to  produce  hard  copys  of  the  printed  results  for 
each  of  the" 
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6000  PRINT 
am  flow" 
6010  PRINT 
6010  PRINT 
6030  PRINT 

opt i on. " 
6040  PRINT 
that  the" 
6050  PRINT 


air  flou  rates.   Upon  completion  of  this-  sensitivity  analysis,  progr 
returns  to  the  option  list  previously  discussed." 

Termination  should  be  done  under  program  control  using  the  final 
This  will  insure  that  the  standard  key  definitions  are  returned  and 
proper  graphics  parameters  are  set  for  the  next  user  of  the  computer 
TAB<32) , "GOOD  LUCK" 


6060  PRINT  LIN 

607O  Pagenuw=6 

6030  El ank  s*16 

6090  IF  P=0  THEN  G0SUB  Paqeprt 

6100  PRINTER  IS  16 

61 10  RETURN 
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APPEND I : 


BOARDS         DATA  INPUT  AND  STORAGE  PROGRAM  FOR  .... 

r*************~«-f  *         THERMELEX  ***■*****-*♦*«■<■-******-»-►*-►-*-►?■■ 

A  SYSTEM  OF  PROGRAMS  FOR  THE  HP  9S45 
TO  PERFORM  THERMAL  ANALYSIS  OF  ELECTRONIC  CIRCUIT  BOARDS 


PREPARED  AT  THE  NAVAL  POSTGRADUATE  SCHOOL  MONTEPE', 
R.  A.  F0LT2  LCDR  USN 
Mill   BOARDS   jIMJ  ^         ..... 


55 
60 
65 

?e 

75 

80  BOARDS  J 

85 

90    ' 

95    OPTION  BASE  1 

100   PRINTER  IS  16 

105   CON  Map, Name*,Pict lbl*C5©3, Case* C 59] , Bdl , Bdh, Sea,  Thick  b, Tern  sol .Kb, Al , Kl 

110   COM  SHORT  X<50>,Y<50>,Tj<50),Pow<50>,Rj_c<59) 

115   COM  INTEGER  I type<50) , Nl <4, 50) , Nxr , Nyr , Nreg 

120   SHORT  L«<50  -  ,Hc  50) 

125   INTEGER  Npin<30 < 

130   DIM  S*C703,Msus*t22] 

135   IF  Map  =  0  THEN  St irt 

140  Start_over:  ON  Hap  GOTO  01 dpict_correct , Bdpi ct , 01  dpi  ct 

145   DISP  "MAP=  M;Map 

150   PAUSE 

155  Starts:  GOSUB  Error 

160  Start:  3*="THERMAL  ANALYSIS  PROGRAM  FOR  ELECTRONIC  CIRCUIT  BOARDS" 

165   ON  ERROR  GOTO  St  arte 

170   GOSUB  Pagehead 

175   PRINT  "   This  program  allocs  the  user  to  input  a  detailed  description  of  an 

e 1 ec  t ron  i  c " 
ISO   PRINT  "circuit  board  and  perform  a  thermal  analysis  to  predict  either  the  s 
teady  state" 

185   PRINT  "temperatures  of  the  components  if  pouzr    levels  are  given  or  the  comp 
onent  power" 

190   PRINT  "level  that  will  result  in  the  steady  state  temperature  specif  led  bv 
the  user." 
195   PRINT  "The  description  can  come  from  a  data  file  previously  created  by  this 

program" 
20O   PRINT  "or  can  be  entered  at  this  time." 
205   AnsJ="H" 

210   INPUT  "DO  YOU  DESIRE  TO  RETRIEVE  A  PREVIOUSLY  STORED  DESCRIPTION  (N  OR  Y>", 
Ans* 

215   IF  An**»"BACK_UP"  THEN  160 

220   IF  (UPC*<Ans*Cl, i3)="N">  OP  <UPC*<Rns*C 1, 13 >="Y" >  THEN  235 
225   GOSUB  Err  in 
230   GOTO  St  art 

235   IF  UPC*<Ans*El, 13)="N"  THEN  Bdpict 
240   GOTO  Oldpict 
245  01dpi cte:  GOSUB  Error 

250  Oldpict:  '  THIS  SECTION  RETRIEVES  A  BOARD  DESCRIPTION  OFF  A  MASS  STORAGE  DEV 
ICE 
255   OH  ERROR  GOTO  01  dpi cte 
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t  1  on  must  hi(>:  D?en 
be  ;ure  to  include 


260   0lde=PI 

265   S*»"B0RRD  DESCRIPTION  FROM  MASS  STORAGE  DEVICE" 
270   GOSUB  Page  he  id 

275   PRINT  "     The  data  file  containing  the  board  desc 
ored  by  " 

230   PRINT  "this  program.   Enter  the  data  file  name  bel 
e  mass  " 

285   PRINT  "storage  unit,  specifier  if  needed,  eg   : T14, : F8, : Y12   etc" 
290   EDIT  'UNDER  WHAT  FILE  NAME  IS  THE  CIRCUIT  BOARD  FILED  (input  or  change  b* 
u) " , Name* 

295   IF  Name* ="BACK_UP"  THEN  166 
300   IF  Name** "NULL*"  THEN  01dpi  ct 
305   ASSIGN  #5  TO  Nar„e* 
310   READ  1*5,  1 

315  READ  #5:  Name*  ,  Pi  ct  1  bl  *,  Case*,  Bdl  ,  Bdh,  Sc  a,N  ,r,  Nyr,  Nreo,  Th  i  ck_b,  Kb,  Al  ,  K  1  ,  Tei 
ol 

320   FOR  1  =  1  TO  Mr  eg 

325   READ  *5; Itype<I > , X< I > , Y< I > , Tj <  I),Pow;i  ),Rj_c<I> 
330   FOR  J=l  TO  4 
335   READ  »5;N1 <J, I> 
340   NEXT  J 
345   NEXT  I 

350  0 1  dp i c  I  _c  orrec  %  \     GOSUB  Bo  ar d_p i  c  t 
355   GRAPHICS 
366   FOR  1=1  TO  Nreg 
365   IF  1=1  THEN  385 
370   IF  KNxr+1  THEN  330 

375   IF  <Itype<I)=10)  AND  C  ItypeU-Nxr )*10:  THEN  400 
380   IF  (Itype<I)*9)  AND  < Itype< I-i >*9)  THEN  400 
335   IF  ItypeC  I  >=0  THEN  GOSUB  Empty 
390   IF  Itype<I>=0  THEN  400 

395   ON  It yp« i  I )  GOSUB  Hd i p 1 4 . Vd i  p 1 4, Hd  i  p 1 6, Vd i p 1 6 , Hd  i p24, Vd i p24, Erro 1 de , Err  o 
,  Hd  i  p40 ,  Vd  i  p40,  Fl  at  16,  Fl  a'  24  ,  F  1  a'  JO  .  F 1  ar  64 
400   NEXT  I 
405   GOSUB  Bd_cha 

410   IF  Ans*="BACK_UP"  THEN  Oldpict 
415   GOSUB  Case_typ 
420   IF  Rns**:'ERCk_UP"  THEN  405 
425   GOSUB  Le  ad  =  _old 
430   IF  Ans*="EACk_UP"  THEN  415 
435   GOSUB  Leads  cha 
440   IF  fins** "BACK  UP"  THEN  425 
445   GOSUB  Tempi n 
450   IF  Rns**"BRCK_UP"  THEN  435 
455   GOSUB  Stow 

460   IF  Rns*="BRCK_UPH  THEN  445 
465   GOSUB  What_now 
470   IF  Rns*="BRCK_UP"  THEN  455 
475   GOTO  St  art _0"€r 
480  Bdpicte:  GOSUB  Error 

485  BdpiCt:!  THIS  SECTION  IS  THE  NORMAL  ENTRY  POINT  FOP  NEW  BOARD  INPUT 
490   ON  ERROR  GOTO  Bdpicte 
495   01de=0 

500   S*»BKEYBORRD  INPUT  OF  CIRCUIT  BOARD  DESCRIPTION" 
505   GOSUB  Page  he  ad 

510  PRINT  "     This  routine  draws  a  blank  circuit  board  and  divides  is  into  r-. 
ons  of  " 

515  PRINT  "interest  as  determined  by  the  user.  There  are  a  maximum  of  50  regi. 
avai ab l e " 
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520  PRINT  "with  one  component-region.   The  circuit  board  is  then  presented  on  th 
e  " 

525  PRINT  "CRT  either  normal  =  i  ze  or  scaled  up  or  down  by  some  integer  if  larger 
or  smaller" 

530  PRINT  "than  I 28mm  by  170mm.  The  specific  components  ire  selected  with  the  us 

er  keys" 

535  PRINT  "in  the  upper  right  corner  of  the  keyboard.   If  you  do  not  have  i  key 

code  " 

540  PRINT  "overlay  for  the  keys  the  program  will  provide  cm    for  you.   Informati 

on  to  be" 

545  PRINT  "used  in  the  thermal  analysis  routine  is  entered  in  this  section  of  th 

e  program. " 

550  PRINT  "You  will  have  the  option  of  recording  the  board  description  onto  mass 

storaoe  .  " 
555  PRINT 

560  PRINT  "Conductivity  of  the  board  will  default  t0  .2942  Wat  t  *,■  n-degK  with  no 
entry. " 
565  PRINT  LINt 1', "NOTE  !  LENGTH  IS  DEFINED  TO  BE  IN  THE  DIRECTION  OF  AIR  FLOW 

570  NAT  Itype=2ER 

575  fins**" NULL" 

5S0  INPUT  "WHAT  IS  THE  LENGTH  OF  THE  BOARD  IN  M ILLIMETERS" , Ans* 

535  IF  Ans*="BACK  UP"  THEN  Start 

590  IF  AnsI="NULL"  THEN  575 

595  Eoi  =ABS<VAL< Ans*>  > 

600  PRINT    LIN<1),TAB< 10), "LENGTH    =";Bdl;"    mm" 

605  Scalex»INT<Bdl     171  '  +  1 

610  IF    EdK171'2    THEN    Scalex».5 

615  Ans*a"NULL" 

620  INPUT  "WHAT  IS  THE  HEIGHT  OF  THE  BOARD  IN  MILLIMETERS" ,  Ar.s* 

625  IF  Rns**"BACK  UP"  THEN  Bdpict 

630  IF  Ans*="NULL"  THEN  615 

635  Edh  =  AES' VAL<Ans*  ■  • 

640  PRINT  LIN<1),TAB< 10), "HEIGHT  =":Edh;"  mm" 

645  Sc al ey*INT<Bdh/121 >♦ 1 

650  IF  Bdh< 121/2  THEN  Scaley=.5 

655  Sca=MAX<Scaley,Scalex)    !  SCALE  FACTOR  TO  FILL  CRT  DISPLAY 

660  IF  Sca=l  THEN  686 

665  IF  Sca<l  THEN  PRINT  LIN<2) , "****NQTE****      DUE  TO  SMALL  SIZE  OF  THE  BOA 

RD  CRT  DISPLAY  WILL  BE  2X  SIZE" 

670  IF  Seal  THEM  PRINT  LIN<2> , "****N0TE****      DUE  TO  LA^GE  SIZE  OF  THE  BOA 

RD  CRT  DISPLAY  WILL  BE  1  ■  " ; Sc a; " S I ZE " 

675  Ans*="NULL" 

680  INPUT  "WHAT  IS  THE  NUMBER  OF  DIVISIONS  '.REGIONS)  IN  THE  HORIZONTAL  DIRECT  10 

N", Ans* 

6S5  IF  Ans**"BACK  UP"  THEN  615 

690  IF  Ans*="NULL"  THEN  675 

695  Nxr  =  ABS'"/AL<Ans  t  >) 

70O  PRINT  LIN<1),TAB<10), "*  Xreg  =";Nxr 

765  AnsJ="NULL" 

710  INPUT  "WHAT  IS  THE  NUMBER  OF  DIVISIONS  ■'.REGIONS;  IN  THE  VERTICAL  DIRECTION" 

,  Ans* 

715  IF  Ans*="BACK  UP"  THEN  675 

720  IF  Ar,sJ  =  "NULL"  THEN  705 

725  Nyr-ABS<VAL<Ans*)) 

730  PRINT  LIN<1),TAB<10),"«  Yreg  =";Nyr 

735  Nreg  =  Nxr*N-,'r  !  Nreg  =  NUMBER  OF  REGIONS 

740  IF  Nreg  51  THEN  770 
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=  " 

;Thi 

ck 

iBUi 

:  t  i  \ 

'IT 

=  " 

;Kb 

ACTERS) 

BOARDS  ..  DATA  INPUT  FOP  THERHELE: 

745   BEEP 

750   PRINT  PAGE, LIN<15>,SPfi<30>, "TOO  MANY  REGIONS" 

755   UAIT  1000 

760   GOSUE  Errin 

765   GOTO  Bdpicx 

770   flns**"NULL" 

775   INPUT  "WHAT  IS  THE  THICKNESS  OF  THE  BOARD ?< mm) " , Rr 

7S0   IF  Ans**"BACK_UP"  THEN  705 

735       IF    flns*=HNULLM    THEN    770 

790      Thick    b=ABS<VAL<Ans*) ) 

795       PRINT    LIN'  1  ;■  ,  TAB<10)  .  "ThTck 

800       Ans*=".2942" 

805   EDIT  "WHAT  IS  THE  THERMAL  CONDUCTIVITY  OF  THE  BORPD"  (WATTS-  M  -deg  K)",Ans* 

810   IF  Rns*="BACK_UP"  THEN  755 

815   Kb«RBS<VRL<Rns*  ' 

820   PRINT  LIN''  1  )  ,  TAB'.  10  •  ,  "Cond 

825   EDIT  "INPUT  fl  SHORT  <<50  CHARACTERS:'  WORD  DESCRIPTION  OF  THE  CIRCUIT  BOARD 

HERE  ",Pictlbl* 

830   IF  Pietlb1*="BRCK  UP"  THEN  300 

835   IF  LEN<Pict  lbl*X51  THEN  350 

840   GOSUB  Errin 

345   GOTO  325 

850   GOSUB  Board_pict 

855   GOSUB  Graphkey 

860   IF  Ans*="BACK__UP"  THEN  Bdpict 

865   GOSUB  Bd_cha 

370   IF  fins**"BRCK_UP"  THEN  850 

875   GOSUB  Case_f,'p 

830   IF  Rns*  =  "BflCi<_UP"  THEN  365 

835   GOSUB  Leads 

890   IF  RnsS*»"BRCK  UP"  THEN  375 

895   GOTO  435  !  END  OF  Bdp 

900  B o  ar  d_p  i  ■: »  ■=  :  GO  S  U  E  E  r  r  o  r 

905  Board~p1ct:  '  ThI--  SECTION  PLACES  ft  BLANK  BOARD  ON  THE  SCREEN 

910   ON  ERROR  GOTO  Board_picte 

915   PLOTTER  IS  "GRAPHICS" 

920   GRAPHICS 

925  N3CALE  O, 10 

930   CSIZE  3 

935   LORG  5 

940   MOVE  26, 139 

945   LABEL  " air  flow >  " 

950   Bdh=Bdh'Sca 

955   Bdl=Edl/Sca 

960   LORG  6 

965   LDIR  PI '2 

970   MOVE  175,70 

975   IF  Sca>l  THEN  LABEL  "THIS  PICTURE  IS  l/"'8cVAL*<Sca>&"  SIZE" 

930   IF  Sc*<l  THEN  LABEL  "THIS  PICTURE  IS  2X  SIZE" 

935   LORG  5 

990   LDIR  O 

995   MOVE  92-LEN<Pict1bl*)/2, 135 

10GO  LABEL  USING  "K";Pictlb1* 

1005  MOVE  0,0 

101O  DRAW  0,Bdh 

1015  DRAW  Bdl , Bdh 

102O  DRAW  Bdl ,0 

1025  DRAW  0,0 
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BOARDS  ..  DATA  INPUT  FOR  THEPMELEX 

1030  LINE  TYPE  3 

1035  Lr=Edl  N*r  !  LP  =  LENGTH  OF  EACH  REGION 

1040  Hr=Bdh-  N;  r  !  HP  =  HEIGHT  OF  EACH  REGION 

1045     FOR  1=1  TO  Nxr-1 

1050     MOVE  I*Lr.O 

1055     DRAW  I*Lr,Bdh 

1060   NEXT  I 

1065     FOR  1=1  TO  Nyr-1 

1070     MOVE  0, I*Hr 

1075     DRAW  Bdl , I*Hr 

1030   NEXT  I 

10S5  LINE  TYPE  1 

1090  LORG  5 

1095  Nr=0 

1100     FOR  J=Nvr  TO  1  STEP  -1 

1105       FOR  1  =  1  TO  N-;r 

1110       Nr=Nr+l 

1115       X<Nr>=Lr*. I-n+Lr  2 

1120       Y<Nr)=<J-l>*Hr+<l+Hr  2) 

1 125       MOVE  X'Nr  > , YCNr  ' 

1130       LfiEEL  USING  "K  "  ;  VRL*  •  Nr  ,• 

1135     NEXT  I 

1140   NEXT  J 

1145  Bdh»Bdh*Sca 

1150  Bdl»Bdl*Sca 

1155   WRIT  1508 

1160   EXIT  GRAPHICS 

1165  RETURN 

1170  Graphkeye:  GOSUB  Error 

1175  Graphkey: !  HERE  15  SUBROUTINE  TO  ALLOW  INPUT  OF  BOARD  DECRYPTION  FROM  KEYS 

DATA  INPUT 
1188  ON  ERROR  GOTO  Graphkeye 
1135  Ans**"Y" 

1190  S**"KEY  CODE  OVERLAY" 
1195  GOSUB  Pagehead 
1200  PRINT  "    This  section  allows  the  components  to  be  defined  and  drawn  on  th« 

screen" 
1205  PRINT  "using  the  keys  in  the  upper  right  ■;orri~r    of  the  keyboard.   The  pr*:>- 
ous  key" 

1210  PRINT  "deff initions  are  nor,  valid  while  entering  the  components  but  .Jill 
be" 
1215  PRINT  "during  later  portions  of  this  program.   As  the  flashing  cursor  moves 

to  each" 
1220  PRINT  "region  in  numerical  order  press  the  corespondi  r.g  key  and  wait  for  tF 
e  curso r  t o " 
1225  PRINT  "appear  at  the  next  region.   All  regions  must  be  defined  even  if  wiU 

k0=>EMPTY.  " 
1230  PRINT 

1235  PRINT  "  NOTE:   k  1  l->k  14  ni~r    to  chip  carriers  or  flat  pacts." 
1240  PRINT 

1245  PRINT  "     If  you  do  not  have  a  plastic  key  code  overlay  filled  in  (HP  part 
*7 120-61 64  •> ,  I  will  make  a  paper  one  for  you." 
1250  Anst="N" 

1255  INPUT  "DO  YOU  NEED  A  PAPER  OVERLAY  <N  or  Y>?",Ans* 
1260  IF  Rns*="BACK_UP"  THEN  RETURN 
1265  IF  UPC t    hns*[ 1 ,  1  ]>  =  "N"  THEN  1320 
1270  PRINTER  IS  0 

1275  PRINT  SPA  <  10.'  ,  "TEAR  OFF  FOR  COMPONENT  DEFF  INI  TION  KEY  CODE  OVERLAY  "  ,  L  I  N  (  2 
) 


98 


drtm  input  fop  thepnele: 


1230 
12S5 
BOVE 
1290 

T 

1295 

1300 

1305 

BELOW 

1310 

1315 

1320 

1325 

1330 

1335 

1340 

1345 

1350 

1355 

1360 

1365 

1370 

1375 

1380 

13S5 

1390 

1395 

1400 

1405 

1410 

1415 

1420 

1425 

1430 

14  35 

1440 

1445 

1450 

1455 

1460 

1465 

1470 

1475 

14S0 

1485 

14  90 

1495 

1500 

1505 

1510 

1515 

1520 

1525 

1530 

1535 

1540 

1545 

1550 


•mpv 
:ko) 


PRINT 
PRINT 
KEVS" 
PRINT  "TERR  HERE 


HOPIZ 
DIP14 


VERT     HORIZ 
D  I  P  1  4    D  I  P  1  6 


VERT 

D I P  li 


HORIZ    VERT 
DIP24    DIP24 


PLACE  A 


■THEN  PRE'! 


:on 


PRINT 
PAUSE 
PRINT 
KEYS 
PRINT 
PRINT 
PRINT 
'  HOT 
ON  KE 
ON  KE 
ON  KE 
ON  KE 
ON  KE 
ON  K  E 
ON  KE 
ON  KE 
ON  KE 
ON  KE 
ON  KE 
ON  KE 
ON  KE 
ON  KE 
ON  KE 
IF  01 
CRRFH 


IF 

IF 

PO 

EN 

RE 

NEX 

Enpty 

Itype 

GOTO 

Hdipl 

Npi  rv'. 

IF  It 

LABEL 

Wed) 

Le<I) 

IF  (I 

Te*p_ 

He<  I  .< 

Led  • 

IF  It 

Drayd 


LIN.'  1  ) 

"|  Ck3)   |  HORIZ  |  VERT   |   FLAT   |   FLAT  |   FLAT  |   FLAT  |     PLACE 

"|         !  DIP40  |  DIP40  |  16  PIN  |  24  PIN |  40  PIN ]  64  PIN  | " 

LIN<4> 
ER  IS  16 
E  THRT  THE  KEY*  IS  THE  CODE  FOR  TYPE  OF  ELEMENT  STORED  IN  I  type  a.' 

Y  *0  GOTO  E'uf.  t 

Y  #1  GOTO  Hdipl4 

Y  #2  GOTO  Vdi pi  4 

Y  »3  GOTO  Hdipl6 

Y  #4  GOTO  VdiplS 

Y  #5  GOTO  Hdip24 

Y  *6  GOTO  Vdip24 

Y  #7  GOTO  14  35 

Y  #3  GOTO  14  35 

Y  #9  GOTO  Hdip40 

Y  #10  GOTO  Vdip40 

Y  #11  GOTO  Fl  it  16 

Y  #12  GOTO  Fl  it  T4 

Y  #13  GOTO  Flit  40 

Y  #14  GOTO  Fl  it  64 
de*PI  THEN  RETURH 
I C  3 

FOR  1=1  TO  Nreg 

IF  I type ( I )»1Q  THEN  Ne  tr 

MOVE  X<I >,Y<I) 

POINTER  X<I),Y<I  »,2 
HIS  IS  AN  ENDLESS  LOOP  TO  ALLOU  FOR  USE  OF  THE  PRIORITY  INTTERUPT  KEYS 

GOTO  L  o  o  p 
:  IF  Reg  cha<>0  THEN  RETURN 
lONreg  THEN  1475 
OldeOPI  THEN  WRIT  100*NregM.5 
INTER  0.0,2 
IT  GPRFHICS 
TURN 
T  I 

:G0SUB  Er&senum 
(  I  )=Npin<  I  >=We<  I  ,'=Le<  I  >=0 
Ha  tr 
4 : GOSUB  Er asenum 

I  >  =  14 

ype<I>=2  THEN  LDIP  PI- 2 

USING  "K";  "DIP  14" 
=We=.25*25.4       !  ACTUAL  WIDTH  IN  r„r„ 
=Le=. 725*25. 4      !  ACTUAL  LENGTH  IN  mm 
type<I)=l)  OR  <Itype<I)=0)  THEN  1545 

dim=We<I)  !  HERE  LENGTH  AND  WIDTH  TRANSPOSED  SUCH  THRT  LENGTH 

*We=Le<I)  !  IS  DEFINED  TO  BE  IN  THE  DIRECTION  OF  AIR-FLOW 

=L*=T~mp_di m 

ypeU>*9  THEN  Itype<I)=l 
ip:!  THIS  SECTION  WILL  DRAW  RLL  THE  DIP  FIGURES  ON  THE  CRT  IF  Npin<23 


BOARDS  ..  DATA  IMFUT  FOR  THERMELE:: 

1555  L»«L*/'Sca 

1560  Me*We  I  ;  i 

1565  MOVE  ;;i!)-L{.2,'rI'-N-  2 

1570  DRAW  X<I>-Le  2.  Y*  I  ■ +U-:  2 

1575   draw  xd^Le  :,iq>nii   : 

1580  DRAW  :;<I'+l_e  2,Y<1>-Ws  2 

1535  DRAW  X<I>-Le  2,Y<D-U«  2 

1590  LDIR  9 

1593  GOTO  Nextr  !  End  of  Draw_dip 

1600  Vdipl4:  Ityp«<  I  >*2 

1605  GOTO  Hdipl4 

1610  Hdipl6:G0SUB  Erasenum 

1515  Npin<I>=16 

1620  IF  Itype<I>*4  THEN  LDIP  PI  2 

1625  LABEL  USING  "K" J "DIP16" 

1630  Ue<I)*Ws=.25*25.4        !  WIDTH  IN  mm 

1635  Le<I>=Le=.325*25.4       !  LENGTH  IN  mm 

1640  IF  ItypedX  4  THEN  Itype<I>=3 

1645  IF  ItypeCI)   4  THEN  Drawdip 

1650  Temp_dim=Me<  I ) 

1655  We(  l">*We*Le<  I ) 

1660  Le<I)*Le=Temp  di« 

1665  GOTO  Drawdip  !  End  of  HdiplS 

1670  Vdipl6: Itype<  I  '=4 

1675  GOTO  Hdiplo 

1630  Hdip24:G0SUB  Erasenuru 

1635  Npin< I, =24 

1690  IF  Itype<I>*6  THEN  LDIP  PI -2 

1695  LABEL  USING  "K";"DIP  24" 

1700  He<I)=We*.55*25.4 

1705  L»CI)=Le*1.25*25.4 

1710  IF  Ityp*a><>6  THEN  Ityp*<I>=5 

1715  IF  Uype<n<>6  THEN  Drawdip 

1720  T*»p  di m=We>  I 

1725  We<l7=We=L*<I> 

17  30  Le<I>=L*»Temp  dim 

1735  GOTO  Drawdip  !  End  of  Hdip24 

1740  Vdip24: It ypc (I)=6 

1745  GOTO  Hdip24 

1758  Hdi'p40:IF  ■' I  MOD  Nxr=0)  OP  <Itype< 1+1 >=10>  AND  <Reg_cha=0>  THEN  Loop 

1755  GOSUB  Er asenum 

1760  I=Ih=I+l 

1765  COSUB  Erasenum 

1770  1=1-1 

1775    Npin<I>=Npin<Ih)=20  !     THIS    IS    FOR    EACH    HALF    OF    THE    40    PINS 

1780    Itype<I>  =  Itypeah>=9 

1735    MOVE    X<I)+Lr/2,Y<I) 

1790  LABEL  USING  "K":"BIP40" 

1795  We< I  »=U« < In >=We=. 6*25. 4 

18O0    L*<I)=Le<Ih>=L*=25.4  !     THIS     IS    FOP    EACH    HALF    IN    EACH    REGION 

1805    We=Wfi  'Sc i 

1310    Le  =  L€  'Sc i 

1815    MOVE    X<I>+Lr/2-Le, Y<I>+We/2 

1320    DRAW    X<I>+Lr/2-Le,Y<n-We/-2 

1325    DRAW    X<0+Lp/2+Le,Y<I>-We/2 

1830    DRAW    K<I>+Lr/2+L«,Y<I)+We/2 

1335    DRAW    X<I)+Lr    2-Le, Y< I )+We/2 

1848    I=Ih 
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BOARDS  ..  DATA  INPUT  FOP  THERMELE:: 

1345  GOTO  Nextr 

1850  Vdip4Q:lF  I+Nxr>Nreg  THEN  1425 

1855  GOSUB  Erasenum 

I860  I=Ih=I>Nxr 

18-55  GOSUB  Erasenum 

1S70  I=I-N<r 

1S75  Mpin<I>=Npin<Ih^*2@     !  THIS  IS  FOR  THE  HALF  OF  THE  40  PINS  IN  EACH  REGION 

1880  [typ«( I >=Ityp*< Ih;=10 

1335  LDIR  PI  2 

1890  HOVE  X<  I  ),Y<I)-Hr/-2 

1895  LABEL  USING  ,,K";"DIP40" 

1900  LDIR  0 

1905  Ue< I >*W*< Ih)=We*. S*25. 4 

1910  Le<I)=LeCIh>*Le=25.4     !  THIS  IS  FOP  HALF 

1915  W*  =  W*.Sca 

1920  Le=Le/Sca 

1925  MOVE  X<I>+We  2,  Y<  I  >-Hr-'2+Le 

1930  DRAW  X<  I  >-We  2,  Y<  I  >-Hr  2  +  Le 

1935  DRAW  X<I)-We/2,Y<I)-Hr/-2-Le 

1940  DRAW  X  I  j-fW*  2.  •  ■  I  ;>-Hr-  2-L-E 

1945  DRAW  X<I>+We  2,  Y<  I  )-Hr>'2+Le 

1950  LDIR  0 

1955  Te»p  dim=Le<I> IHERE  TRANSPOSE  Li    AND  Ws  TO  MAINTAIN  LENGTH  DEFF  WITH  AIP 

1960  LeClT=Le< Ih)*We<I) 

1965  U«< I )=Ws < Ih)=Temp_di m 

1970  GOTO  N«xtr 

1975  Flat 16: I type" 1=11 

1980  GOSUB  Erasenum 

1935  Npintl)»16 

1990  LABEL  USING  "K": "16" 

1995  Le«W*=Me<I>=Le<I)=4.57 

2000  GOTO  Drawdip 

2005  Flat  24:  1 1  ype-:  I  .' =  1  2 

2018  GOSUB  Eriienum 

2015  Hpin<I>=24 

2020  LABEL  USING  "K"; "24" 

2025  L««We»We<I)=Le<I)=7.75 

20  JO  GOTO  Drawdip 

2035  Flat  40: Itype<I>*13 

2040  GOSUB  Erasenum 

2045  Npin<I)=40 

2050  LABEL  USING  "  K  "  ; "40" 

2055  L«=Ue=Ue< I ^=Le* I  >  =  12.  19 

2060  GOTO  Draudip 

2065  Fl  it  64:  I  type'  I  .'  =  14 

2070  G01UB  Erasenum 

2075  Npin<I>=64 

2080  LABEL  USING  "K";"FLAT64" 

2035  Le=We=We<I>=Le< I)=18.29 

2090  GOTO  Drawdip  !  END  OF  COMPONENT  DEFFINITIOHS 

2095  ! 

2100  Erasenum: •  HERE  WE  MOVE  THE  REGION  NUMBER  TO  THE  UPPER  LEFT  CORNER 

2105  CSIZE  3 

2110  MOVE  X(I),Y(I) 

2115  PEN  -1 

2120  LABEL  USING  "K  ";  I 

2125  CSIZE  2.5'Sca 

2130    IF    Reg    cha=0    THEN    PEN     1 
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BOARDS  ..  DATA  IMPur  FOR  ThFRMELE:; 

2135  MOVE  X<I>-.4*Lr,Y<I)+.4*Hr 

2140  LABEL  USING  "K"J I 

2145  MOVE  X< I),Y<I) 

2150  CSIZE  3.4/Sca 

2155  RETURN  \     END  OF  ERfiSENUM 

2160  ' 

2165  Bd_chae: GOSUB  Error 

21?0  Bd"~cha: !  THIS  ROUTINE  ALLOWS  FOR   CHANGE  OF  TYPE  FOR  SPECIFIED  COMPONENTS 

2175  ON  ERROR  GOTO  Bd_chae 

2  ISO  3*  =  " CHANGE  COMPONENTS" 

2185  GOSUB  Pagehead 

2190  PRINT  "     This  section  allows  corrections  to  the  components  on  the  circuit 

board. " 
2195  PRINT  "However;  you  may  not  change  the  circuit  board  itself.   You  will  fir; 
t  remove" 

2200  PRINT  "ill  unwanted  components  Cup  to  10)  by  entering  the  region  number  on 
the  pi  c ture. " 
2205  PRINT  "When  removing  any  component  which  requires  two  spaces,  use  the  lower 

regi  on  num. " 
2210  PRINT  "Entering  a  0  or  pressing  CONT  with  no  entry  wi 11  allow  program  to 
cont  i  nue" 

2215  PRINT  "with  the  section  that  allows  deffinition  of  components  in  all  empty 
spaces.  " 
2220  PRINT 

2225  PRINT  "Remember  to  allow  space  for  the  larger  component  3  that  require  two  r 
eg  ions.  " 
2220  ltt=0 

2235  Again:  Itt=Itt+l 
2240  Ans*="0" 

2245  INPUT  "ENTER  THE  REGION  NUMBER  TO  REMOVE  COMPONENT  < DEFAULT* >0  =>N0NE> . " , fit 
s* 

2250  IF  Ans*="BACKJJP"  "HEN  RETURN 
2255  R«g_cha*VAL<flns*) 
2260  IF  Reg_cha=0  THEN  Rdd_comp 

2265  IF  CReg_cha>0)  AND  <Reg_cha<Nreg+l )  THEN  2230 
2270  GOSUB  Err  in 
2275  GOTO  Again 
2280  Reg_cha' Itt)=Reg_cha 
22S5  GOSUB  Erase_comp 
2290  Reg_cha=0 
2295  IF  Ftt<  =  10  THEN  Again 
2300  S*="MAXIMUM  NUMBER  OF  REGIONS" 
2305  GOSUB  Pagehs ad 

2310  PRINT  "    There  are  10  changes  to  accomplish  at  this  point  and  to  avoid  pre 
blems  you" 

2315  PRINT  "must  redefine  the  beard." 
2320  PRINT  LINO), SPR<15),  "Press  CONT  " 
2325  GOSUB  Add_comp 
2330  GOTO  Bd  cha 

2SS5  Erase_con,p:  !  THIS  SECTION  REMOVE  UNWANTED  COMPONENT  TO  ALLOW  CHANGE 
2340  GRAPHICS 
2345  PEN  -1 
2350  I=Peg_cha 

2355  IF  Itype<Reg_cha)=8  THEN  2365 

2360  ON  1 1 ype  <  I  >  GOSUB  Hd 1 p 1 4 . Vd -pi4. Hd 1 p 1 6 . Vd  i  plS, Hdi  p24, V d  i  p24, Err  o 1 de , Err  ol ds 
,  Hd  1  p40  ,  Vd  1  p40  ,  F  1  at  1  6  ,  F  1  at  24  ,  F  1  at  40  ,  F  1  at  6  4 
2365  PEN  1 
2370  GOSUB  Empty 
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BOARDS  ..  DATA  [NPUT  FOR  thermele:: 

2375  IF  I<2  THEM  2409 

2380  IF  <Itype<I-l>=9)  AND  CItype<I)=9)  THEN  2399 

2335  GOTO  2400 

2390  1=1-1 

2395  GOSUB  Empty 

2400  IF  KHxr+1  THEil  2420 

2405  IF  Itype<I-Nxr)Ol0  THEN  2428 

24  10  I*I-Nxr 

2415  GOSUB  Empty 

2420  WAIT  500 

2425  EXIT  GRAPHICS 

2430  RETURN 

2435  fldd_Comp:  '  THIS  SECTION  ALLOWS  ADDITION  OF  COMPONENT:'.  TO  EMPTY  SPACE? 

2440  IF  Itt<2  THEN  2475 

2445  IF  01de=PI  THEN  GOSUB  Graphkey 

2450  FOR  1=1  TO  Nreg 

2455  IF  Itype<I)O0  THEN  2470 

2460  POINTER  X<I),Y<I>,2 

2465  GOTO  2460 

2470  NEXT  I 

2475  FOR  1=1  TO  15 

2430  OFF  KEY  *  I 

2435  NEXT  I 

2490  EXIT  GRAPHICS 

2495  RETURN  !  END  OF  ADDjCOMP 

2500  ! 

2505  Case_type: GOSUB  Error 

2510  Case  typ!  !  THIS  SECTION  ALLOWS  THE  USER  TO  SPECIFY  THE  CASE  TYPE  FOR  DIPS 

2515  ON  ERROR  GOTO  Case  type 

2520  IF  v01de=PI>  ANT  <Ttt<2)  THEN  GOTO  Case_cha 

2525  !  BUT  FIRST  MUST  TELL  THE  USER  TO  CHANGE  THE  KEY-CODE  OVERLAY 

2530  3*= "CHANGE  OF  THE  SPECIAL  FUNCTION  KEYS" 

2535  GOSUB  Pagehead 

2540  PRINT  LIN<10),SPA<10), "REMOVE  THE  KEY  CODE  OVERLAY  FOP  COMPONENT  INPUT" 

2545  PRINT  LIN<1>,SPA<10), "REPLACE  THE  GENERAL  PURPOSE  KEY  CODE  OVERLAY" 

2550  DISP  SPAC25),  "PRESS  CONT  WHEN  READ'.'" 

2555  PAUSE 

2560  S*="TYPE  OF  DIP  PACKAGE" 

2565  GOSUB  Pagehead 

2570  PRINT  "     Dual  inline  packages  CDIPs)  generally  cor,,-:-  in  either  an  injectio 

n  mo 1 de d " 

2575  PRINT  "plastic  case  or  a  ceramic  sandwich  case.   The  style  of  case  construe 

t i on  affects" 

2530  PRINT  "the  heat  transfer  and  must  be  known  for  the  thermal  model." 

2535  PRINT 

2590  PRINT  "The  methods  of  specifying  the  types  of  cases  are  listed  below." 

2595  PRINT  TAB< 15); "THESE  APPLY  TO  THE  DIPs  ONLY"  . 

2600  PRINT  LIN' .2) 

2605  PRINT  TAB<10>;"1.  HLL  CERAMIC" 

2610  PRINT  TAB<10>;"2.  ALL  PLASTIC" 

2615  PRINT  TAB<10);"3.  MAJORITY  CERAMIC  'USER  SPECIFY  WHICH  ARE  PLASTIC)" 

2620  PRINT  TAB<10>;"4.  MAJORITY  PLASTIC  (USER  SPECIFY  WHICH  ARE  CERAMIC)" 

2625  PRINT  TAB<10>;"5.  USER  SPECIFY  CASE  STYLE  FOR  EACH  hi?    COMPONENT" 

2630  AnsS="  1" 

2635  INPUT  "ENTER  YOUR  CHOICE  <  1  .  2 ,  3  ,  4  ,  5  >  "  ,  Ar,i  X 

2640  IF  Ans*<>"BACKJJP"  THEN  2S70 

2645  PRINT  LIHC20),SPA<20), "UNABLE  TO  BACK  UP  WITHOUT  LOSING  THE  BOARD" , LIN<5) 

2650  BEEP 
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2655    DISP    "Preii    CONT    when    ready    to    continue" 

2660    PAUSE 

2665  GOTO  Case_typ 

2670  Rns=VRL<Ans*> 

2675  IF  (Ans<l)  OR  <flns>3)  THEN  Case_type 

2630  ON  Ans  GOSUB  Cere,Plas, Mcere,Mpl as, Any 

2635  Case_cha:S**"CORRECTIONS  TO  CASE  STYLE  FOP  DIP  PACKAGES" 

2690  GOSUB  Pagehead 

2695  PRINT  "     On  the  screen  below  is  a  line  that  represents  the  case  type  for 

all  regions." 

2700  PRINT  "Those  regions  that  contain  a  plastic  cased  component  are  shown  in  in 

werse  y  i  d  i  o  " 

2705  PRINT  "<1)  while  ill  others  are  shown  in  normal  video  <  1 )  .  " 

2710  PRINT 

2715  PRINT  "To  make  corrections  enter  the  region  #  that  is  incorrect, if  no  chang 

e s  no  ent  r y .  " ,  L I N  <  1 ) 

2720      FOR  1=1  TO  Nreg 

2725      IF  I  NOD  10=1  THEN  PRINT  SPRUO); 

2730      IF  <CaseXC I,  I]  =  "0"  <  AND  <I>  =  19)  THEN  PRINT  I;SPA<1>; 

2735      IF  <Case*CI,  I]  =  "0  '  ■  AND  <:  I  1  0  .'  THEN  PRINT  SPA<  1  )  ;  I ;  SPA<  1)  ', 

2740      IF  (CaseSCI,  I3  =  "l")  AND  ■  I  >  =  1  0  :■  THEN  PRINT  CHR*  <  129  )  ;  I ;  CHRS  <  123)  ;  SPA  <  1  ) 

2745      IF  (CaseSC  I,  I  ]  =  "  1"  >     AND  •'  I  <  1  0  >  THEN  PRINT  C  HR.T  •:  1  29  ,  ;  "  "  ;  I  ;  C  HRf  '.  1  28  >  ;  SPA 

ci>; 

2750      IF  I  MOD  10=0  THEN  PRINT  LINC1) 

2755      NEXT  I 

2760      PRINT 

2765  Ans»="NULL" 

2770  INPUT  "ENTER  REGION  *  TO  CHANGE  OP  PRESS  CONT  WITH  NO  ENTRY  IF  0K",Ans* 

2775  IF  (flns**"NULL")  OR  CAns*="BACKJJP">  THEN  RETURN 

2730  Pegl*VflL« Ans* > 

2735  IF  (Regl>0)  AND  <Regl<=Nreg>  THEN  2300 

2790  GOSUB  Err in 

2795  GOTO  Case_cha 

23O0  CaseSCRegl, Regl 1=VAL*< INK < VAL< Case* C Reg  1 , Regl3>+1 >  MOD  2) > 

2805  GOTO  Case_cha  !  END  OF  Case_cha 

2810  ' 

2815  Cere: '  ALL  CERAMIC  CASES 

2320  FOR  1=1  TO  Nreg 

2825  CaseJC I, I]="0" ' 

2830  NEXT  I 

2835  RETURN  I  END  OF  CERE 

2340  ! 

2845  PI  as:'  ALL  PLASTIC  CASES 

2850  FOP  1=1  TO  Nreg 

2855  CaseSCI , I ]="1" 

2860  NEXT  I 

2865  RETURN  ! END  OF  PI  as 

2870  ! 

2875  Mcere:   S*="MAJORITY  CERAMIC  ...  SPECIFY  WHICH  ARE  PLASTIC" 

2830  GOSUB  Page  he  ad 

2335  FOR  1=1  TO  Nreg 

2390  Cases  CI , I]="0" 

2895  NEXT  I 

2900  Ans*="NULL" 

2905  INPUT  "ENTER  THE  NUMBER  OF  PLASTIC  COMPONENTS" , Ans* 

2910  IF  Ans*»"BACi  JP"  THEN  Case_typ 

2915  Rns*VAL<flns*> 
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292(3  IF  CAns>0)  FIND  <  AnssNr  eg+ 1  )  THEN  2935 

2925  GOSUE  Err  in 

2930  GOTO  Mcere 

2935  IF  Ans=Nreg  THEN  P 1  as 

2940      FOR  1=1  TO  Ans 

2945      Ans*="NULL" 

2950      INPUT  "ENTER  THE  REGION  NUMBER  THAT  CONTAINS  A  PLASTIC  COMPONENT ", An  =  $ 

2955      IF  Ans*="BfiCK  UP"  THEN  Mcere 

2960     Regl=VAL<Ans*1> 

2965      IF  (Regl>0>  AND  (Regl  Nreg+1)  THEN  29:30 

2970      GOSUB  Errin 

2975      GOTO  2950 

2930     Case*CRegl,Regl]="l" 

2935      NEXT  I 

2990  RETURN  !  END  OF  Mcere 

2995  ! 

3000  Mplas:  3*= "MAJOR  I  TV  PLASTIC  ...  SPECIFY  WHICH  APE  CERAMIC" 

3005  GOSUB  Pagehead 

3010      FOR  1*1  TO  Nreg 

3015      CaseSCI, I]="l" 

3020      NEXT  I 

3025  fln*$a"NULL" 

3030  INPUT  "ENTER  THE  NUMBER  OF  CERAMIC  COMPONENTS ", Ans $ 

3035  IF  fins*«" BACKUP"  THEN  Case  typ 

3040  An**INT<VAL<Ans*)) 

3045  IF  (fins    0)  AND  (Ans<Nreg+l>  THEN  30-50 

3050  GOSUB  Err  in 

3055  GOTO  Mpl is 

3060  IF  Ans*Nreg  THEN  Cere 

3065      FOR  1=1  TO  Ans 

3070      Rns**"NULL" 

3075      INPUT  "ENTER  THE  REGION  NUMBER  THAT  CONTAINS  A  CERAMIC  COMPONENT ". Ani f 

3030      IF  Ans*="BACKJJP"  THEN  Mplas 

3035      IF  CR«gl>0)  AND  CRegKNreg+l)  THEN  3100 

3090      GOSUB  Errin 

3095      GOTO  3075 

31O0     Case*tRegl,Regl]=u0" 

3105      NEXT  I 

3110  RETURN  !  END  OF  Mplas 

3115  ! 

3120  Any: i  EACH  REGION  MUST  HAVE  IT'S  DIP  SPECIFIED 

3125      FOR  1=1  TO  Nreg 

3130      IF  CltypeCIXl)  OP  <  I  typed )  >10)  THEN  3160 

3135      DISP  "IDENTIFY  CASE  TYPE  FOP  REGION  # " ; I ; " ( ENTER  1  FOR  PLASTIC  AND  0  FO 

R  OTHER)"J 

3140     INPUT  Case*[ 1,1] 

3145      IF  CCase#EI, n="l")  OR  (Case*C I , I ]="0")  THEN  3160 

3150      GOSUB  Errin 

3155      GOTO  3135 

3160    -  NEXT  I 

3165  RETURN  ! END  OF  Any 

3170  ! 

3175  Mhat_now: '  THIS  SECTION  IS  THE  FINAL  SECTION 

3130  S*="WHAT  NOW  ?" 

3185  GOSUB  Pagehead 

3190  PRINT  "    You  hive  completed  one  cycle  through  BOARDS  and  have  the  follow 

i  ng  opt i ens: " 

3135  PRINT  LINa>,TAB<10>,  "1.  PERFORM  THERMAL  ANALYSIS  ON  "jName*;"  USING  THERM 

L.  " 
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3260  PRINT  LIN<1),TAB<10), "2.  INPUT  ANOTHER  CIRCUIT  EOAPD  DESCRIPTION  FROM  KEYED 

ARD." 

3205  PRINT  LIN<1 ) , TAB( 10), "3.  INPUT  ANOTHER  CIRCUIT  BOARD  DESCRIPTION  FROM  MASS 

STORAGE. " 

3210  PRINT  LINM  ),TAEUO,',  "4.  TERMINATE." 

3215  ! 

3220  Ar»s*  =  "l" 

3225  INPUT  "YOUR  CHOICE  ( 1 , 2, 3, 4> " , fins* 

323©  IF  Ans*=MBflCK_UP"  THEN  RETURN 

3235  ftns*ABS< INT(VflL(flnst) > ) 

32-40  IF  Ans<  5  THEN  3255 

3245  GOSUB  Err  in 

3250  GOTO  3220 

3255  ON  fins  GOTO  Therml , Redo, Redo, Ends 

3260  Therml:DISP  "WORKING  LOADING  THERML  TO  PERFORM  ANALYSIS  OF  ";Name$ 

3265  LOAD  "THERML", 35 

3270  STOP 

3275  ' 

3280  Redo:  Hap«Ans 

3285  RETURN 

3280  ! 

3285  EnderPRINT  PAGE.  LIN-  20''  ,  SPA  ■  10  >,  "NORMAL  TERMINATION  " 

3300  DI3P  "WORKING  LOADING  STANDARD  KEY  DEFF INIT IONS" 

3305  GCLEAR 

3310  LOAD  KEY  "STDKEY" 

3315  END  !  END  OF  What_r.ow 

3320  ! 

3325  Pagehead: !  THIS  SECTION  PLACES  DESIRED  HEADING  ON  A  BLANK  CRT 

3330  PRINTER  IS  16 

3335  PRINT  PAGE,TAB(34-LEN<S*>  2s"***  M ; CHR*< 132> ; S*J CHR*< 123) ; "  ***",LIrK2> 

3340  RETURN  !  END  OF  PAGEHEAD 

3345  ! 

3350  Errln:!  THIS  SECTION  ALERTS  THE  USER  TO  AN  ATTEMPT  TO  INPUT  BAD  DATA 

3355  BEEP 

3360  DISP  "INPUT  OUT  OF  RANGE TRY  AGAIN" 

3365  WAIT  1500 

3370  EEEP 

3375  WAIT  1000 

33S9  BEEP 

3385  RETURN 

3380  !  END  OF  Err in 

3385  i 

3400  Leads  e: GOSUB  Error 

3405  Lead.-:!  THIS  SECTION  ALLOWS  THE  INPUT  OF  THE  NUMBER  OF  LEADS  BETWEEN  REGION 

S 

3410  ON  ERROR  GOTO  L^  adse 

3415  S**"ELECTRICAL  LEADS  OR  OTHER  CONDUCTION  PATHS" 

3420  GOSUB  Pagehead 

3425  PPINT  "     A11  circuit,  boards  have  thermally  conductive  materials  applied  t 

o  their  " 

3430  PRINT  "surface  either  in  the  form  of  electrical  leads  or  as  thermal  conduct 

i  on  paths . " 

3435  PRINT  "Heat  will  be  transferred  through  these  materials  between  the  regions 

3440  PRINT  "if  there  are  conduction  paths,  (rails)  present,  they  will  haue  a  much 

1 arger " 
3445  PRINT  "effect  than  the  electrical  leads  and  the  electrical  leads  will  be  ne 
g  1  e  c  t  e  d  " 
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3450  PRINT  "in  the  thermal  model." 

3455  Ans*="N" 

3460  INPUT  "ARE  THERE  METAL  CONDUCTION  PATHS  PRESENT  -,  N  *r    Y)" ,  fins* 

3465  IF  Ans**"BACK  UP"  THEN  RETURN 

3478  IF  UPC*(Ans*[T, 1])="N"  THEN  Leads_size 

3475  IF  UPC*(Ans*Cl, i])»"Y"  THEN  3490 

3480  GOSUB  Err  in 

3485  GOTO  Leads 

3490  S*  =  " CONDUCT  I  ON  RAILS" 

3495  Al*="CPU's" 

3500  GOSUB  Pagehead 

3505  PRINT  "     Therru.il  conduction  rails  are  Much  wider  and  thicker  than  the  ele 

c  t  r  i  c  a  1  " 

3510  PRINT  "leads  plated  onto  the  circuit  board.   The  thickness  is  generally  uni 

form,  " 

3515  PRINT  "however;  the  rails  are  typically  not  of  uniform  width.   It  is  theref 

ore" 

3520  PRINT  "necessary  to  define  a  concept  of  Conduction  Path  Units  (CPU's).   The 

rail  width" 
3525  PRINT  "will  be  modeled  in  terms  of  an  i nt  eger  number  of  CPU's.   You  will 
be  asked  to" 

3538  PRINT  "input  the  thickness  (mm)  of  a  CPU  and  the  width  will  automatically  d 
efault  to" 

3535  PRINT  ".1  mm  such  that  a  rail  of  width  1.7  mm  can  be  modeled  as  17  CPU's." 
3540  PRINT 

3545  PRINT  "     For  your  convenience  the  thermal  conductivities  (watts^M-C)  of    t 
hree  all oys" 

3550  PRINT  "commonly  used  as  material  for  conduction  rails  are  given  below:", LIN 
(i) 

3555  PRINT  "A!  5052=  133.2  W- M-C  AL  5101  =  215.7  W-M-C  Cu  113  ■  339 

.8  W-'M-C" 

3560  PRINT 

3565  01de=0 

3570  Ans*="l" 

3575  INPUT  "ENTER  THE  THICKNESS  OF  THE  THERMAL  CONDUCTION  PATHS  or  CPU's  (mm)", A 

ns$ 

3580    IF    Ans**"BACK_UP"     THEN    Leads 

3535    T1*ABS<  '/AL'  Ans*  i  i 

3590  PRINT  "THICKNESS  OF  CPUs  =  ";Tl;"mm" 

3595  Ans-r=". 1" 

3680  !  INPUT  "ENTER  THE  WIDTH  OF  THE  CPU's  (mm)",Ans* 

3605  IF  AnsJ="EACK_UP"  THEN  3570 

3610  HI =ABS<VAL(flns«) ) 

3615  fll=-Wl*Tl   i *****#**#  NOTE  THAT  THIS  AREA  IS  NEGATIVE  FOP  CPU's 

3620  PRINT  "WIDTH  OF  CPU's      =  ";Ul;"mm       AREA  FOP  HEAT  TRANSFER  PER  CPU  =" 

;  -A  1  ;  "  turn "2  " 

3625  Ans*="133.2" 

3630  INPUT  "ENTER  THE  THERMAL  CONDUCTIVITY  OF  THE  CONDUCTION  PATHS  ( W/M-C) " , Ans* 

3635  IF  Ans**"BACK_UP"  THEN  3570 

3640  Kl=ABS<VAL<Ans$) ) 

3645  PRINT  "THERMAL  CONDUCTIVITY  = " ; K 1 ; "W at t s-  M-De g  C" 

3650  UAIT  1500 

3655  GOTO  Leads  enter  !  END  OF  CPU 

3660  ! 

3665  Leads_sire:S*="NUMBER  OF  LEADS" 

3670  GOSUB  Page  he  ad 

3675  Ali="LEADS" 

3630    PRINT    "  This    section    of    the    program    allows    the    entry    of    the    size    and    con 

duct  wily    of" 
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3635  PRINT  "leads  (traces)  on    the  surface  of  the  circuit  board.   These  act  as  bo 

th" 

3690  PRINT  "electrical  and  thermal  conductors  between  the  regions.   When  enter  in 

g  the  width" 

3695  PRINT  "figure  an  average  width  for  all  the  leads." 

3700  PRINT  LINO.',''    The  thermal  conductivity  will  default  to  pure  copper  '334 

Uatts-'M-Beg  C>" 

3705  PRINT  "with  no  entry." 

3710   01de=0 

3715   Rns*»".033" 

3720   INPUT  "ENTER  THE  THICKNESS  OF  THE  LEAD?  (TRACES)  ON  THE  CIRCUIT  BOARD  (mm) 

", AnsS 

3725   IF  Ans*="BACK_UP"  THEN  Leads 

3730   Tl«flBS' VAL<Ans* ' i 

3735   PRINT  LIN< 1 ) , "LEAD  THICKNESS       =";Tl;"mm" 

3740   Ans-r="l" 

3745   INPUT  "ENTER  THE  AVERSE  WIDTH  OF  THE  LEADS  '.TRACES:'  ON  THE  CIRCUIT  BOAR 

D  <mm) " , Ans* 

3750   IF  Ans*="BACK  UP"  THEN  3715 

3755   Ml=RBS<VAL<Ans*>) 

3760   A1*W1*T1 

3765   PRINT  "AVG  WIDTH  OF  LEADS   =  ";Wl;"mm       AREA  FOR  HEAT  TRANSFER  PER  LEAD 

=  "  ;  Al  ;  "mm-x2  " 
3770   Ansf="384" 

3775   INPUT  "ENTER  THE  CONDUCTIVITY  OF  THE  TRACES  'DEFAULT  COPPER  =  334W  at  t.  s-'  M-d 
egK)", Ans* 

3780   Kl =  ABS<  VAL< Ans* >  ) 

3735   PRINT  "CONDUCTIVITY  OF  THE  LEADS  =  "  ;  h' 1  ;  "  W at t = - N-C " 
3790   WAIT  1500 

3795  Leads_*nter:  S*="NUMBERS  OF  "iAlfi"  BETWEEN  REGIONS" 
3300  GOSUE  Pagene  ad 

3805  PRINT  "     This  section  of  the  program  allows  the  entry  of  the  number  of  th 
e  ";A1* 

3810  PRINT  "that  cross  each  of  the  internal  region  boundaries.   It  these  are  on 
both  side-" 

3815  PRINT  "of  the  board  add  both  numbers  together.   The  flashing  cross  will  mov 
e  to  the" 
3820  PRINT  "appropriate  location  on  the  screen  but  the  graphics  picture  will  not 

remai  n. " 
3825  PRINT  "If  a  hard  copy  of  the  graphics  is  needed  for  a  guide,  recall  that  k 3 

will" 
3830  PRINT  "provide  one.   The  picture  will  return  to  the  screen  while  the  number 

of  ";Aif 
3835  PRINT 

be  entered. " 
3340  PRINT  LIN.ri 
ING  ENTERED' " 

3345  PRINT  LINU),  "There  will  be  correction  opportunities  later." 
3858  DISP  "PRESS  CONT  WHEN  READY  TO  START  WITH  REGION  #1" 
3855  PAUSE 
3360   GOTO  Leads_in 

3365  Leads  label:   '  HERE  THE  NUMBERS  OF  LEADS  OR  CPU'S  IS  WRITTEN  ON  THE  SCREEN 
3870  CSIZE  2.5 
3875  GRAPHICS 
3880  LABEL  USING  "k";Nl 
3885  PEN  -1 
3390  LDIP  -PI 
3895  LABEL  USING  "K":N1 
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3988  PEN  1 

3905  LDIR  9 

3910  IF  Olde-  PI  THEN  WRIT  500 

3915  RETURN 

3920  Leads_old:CSIZE  2.5 

3925  S#*"NUM£ERS  OF  LEADS" 

3930  A1J="LEADS" 

3935  IF  AKO  THEN  Al*»"CPU's" 

3940  GOSUB  Pagehead 

3945  PRINT  "     This  section  of  the  program  displays  the  numbers  of  leads  cros: 

ng    ■ 

3950  PRINT  "boundaries  of  regions  by    placing  numbers  on  the  sides  of  the  regioi 

that " 
3955  PRINT  "represent  the  informations  stored  in  ";Name*;".   You  will  be  allow 

to  nak €  " 
3960  PRINT  "corrections  to  that  information  and  re-store  in  on  mass  storage  if 
equi red. " 
3965  PRINT 

3970  i  PRINT  "It  you  do  not  desire  to  view  the  data  concerning  number  of  leads 
enter  any  " 

3975  !  PRINT  "number  before  you  press  C0NT  . " 
3988  DISP  "PRESS  CONT  WHEN  READY  TO  VIEW  DATA  " 
39S5  Leads  in:  FOP  1=1  TO  Nreg 
3990   J=l   !  HERE  I?  BOTTOM  OF  THE  REGION 
3995   GRAPHICS 

4000   IF  I+Nxr<*Nreg  THEN  4025 

4005   Nl  < J,  I  '=N1 =0    i  HERE  IS  BOTTOM  OF  THE  BOARD  WHERE  CONNECTORS  WILL  GO 
4010   MOVE  X<I),Y<I>-.45*Hr 
4015   GOSUB  Leads  1 abel 
4020   GOTO  4165 

4025   IF  (Itype<IX>10)  OR  < I type< I+Hxr ><>10)  THEN  4  100 
4030   MOVE  X<I)+.32*Lr,Y<I)-.45*Hr 
4035   IF  01de< >PI  THEN  4050 
4040   Nl =  N1  <  J,  I  ) 
4045   GOTO  4080 

4050   POINTEF  X<I)+.32*Lr,Y<I>-.45*Hr,2 
4055   WAIT  2000 

4060   DISP  "ENTER  THE  NUMBER  OF  ";A1*;"  BETWEEN  REGIONS   " ; CHR*<132) ; I ; CHR$<12: 
;  "  AND  ";CHR*<132>;  I+Nxr;CHR*<123>  j 
4065   INPUT  Nl 
4070   Nl =ABS<N1 ) 
4075   Nl CJ+2, I+Nxr)=Hl ( J, I >=N1 
4080   GOSUB  Leads_label 

4035   MOVE  X<I+Nxr)+.32*Lr, Y< I+Nxr>+. 42*Hr 
4090   GOSUB  Leads_l abel 
4095   GOTO  4  165 
4100   MOVE  X(I),Y<  I  )-. 45-Hr 
4105   IF  OldeOPI  THEN  4  120 
4110   N1=N1(J,I> 
4115   GOTO  4  150 

4120   POINTER  X<I>,Y<I)-.45*Hr,2 
4125   WAIT  2000 

4130   DISP  "ENTER  THE  NUMBER  OF  ";A1*;"  BETWEEN  REGIONS   " ; CHR*< 132) ; I ! CHR*( 12: 
; u  AND  ";CHR*<132>; I+Nxr;CHR*(123); 
4  135   INPUT  Nl 
4140   N1=ABS(N1> 
4145   NKJ,  I)»HKJ+2,  I+Nxr)=Nl 
4  150   GOSUB  Leads  1 abel 
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4155  MOVE  X<I+Nxr),Y<I+Nxr)+.42*Hr 

4160  GOSUB  Leads_label 

4  165  J  =  2   !  HERE  IS  FOR  RIGHT  SIDE  OF  REGION 

4170  IF  I  NOD  Hxr<>8  THEN  4  195 

4175  NK2,I)=N1=0   !  HEPE  IS  THE  RIGHT  SIDE  OF  THE  BOARD 

4180  MOVE  X<I)+.42*Lr,Y<I> 

41 35  GOSUB  Leads  1 abel 

4190  GOTO  4340 

4195  IF  (Itype<I><>9)  OR  C Itype< 1+1 ><>9J  THEN  4270 

4200  MOVE  X<I)+.42*Lr, YCI)-.32*Hr 

4205  IF  01de<  >PI  THEN  4220 

4210  N 1 =  N 1  <  J ,  I > 

4215  GOTO  4250 

4220  POINTER  X<I)  +  .42*Lr,Ya>-.32*Hr,2 

4225  WRIT  2000 

4230  DISP  "ENTER  THE  NUMBER  OF  ";R1*;"  BETWEEN  REGIONS   " ; CHR*< 132) J  I ; CHRf < 123) 

;  -  and  ";CHR*<132>;  l  +  i;CHR*ci28>-; 

4235  INPUT  Nl 

4240  N1=ABS<N1) 

4245  Nl  (  J,  I  >  =  M  1  <  1+2,  I  +  1  :•  =H  1 

4250  GOSUB  Leads_label 

4255  MOVE  X<I  +  l>-.42*Lr,Y<I  +  n-.32*H> 

4260  GOSUB  Leads  1 abel 

4265  GOTO  4335 

4270  MOVE  X<I>+.42*Lr,Y<I) 

4275  IF  01de<  >PI  THEN  4290 

4230  N1=NUJ,I> 

4285  GOTO  4315 

4290  POINTER  X<I>+.42*Lr,Y<I>, 2 

4295  WAIT  2OO0 

4300  DISP  "ENTER  THE  NUMBER  OF  ";A1*;"  BETWEEN  REGIONS   " ',  CHR*<  132) ;  I J  CHR*<  123.1 

;"  AND  ■;CHR*<132J; I+l;CHR*<123>; 

4305  INPUT  Nl 

4310  Nl«flBS<N1) 

4315  GOSUB  Leads_1abe1 

4320  Ml  <  J  +  2,  I  +  1>=NKJ,  I  >=N1 

4325  MOVE  X<I+l)-.42*Lr,Y<I) 

43  30  GOSUB  Leads_l abe 1 

4335  J=3   !  HERE  IS  FOP  TOPS  OF  EACH  REGION 

4340  IF  I>Nxr  THEN  4360 

4345  Nl <J, I)=N1=Q 

4350  MOVE  X<I>,Y<I)+.42*Hr 

4355  GOSUB  Leads  1  abel 

4360  J=4   !  HEPE  IS  FOP  LEFT  SIDE  OF  REGION  OR  BOARD 

4365  IF  <I-1>  MOD  Nxr=0  THEN  4375 

4370  GOTO  4390 

4375  Nl<J,I<=Nl=y 

4330  MOVE  X<I)-.42*Lr,Y(I> 

4335  GOSUB  Leads_label 

4390  NEXT  I 

4395  GRAPHICS 

4400  MOVE  160,130 

4405  LOPG  6 

4410  CSI2E  3 

4415  LABEL  USING  "K";"*  OF  "««fll* 

4420  DRAW  X<Nxr),Y<Hx-r)+Hr/2 

4425  POINTER  O.O.O 

4430  WAIT  5000 
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4435   EXIT  GRAPHICS 

4440   RETURN 

4445  Leadi_chi«:  GOSUB  Error 

4450  Leid:-_chi:  S*="CORRECTIOHS  TO  MUMEER  OF  "::.hU 

4455  ON  ERROR  GOTO  Leadi_chaf 

4460  Check=0 

4465   GOSUB  Pageheid 

44  70   PRINT      You  may  now  make  corrections  to  the  numbers  of  " ;  A 1  X  ;  "  crossing 

the  boundar i es. " 

4475   PRINT  "Recall  that  k3  will  provide  a  hard  copy  of  ths  graphics  if  needed 

uhi 1 e  k2  " 
44S0   PRINT  "will  return  the  graphics  picture  to  the  scrsen  for  visual  checking. 

4485   PRINT 

4490   PRINT  "In  responce  to  ths  prompts  below,  INPUT  Reg  * ,  Reg  *,  CORRECT  »  OF 

"; A1*,LIN< 1 ; 

4495   PRINT  "IF  THERE  ARE  NO  CHANGES  PRESS  CONT  KEY  WITH  NO  INPUT." 

4500   LORG  5 

4505   IF  Check >0  THEN  EEEP 

4510   IF  Check>0  THEN  PRINT  L I N  ■,  3  >  ,  "  CHECK  O'-'ER  THE  FIGURE,  YOU  HAVE  MADE  AN  IL 

LOGICAL  CHOICE  OF  REGIONS." 

4515   Check«Check+l 

4520   Regl=0 

4525   DISP  "REG  *  r  REG  *  j_    CORRECT  *  OF  ";AU; 

4530   INPUT  Regl,Reg2,Newnl 

4535   IF  Regl=0  THEN  RETURN 

4540   IF  <Regl>Nreg)  OR  CReg2>Nreg)  THEN  4465 

4545   IF  Reg2>Regl  THEN  4565 

4550   Dum=Peg2 

4555   Reg2=Regl 

4560   Regl=Dur„ 

4565   J»0 

4570   IF  <Regl+l=Reg2)  AND  'Reg!  MOD  Nxr  >0  •  THEN  J=2 

4575   IF  Regl+Nxr=Reg2  THEN  J=l 

4530   IF  J<  >0  THEN  4655 

4585   EXIT  GRAPHICS 

4590   S**"ERRQR  IN  CORRECTIONS" 

4595   GOSUB  Pagehead 

4600   PRINT  LIN<5>,SPfiC5), "  THOSE  TUG  REGIONS  DO  NOT  CONNECT  TRY  AGAIN" 

4605   EEEP 

4610   WAIT  2500 

4615   GOTO  Leads  cha 

4620   !  HERE  IS  THE  CORRECTION  SCHEME 

4625  Leads  erase:  FEN  -1 

4630  GOSUB  Leads_1abe1 

4635  Nl=Newnl 

4649  PEN  1 

4645    GOSUE    Leads_l abel 

4650  RETURN 

4655  IF  J=2  THEN  4745 

4660  IF  (Itype<Regl)<>10)  OR  CItype<Reg2)  >16)  THEN  4705 

4665  MOVE  X<Regl)  +  .32*Lr,Y<RegO-.45*Hr 

4670  Nl*NKJ,Regn 

4675  GOSUE  Leads  erase 

46S0  MOVE  X<Reg27+.32*Lr,Y<Reg2>+..42*Hr 

4635  N1*NUJ+2,Reg2) 

4690  GOSUB    Leads_erase 

4695  HI ( 1 ,Regl >»N1 < 1 . Reg2)*Newnl 
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4700   GOTO  Leads  change 

4705   MOVE  X<RegT>,Y<Regl>-.45*Hr 

4710   Nl =H1 < J.Regl ) 

4715   GOSUB  Leads  erase 

4720   MOVE  X<Reg2>,Y<Reg2)+.42*Hr 

4725   N1«N1 <J+2,Reg2) 

4730   GOSUB  Leads_erase 

4735   Ml  O,  Reg!  >-NK  J+2,  Reg2>=Newnl 

4749   GOTO  Leads_cha 

4745   !  HERE  WE  APE  TO  CORRECT  THE  RIGHT  SIDE  OF  REGIONS 

4759  IF  <Itype<Regl><>9)  OR  <  Itype<Reg2><>9:>  THEN  4795 
4755   MOVE  X<Regn  +  .42*Lr,Y<Regl>-.32*Hr 

4760  Nl«NKJ,Regl) 
4765   GOSUB  Leads  erase 

4770   MOVE  X<Reg2")-.42*Lr,Y<Reg2>-.32*Hr 

4775   N1*NUJ+2,Reg2) 

4730   GOSUB  Leads  erase 

47S5   NKJ.Regl  >*N1  •  J  +  2,  Reg2)sNewnl 

4790   GOTO  Leads_change 

4795   MOVE  X<Regl>  +  .42*Lr,Y*!Regl) 

43O0   M1=N1 < J,Regl> 

4805   GOSUB  Leads  erase 

4810   MOVE  X<Reg2T-.42*Lr,Y<Reg2) 

4815   Nl=NKJ+2,Reg2 

4320   GOSUB  Leads_eras« 

4825   N1 < J, Regl )=N1 <J+2,Reg2 '=Newn1 

4830   GOTO  Leads^cha 

4335  Tempi n:    !  THIS  SECTION  INPUTS  THE  KNOWN  TEMPERATURES  OR  POWERS  OF  ELEMENTS 

4840  EXIT  GRAPHICS 

4345  S*=" TEMPERATURES  OR  POWER  LEVELS  OF  COMPONENTS" 

4350  GOSUB  Pagehead 

4355  PRINT  "    The  thermal  model  used  by  this  program  assumes  each  component  to 

be  a  heat  " 

4860  PRINT  "source  tor  which  the  user  specifies  either  the  maximum  junction  temp 

erat ure  or-  " 

4365  PRINT  "the  rate  of  heat  generation  within  that  component.   When  the  maximum 

junct  i  on" 
4870  PRINT  "temperature  is  specified,  the  maximum  stead'-'  state  power  levels  are 
calcul  ated.  " 

4975  PRINT  "When  the  rate  of  power  di si  pat  ion  is  specified,  the  steady  state  jun 
ct i on" 

4330  PRINT  "temperatures  are  calculated." 
4835  PRINT 

4890  PRINT  "The  program  uses  the  component  surface  temperature  in  the  thermal  mo 
del  and  " 

4395  PRINT  "therefore  requires  a  case  to  junction  thermal  resistance  <Rj_c>;  how 
ever, " 

49O0  PRINT  "if  rero  is  specified  then  the  surface  'temperature  is  assumed  to  be  t 
he  same  as" 

4905  PRINT  "the  junction  temperature.   When  components  span  two  regions,  enter  h 
alf  the  " 

4910  PRINT  "component  power  for  each  region." 

4915  PRINT  LIN<1>,"     AH  entries  must  be  i  r,  Deg  C  or  Watts  and  deg  C  'Watts." 
4920  IF  OldeOPI  THEN  4945" 

4925  DISP  "PRESS  CONT  WHEN  READY  TO  VIEW  THE  DATA  FROM  "  J  CHR*  <  132  > ',  Name*;  CHRf  < 
123);"  FOR  POSSIBLE  CHANGES" 
4930  PAUSE 
4935  GOTO  Temp  cha 
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4940  ! 

4945  Ans*="l" 

4950  NAT  Pow=ZER 

4955  MAT  Tj=ZER 

4960  NAT  Rj_c=ZER 

4965  INPUT  "DO  YOU  DESIRE  TO  SPECIFY  POWER  LEVELS  <1)  OR  TEMP'S  <2)  ?",Rns* 

4970  IF  Rns*»"BACK  UP"  THEN  RETURN 

4975  Rns*VRL<RnsS) 

4930  IF  <Ans=l)  OP  <flns»2)  THEN  4995 

4935  GOSUB  Err  in 

4990  GOTO  Temp  in 

4995  IF  flns«l  THEN  Pow_in 

5000  Te»p_in:Tem_sol=0   '  HERE  INPUT  TEMPERATURES  WILL  SOLVE  FOP  POWERS  LATER 

5005  3*=" INPUT  OF  JUNCTION  TEMPERATURES" 

501O  GOSUB  Pagehead 

5015  PRINT  "     You  are  nou  entering  junction  temperatures  (deg  C  >  and  junction 

to  case" 

502O  PRINT  "thermal  resistance  <deg  C'Watt).   Do  not  use  kO  <Back_up)  option  u 

h  i  1  e  entering" 

5025  PRINT  "the  dat a. " 

5030  PRINT  LINO) 

5035  PRINT     I  T  j(  I  >  P  j  ,;  <  I  >    " 

5040       FOR  1=1  TO  Nreg 

5045       IF  [type<I)*0  THEN  5863 

5050       DI3P  "FOP  ELEMENT  IN  REGION "; I ; "ENTER   Tjunc , Rj_c " ; 

5055       INPUT  Tj<I),Rj_c  I) 

5068       PRINT  SPR<  1  <  .  I  ,  T  j  I  :i)  ,  F  j_c  <  I  > 

5065       Tj<I)=Tj<I)+273 

5070       NEXT  I 

5075   GOTO  Temp_cha  !  END  OF  TEMP_IN 

5O30  Pow_in:  Tem  sol*l    '  HERE  INPUT  POWER  LEVELS    WILL  SOLVE  FOR  TEMPERATURES 

5035   Sf=" INPUT  OF  COMPONENT  POWER  LEVELS" 

5090   GOSUB  Pagehead 

5095   PRINT  "     You  are  now  entering  component  power  dissipation  (Watt)  and  jun 

c  t  i  o  n  to" 

510O   PRINT  "ca  =  e  thermal  resistances  (deg  C/Watt).  Do  not  uss  k@  <Back_up  whi 

1 e  entering" 

5105   PRINT  "the  dat.  a.  " 

5110   PRINT  LIN'  :  i 

5115   PRINT     I  Pow<I)  Rj  c<I>     " 

5120       FOR  1=1  TO  Nreg 

5125       IF  ItypeU)«8  THEN  5140 

5130       DISP  "FOP  ELEMENT  In  REGION  #" J  I J "ENTER  Pow< I ) , Rj_C < I ) " J 

5135       INPUT  Powa),Rj_C<I) 

5140       PRINT  SPRU),  I,Pow<I),Rj  c<I) 

5145       TjCI>=273 

5150       NEXT  I 

5155   GOTO  Ternp_cha.  !  END  OF  P0W_IN 

5160   ! 

5165  Temp_chae!  GOSUB  Error 

5170  T*mp_cha: !  THIS  SECTION  ALLOWS  CHANGES  TO  ThE  TEMPS  0R  POWER  LEVELS 

5175  Sf="DATA  FOR  "fcPictlbl* 

5180  GOSUB  Pagehead 

5185  PRINT  "     The  data  listed  below  are  the  current  values  for  the  variables 

specified. " 
5190  PRINT  "REG  MM  *     Tiunc  -'BegC-'     Power  (Watts)   Po  ase-.i(W  'C  '  " 
5195       FOP  1=1  TO  Nreg 
520O       PRINT  TRB<2),  I J  TRB<17)  ,  Tj  (  I  >-273; TAB<33) , Pow<  I  )  ;  TAB':  43  <,  R  j_c  '  I) 
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5205       NEXT  I 

5210  PRINT  LIN<2) 

5215  IF  Skip=PI  THEN  GOTO  5350 

5220  PRINT  "Use  the  DISPLAY  up-arrow  OR  down-arrow  to  move  the  data  list  as  requ 

ired" 

5225  Check=0 

5230  flns*="0" 

5235  INPUT  "ANY  CHANCES  ?  INPUT  REG  *  IF  YES  OR  PRESS  CONT  FOR  NO  CHANGES" ,  flriS 

5240  IF  Rns**"BACK  UP"  THEN  Temp  in 

5245  R«gl=INT<VAL<Ans*>> 

5250  IF  CRegl>-l)  AND  CRegKNreg+l  >  THEN  5265 

5255  GOSUB  Erf  in 

5260  GOTO  5230 

5265  IF  <Regl*0)  RND  CCheck=3)  THEN  Hard 

5270  IF  (Regl*0)  RND  (CheckOO)  THEN  Temp_cha 

5275  DISP  "TO  CHANGE  VALUES  IN  REGION  #";Regl; 

52S0  Check=l 

5235   IF  Tem_sol=i  THEN  DISP  "ENTER   Power  ,  Rj_c"| 

5290   IF  Tem~sol=0  THEN  DISP  "ENTER   Tjunc  ,  Pj_c"; 

5295   IF  Tem_sol=l  THEN  INPUT  Pow<Reg  1  > ,  R.j_c  (Regl  > 

5300   IF  Tem_sol=0  THEN  INPUT  T  j  v  Reg!  >  ,  R.j_c  <Reg  1 ) 

5305   IF  Tem~sol=0  THEN  Tj  <Regl  )=T.j  CRegl  >+273 

5310   PRINT  LIN<l),TRB<2>,Regl;TAB<17),Tj CRegl >-2?3; TAB < 33) , PowCRegl > ; TAB < 48) ,Rj 

_c<Regl>;"    ***  CHANGE  ***" 

5315       GOTO  52  30  !  END  OF  TEMP  CHA 

5320  Hard: '  THIS  SECTION  PRINTS  A  COPY  OF  THE  INPUT  DATA  IF  REQUESTED 

5325  Ans*="N" 

5330  INPUT  "DO  YOU  DESIRE  A  PRINTED  COPY  OF  THE  DATA  ABOVE  (N  or  Y>",Ans* 

5335  IF  Ans*="BflCK_UP"  THEN  Temp_cha 

5340  IF  flns**"N"  THEN  RETURN 

5345  PRINTER  IS  0 

5350  Skip=PI 

5355  GOTO  5190 

5360  PRINT  TRB<39-LEN<"THE  ABOVE  DATA  IS  FOP  "&Pic  t  1  b  1  S  >  ■  2  )  ,  "THE  ABOVE  DATA  IS  F 

OR  " ;Pict lbl s . LIN ■  1  ' 

5365  PRINT  "BOARD  LENGTH  (defined  along  air  f 1 ou)*" ; Bdl J "mm  HEIGHT  ="; 

Bdh; "mm" ,  LINO) 

5370  PRINT  "BOARD  THICKNESS*" ; Thi ck_b; "mm  CONDUCTIVITY  =" 

;Kb; "Watts'M-K" , LIN<2> 

5375  PRINT  "THE  BOARD  MODEL  ASSUMES  ";A1*;"  AS  CONDUCTION  F-THS  WITH  AN  AREA  OF 

";AB3(A1  >; "  mm- 2".LIN'  1  i 

53S0  PRINT  "THERMAL  CONDUCTIVITY  OF  THE  ";A1*;"  =  "  ;  K  1  ;  "  Watts/-M-C" 

5385  PRINT  " 

",LIN<2) 


5390  PRItlTER  IS  16 

5395  RETURN  .  !  END  OF  Hard 

5400  ! 

5405  St  owe:  GOSUB  Error 

5410  Stow:    !  THIS  SECTION  PLACES  THE  DESCRIPTION  OF  THE  CIRCUIT  BOARD  ON  TAPE 

5415  ON  ERPOR  GOTO  St  owe 

5420  S*="PECGPD  BOARD  DESCRIPTION  ON  MASS  STORAGE" 

5425  GOSUB  Pagehead 

5430  PRINT  "     You  may  record  all  the  data  concerning  the  circuit  board  on  any 

avai  1 abl e" 

5435  PRINT  "mass  storage  device.   This  allows  any  user  to  retrieve  the  descripti 

on  at  some" 

5449  PRINT  "later  time  without,  the  need  to  input  all  the  details.   This  option  o 

ccurs  both" 
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5445    PRINT    "before    and    after    the    thermal     analysis.       Enter    desired    data   file    rn 

me    bel ow, " 

5450  PRINT  "be  sure  to  include  the  mass  storage  unit  specifier  if  not  the  defaul 

t  .  " 

5455    PRINT    "For    example  :T14     ,     :F8    ,     Yt2    etc." 

54(50    Rns*='*Y" 

5465  INPUT  "HO  YOU  DESIRE  TO  RECORD  THE   DESCRIPTION  DATA  CV  or  N)?",Ans* 

5470  IF  Ans**"BACK_UP"  THE'I  RETURN 

54:30  IF  UPC*<Ans*Cl,13)="N"  THEN  RETURN 

54S5  IF  OldeOPI  THEN  5505 

5490  PRINT  LIN<5),"THE  PRESENT  DESCRIPTIVE  TITLE  FOR  THAT  BOARD  IS  ";CHR*U32);P 

i  ct Ibl *;  CHR*',  12S) 

54?5  EDIT  "CHANGE  THE  TITLE  OR  FRESS  CONT  WITH  NO  ENTRY  FOR  NO  CHANGE.  ",  Pi  cr.  1  b 

1$ 

5500  IF  Pictlbl**"BRCK_UP"  THEN  Stow 

5505  Rns*=Name* 

5510  EDIT  "UNDER  WHAT  NAME  DO  YOU  DESIRE  TO  STORE  THE  DATA  (change  bel ow)?" , Name 

$ 

5515  IF  Na»e*=MBACK_UP"  THEN  5490 

5520  IF  Name*=Ans*  THEN  5560 

5525  Ans*="Y" 

5530  DISP  "DOES  A  DATA  FILE  WITH  AT  LEAST" ; 4S*Nreg+S00; "  BYTES  EXIST  UNDER  THAT 

NAME  <Y  or  N > ?  "  ; 

5535  INPUT  fins* 

5540  IF  Bns**"BACK_UP"  THEN  5505 

5545  IF  UPC*<flns*Cl,  U)  =  "Y"  THEN  5560 

5550  DISP  "WORKING  CREATING  DATA  FILE  FOR  ";Name*;"  THAT  IS" ; 40*Nreg+600;  "BYTE 

S  IN  SIZE" 

5555  CREATE  Name*, 1 , 40*Nreg+S00 

5560  ASSIGN  #5  TO  Name* 

5565  PRINT  LINC3), "WORKING  WRITING  BOARD  DESCRIPTION  OF  " ; Pi ct 1 bl *J "  ON  MASS  S 

TOPAGE" 

5570  DISP 

5575  READ  #5, 1 

5530  PRINT  »5;  Name*,  Pi  ctlbl*,  Case*,  Bdl  ,Bdh,Sc  a,Nxr,Nyr,  Nreg,  Thi  ck_b,  Kb,  A1  ,  Kl  ,  Ter„ 

_sol 

5535  FOR  1=1  TO  Nreg 

5590  PRINT  #5;  Itype<I),X<I),Y<D,Tj<I  >,Pow<  I  > ,  R  j_c  <  I ) 

5595  FOR  J=l  TO  4 

5600  PRINT  #5;N1 < J, I) 

5605  NEXT  J 

5610  NEXT  I 

5615  PRINT  *5;END 

5620  ASSIGN  *Z    TO  * 

5625  PRINTER  IS  8 

5630  PRINT  LIN<2),UTHE  CIRCUIT  LISTED  BELOW  IS  STORED  UNDER  THE  FILE  NAME  "JCHR* 

(132); Name*; CHR*< 128) 

5635  PRINT  LIN<2> ,SPA<25) , Pi ct 1 bl *, LIM<2) 

5640  PRINT  SPfi<25),"SAVE  FOR  YOUR  RECORDS" 

5645  PRINT  " 

__"  ,  L  I  N  ■  2  :• 


5650  PRINTER  IS  16 

5655  DISP 

5660  RETURN 

5665  Error:! 

5670  BEEP 

5675  WAIT  300 

5680  IF  ERRN=56  THEN  Err  name 
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5635  PRINT  LIN<20>,SPfl<lO>, "ERROR  NUMBER" ; ERRN; "HAS  OCUPPED  IN  LINE" ; ERRL; " .   PR 

ESS  CONT  WHEN  READY" 

5690  DISP 

5695  BEEP 

5700  PAUSE 

5705  RETURN 

5710  ' 

5715  Err_nime: !  THIS  SECTION  FOR  IMPROPER  FILE  NAME 

5720  PRINTER  IS  16 

5725  PRINT  PAGE 

5730  Msusfg'TErrilJLT  MASS  STORAGE" 

5735  FOR  1=2  TO  LEM<Naroef) 

5740  IF  Name* CI ,  I  ]  =  ": "  THEN  5755 

5745  NEXT  I 

5750  GOTO  5770 

5755  Msus*«N*me*CI] 

5760  CAT  rtsus* 

5765  GOTO  5775 

5770  CAT 

5775    PRINT    LIN<2>,CHR*<I32>;Name*Cl, I-1];CHR*<123>; "     is    NOT    on    "JMsus*;"    with    th 

at    spel  1  i  rig " 

5730  PRINT  LINO), "CHECK  OVER  THE  DIRECTORY  ABOVE  FOR  CORRECT  NAME 

5785  DISP  "PRESS  CONT  WHEN  READY" 
5790  PAUSE 
5795  RETURN 
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THERMAL  ANALYSIS  PROGRAM  FOR 


•  ****- 


>♦****#*#********♦*        thermele: 

a  system  of  programs  for  the  hp  9845 
to  perform  thermal  analysis  of  electronic  circuit  boards 


PREPARED  AT  THE  NAVAL  POSTGRADUATE  SCHOOL  MONTEREY,  CA 
R.  A.  FOLTZ  LCDR  USN 
Mill   THERML   Mill  ........ TiiME  1 


:k  b, Tern  so  1 , Kb, Al , Kl 


,R1<4, 5@) , R#l 


Tb' 30: 


so 

90 

100 

1  10 
120 
130 
140 
150 
160 
170 
ISO 
190 

200   OPT  I  OH  BASE  1 
210   PRINTER  IS  1-; 

2  20   DIM  S*C?0],Type*'  0:  14  •  , Msus  1 1221 ,  Y 1 b 1  I C 23 ] 
230   COM  Map. Name*, Pict I bl *C50] , Case* C 50] , Bdl . Bdl 
240   COM  SHORT  X<50),Y<50),Tj<:50),Pow<50),Rj  c<50) 
250   COM  INTEGER  Itype<50) ,  NK4,  50>  ,  Nxr,  Nyr,  Nreg 
260   INTEGER  Npin< 15) 

2?0   SHORT  fie  C50  ) ,  ?e_:  ond  <  30 ) ,  Re_c  onw  <  50  ) ,  Rtot  e_b  <  50  - ,  Pr  op. 
0) , Rb_conv<50) 

2S0   SHORT  Wideset  14) , Lenset < 14) 

290   SHORT  Le<50),We<50>  He<50> , Te<50) , Tai r<50> ,  FK50, 50> , B< 
300   Map* INT < Map) 

310   IF  <Map=0)  OR  (Map>3)  THEN  Oldpict 

320  Start_over:   ON  Map  GOTO  01  dpi ct , Thermal , Thermal 

330   GOTO  Oldpict 

340  Oldpicte!  GOSUB  Error 

350  Oldpict:  !  THIS  SE:t'::.  RETRIEVES  n  BOARD  DESCRIPTION  OFF  A  MASS  STORAGE  DEV 

ICE  FOR  THE  PURPOSE  OF  DEBUGGING  THERMAL 

360   ON  ERROR  GOTO  Oldpict* 

370   01de*PI 

380   S*-'"BQfiRD  DESCRIPTION  FROM  MASS  STORAGE  DEVICE" 

390   GOSUB  Pagehead 

400   PRINT  "     You  have  chosen  to  input  the  circuit  board  description  in  THERM 

L  di  rect 1 y" 

4  10   PRINT  "from  a  mass  storage  device.   This  program  in  the  THERMELEX  System  wi 

1 1  NOT  " 

420   PRINT  "allow  graphical  data  checking  and  while  taster,  there  is  the  c nances 

that  the" 
430   PRINT  "data  is  incorrect.   It  you  decide  that  it  would  be  better  to  check  t 
he  data" 

440   PRINT  "press  K0  <Back_up)  and  BOARDS  will  be  loaded  from  the  DEFAULT  ma 
ss  storage  un it." 
450   PRINT 

460   PRINT  "     The  data  Pile  containing  the  beard  description  must  have  been  st 
ored  by  " 

470   PRINT  "this  program.   Enter  the  data  file  name  below,  be  sure  to  include  th 
e  mais  " 


430   PRINT  "storage 
490   PRINT  LIN<1), " 


f  needed. 


T14, : F3, 


Yl; 


Do  NOT  use  quotes' 
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5a0  EDIT  "UNDER  WHAT  FILE  NAME  IS  THE  CIRCUIT  BOARD  FILED  (change  or  enter  belo 

w)  " , Name* 

510  IF  Name*»"BACK  UP"  THEN  53d 

520  GOTO  550 

530  Map=l 

540  LOAD  "BOARDS", 1 

550  ASSIGN  #1  TO  Name* 

568  READ  #1,1 

570  READ  #1 ;  Name*,  Pict  lblt-.Casef,  Bdl  ,  Bdh,  Sca,Nxr  ,Nyr,  Nreg,  Thi  ck_b,  Kb,  fll  .  Kl  ,  Tem_ 

sol 

580  FOP  1=1  TO  Nreg 

590  READ  #1  ;  [type  •  I  >,X(I),YCI  > ,  T  j  <  I>,Pew<I  >,  R  j_c  <  I) 

600  FOR  J=l  TO  4 

610  READ  ilJNl <J, I) 

620  NEXT  J 

630  NEXT  I  !  END  OF  OLDF'ICT 

640  ! 

650  GOTO  Thermal 

660  Therm  ale:  GOSUB  Error 

670  Thermil :    !  THIS  SECTION  IS  THE  MAIN  INPUT  AND  CALLING  ROUTINE 

6  S  0  0  N  ER  P  0  F  G  0  T  0  T  h  e  r  m  a  1  e 

650  £*=" THERMAL  ANALYSIS  OF" 

700  GOSUB  Pagehead 

710  PRINT  TAB'.  3?-LEN<Pict  lbl*  i  2  > ;  CHR*<  133)  ;  Pi  ct  1  bl  *;  CHR*<  123);LIN<2) 

720  PRINT  "    This  section  assumes  you  hive  completely  and  correctly  described 

the  board" 

730  PRINT  "itself.   You  will  be  asked  questions  concerning  ONLY  -.he-  environment 

740   PRINT  "     The  first  questions  are  concerned  with  the  cooling  air  supply. 

Recall  the" 

750   PRINT  "direction  of  air  flow  on  the  graphics  picture  is  assumed  to  be  left 

to  right . " 

?€■*      PRINT  "The  clearance  between  the  boards  is  used  to  determine  the  velocity  o 

f  the  " 

770       PRINT     "cooling     ai r . " , L I N < 2 ) 

760      flns$*u20" 

790      Map=6 

800   INPUT  "ENTER  THE  INLET  TEMPERATURE  OF  THE  COOLING  A  IP  •  'deg  CV'.Ans* 

810   IF  Rns$»"BACK_UP"  THEN  Start_over 

820   Tair*VALCAns#) 

830   PRINT  "     INLET  AIR  TEMP  T a i r= " ; T ai r ; " de g  C" 

340   Tai r =Tai r+2?3    !  ALL  CALCULATIONS  DONE  IN  ABSOLUTE  TE.  P 

850   Ansi=".0005" 

360   INPUT  "ENTER  THE  AIR  SUPPLY  PEP  BOARD  <MA3^SEC) " , Ans* 

870   IF  Ans*="8ACKJJP"  THEN  Thermal 

880   Fai  r=VAL'  Rns*) 

890   PRINT  "     FLOW  PATE  OF  A  I P= " ; F ai r j  " M  3-  Sec " 

900   flnst»"15.24" 

910   INPUT  "ENTER  THE  DISTANCE  FROM  THE  FACE  OF  THE  BOARD  TO  THE  NEXT  OBJECT  (mm 

)"  ,fin±$ 

920   IF  Ans*="BACK_UP"  THEN  7S0 

930   Zb*VAL<Ans*> 

940   PRINT  "     BOARD  SPACING=";Zb; "mm" 

950   Zb=Zb*.0Ol 

960   Ans*="NULL" 

970   INPUT  "ALL  OK?.. PRESS  CONT  ANY  ENTRY  FOLLOWED  BY  CONT  HILL  ALLOW  REENTP 

Y  OF  ALL",Ar,if 

980   IF  fins*»"NULL"  THEN  1000 


US 
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990   GOTO  Thermal 

1000  S*»"CQNVERGENCE  CRITERIA" 

1910  GOSUB  Pagehead 

1020  PRINT  "    Conyerg?nce  is  indicated  by  successive  itterations  that  result  in 

element" 
1030  PRINT  "temperatures  that  differ  only  by  some  small  amount.   Each  element  te 
mperature" 

1040  PRINT  "is  compared  to  that,  obtained  in  'he  previous  ltteration  and  if  the  1 
ar  g  e  s  t  " 

1050  PRINT  "difference  is  less  than  a  maximum  specified  error,  results  are  print 
ed.  " 

1060  PRINT 
1070  PRINT  "    Typically  tuo  or  three  i iterations  result  in  a  maximum  difference 

on  the" 
1030  PRINT  "order  of  one  degree  centigrade  when  solving  for  temperatures  and  fou 
r  or  five" 

1030  PRINT  "itterations  will  result  in  a  maximum  difference  on  the  ord=r  of  .1  W 
att  uhen  " 

1100  PRINT  "solving  for  powers." 
1110  IP  Tem_sol=l  THEN  Ans*=".5" 
1120  IF  Tem_sol=0  THEN  Mn;{="l" 

1130  IF  Tem_sol=l  THEN  INPUT  "ENTER  THE  MAXIMUM  DIFFERENCE  BETWEEN  ITTERATIONS  (. 
deg  C)  DEFAULT*. 5", fins* 

1140  IF  Tem_sol=0  THEN  INPUT  "ENTER  THE  MAXIMUM  PERCENT  CHANGE  BETWEEN  ITTERATIO 
NS  (Watt)  DEFAULT- lJt",Ans* 
1150  IF  Ans*»"BACK_UP"  THEN  Thermal 
1 160  Errmax»ABS<VAL< AnsS) ) 

1170  GOSUB  Database  !  THE  FOLLOWING  LINES  FORM  THE  MAIN  CALLING  ROUTINE 
1180  GOSUB  CalcO 
1190  GOSUB  Calcl 
1200  GOSUB  Calc  air 
1210  GOSUB  Calc_t 
1220  GOSUB  Calc2 
1230  GOSUB  Debug 
1240  GOSUB  Solve 
1250  GOSUB  Units 
1260  IF  Bomb=l  THEN  1280 
1270  GOSUB  Output 
1280  GOSUB  What  now 
1290  IF  Map-8  THEN  GOTO  Thermal 
1300  GOTO  Oldpict 
1310  ! 

1320  !  END  OF  THE  MAIN  CONTROL  SECTION  OF  THE  PROGRAM  THERML 
1330  ! 

1340  !*#*******DATABASE******** 

1350  Database:  !   IN  THIS  SECTION  MANY  OF  THE  CONSTANTS  USED  IN  THE  CALCULATIONS 
1360  !   ARE  READ  IN  FROM  THE  DATA  LINES  BELOW 

1370  DISP  "WORKING  ON  NON-CHANGING  PARAMETERS" 
1330  PRINT  PAGE 

1390  Axpin=4.3E-7  JX-SECTIONAL  AREA  FOR  PINS  <MA2> 

1400  Aspin=lE-5  ! SURFACE  AREA  FOR  PINS     <M~2> 

1410  Beta»3.33E-3  !V0L  COEFF  OF  EXPA.  AIR    '1-deq  K>   AT  300  deg  K 

1420  Cpa*l.006E3  'SPECIFIC  HEAT  OF  AIR      <  UATT-i EC-'Kg-deg  K> 

1430  Dis=.001  'AVG  DISTANCE  FROM  BOTTOM  OF  DIF  TO  BOARD 

1440  Epsb=.3  IEMMISIVITY  OF  THE  SURFACE  OF  BOARD 

1450  Epse=.9  IEMMISIVITY  OF  THE  DIP  SURFACE 

1460  G-9.31  ! GRAVITY  <M/Sec-2) 

1470  Gnu=1.634E-5  'KINEMATIC  VISCOSITY  AIR   (Mt^/Sec) 
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1480  Hec=.003  ! HEIGHT  OF  CERAMIC  DIP  PACKAGE  -****** 

1490  Hep».eeS  (HEIGHT  OF  PLASTIC  DIP  PACKAGE  ******* 

1300  Ks  =  59  !  THERMAL  COHD  OF  STEEL     UJATTS-- M-d*g  K> 

1510  Kair*.826  ! THERMAL  COND  OF  AIR      <WfiTTS/M-d*g  K)   300  deg  K 

1320  Kpin»334  (THERMAL  COND  OF  PINS      < WATTS'M-deg  K) 

1530  Lpin=,0025  '  AVG  LENGHT  OF  PINS        CM) 

1540  Mu=1.983E-5  (DYNAMIC  VISCOSITY  300  K   (Kg/M-Sec) 

1550  Pr=.703  IPRRNDTL  NUMBER  <AT  SOOdeg  K^ 

1360  Rho= 1.1774  (DENSITY  OF  AIR  AT  300  K   <.Kg.-M-3> 

1570  Sig=5.67E-3  ! STEFFAH-BOLTZMAN  iW/M~2-dsg  KA4> 

1530  ltt=0 

1590  RESTORE  1630 

1600  FOR  1=1  TO  14  i  THIS  READS  THE  CASE  WIDTHS  FOP  EACH  TYPE 

1610  READ  Uides«t<I) 

1620  NEXT  I 

1630  DATA  .25, .725, .25, .325, .55, 1.25,0,0, .6, 1, . 18, .35, .43, .72     (  ALL  IN  INCHES 

1640  FOP  1=1  TO  14  i  THIS  READS  THE  CASE  LENGTHS  FOR  EACH  TYPE 

1650  READ  L»ns*t<I) 

1660  NEXT  I 

1670  DATA  .725, .25,-825, .25, 1,-55,0,0,1, .6, . 13, .35, .43, .72        !  ALL  IN  INCHES 

1630  MAT  Wid«set=<25.4)*Wides*t  !  CONVERT  TO  mm 

1690  MAT  L*nset  =  <25.4)*l_enset 

1700  FOR  1=1  TO  14  (THIS  READS  THE  *  OF  PINS  FOR  EACH  TYPE 

1710  READ  Npir.(I> 

1720  NEXT  I 

1730  DATA  14,14,16,16.24.24.0,0.20,20,16,24.40,64 

1740  FOP  1=0  TO  14  !  THIS  i?EADS  THE  CASE  LABELS  FOP  EACH  TYPE 

1750  READ  Type*<I) 

1760  NEXT  I 

1770  DATA  EMPTY, DIP  14. DIP  14, DIP  16, DIP  16, DIP  24. DIP  24 , NULL , NULL , D I P  40,DIP  4 

0,FLAT16,FLAT24.FLAT4y,FLAT64 

1730  FOR  1=1  TO  Nreg 

1790  IF  It,..p*'I>=0  THEN  1330 

1800  L«<I>=Lenset<  Ityp*<I>) 

1810  We(  I  >*Wideset  •'  Ityp*(  I  > ) 

1820  GOTO  1340 

1830  W«<I>»Le<I>-0 

1840  NEXT  I 

1850  IF  Terr._so1  =  l  THEN  MAT  T»*<300)  !  INITIAL  GUESS  FOR  TEMP  CASE  =  27  deg  C 

1860  IF  Te»"*o1=0  THEN  MAT  ?ow=  .23.'     !  INITIAL  GUESS  FOR  POWER  =  .25  Hi!U 

1870  MAT  Le=( . 0O1  '*Le 

1830  MAT  We=C.001)*We 

1890  MAT  Xa<.091)*X 

1900  MAT  Y=<.001 >*Y 

1910  MAT  A*=L*.We 

1920  MAT  He«<Hec)      !  ASSUME  ALL  CERAMIC  MAKE  CORRECTIONS  AS  NEEDED  IN  CilcO 

1930  Bdl=Bdl*,001 

1940  Bdh=Bdh*,001 

1950   Areg=Bdl *BdlvNreg 

1963  Hr=Boh/Nyr 

1970  Lr=Bdl/Nxr 

1930  Thick  b»Thick  b*. 001 

1990  Al=A1*lE-6 

2000  RETURN  !  END  OF  DATABASE 

2010  ! 

2020  Cal0e:GO3UB  Error 

2030  CalcO:'  ON  ERROR  GO  TO  CALCOE 

2040    Poutot=Ni"g  =  H*ii,.g  =  U.r  avg  =  0 
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2858 
2868 
2878 
2838 
2898 
2188 
2118 
2128 
2138 
2140 
2158 
2168 
2178 
2138 
2198 
2288 
2219 
2228 
2238 
2248 
2258 
2268 
2278 
2288 
2298 
2398 
2318 
2328 
2338 
2348 
2358 
2369 
i  Kfl 
2378 
2338 
2398 
2498 
2418 
2429 
2438 
2449 
2459 
2468 
2479 
2439 
2498 
2588 
2510 
2528 
2538 
2549 
2558 
2569 
2579 
2539 
2598 
2609 
2619 


FOR  1=1  TO  Nreg 

IF  Cas*tCI,I]="l"  THEN  He<I)»Hep 

IF  Itype<I)=0  THEN  2120 

Poutot»Powtot +Pow< I) 

H«av»g=Meavg+<W*<I>  +  .001>-'Nreg   !  RVG  WIDTH*  INCLUDING  PINS 

HeawgssHtfayg+CHeCD+LpinJ^Nreg   !  flVG  HEIGHT  INCLUDING  PINS 

Navg*Navg+l«'Nyr  !  NUMBER  OF  ELEMENTS  IN  AN  RVG  CROSS  SECTION 

NEXT  I 

Rair=Bdh*Zb-Navg*W« &vg*Heavg 

Pen  «»2*Bdh+2*Navg*Heavg+2*Zb 

Fr=l+5*<P*rim-2*Bdh-2*Zb)<Perim 

Dh=4*Aai  r/Per i m 

Vair  =  Fair/fliir 

Re=Vai r*Dn- Gnu 

IF    ReHOOu    THEN    Hbf ».  023*Kai  r/'Dh»R«A.  8 


ROUGHNESS  FACTOR 


IF  Re<  =  1808  THEN  Hbf  *5. 40*Kai  r*Fr/'Dh 

Rl_hor«Lr/'CKl*flBS<fl1  O 
Rl_ver=Hr  'CK1 *H8S<fll  > ) 
Rb_hor=Lr '<Kb*Hr*Thi ck_b) 
Rb~v«r»Hr/'<Kb*Lr*Thi  ck_b> 

RETURN 


TURBULENT  AT  1000  DUE 
TO  MANY  TRIPS  OF  C0MP3 
LAMINAR  FOR  1  Gz>.85 


!  END  OF  CALCi 


Calcl*:GOSUB  Error 

Calcl:  !  THIS  SECTION  COMPUTES  SOME  OF  THE  NON-CHANGING  PARAMETERS 
ON  ERROR  GOTO  Calcl* 
FOR  1=1  TO  Mreg 

1  BELOW  HERE  WE  CALCULATE  THE  CONVECT IVE  LOSSES  FOP  EACH  BOARD  REGION 
IF  XCD'DhMO  THEN  2390  !  OUTSIDE  OF  THE  DEVELOPMENT  REGION 

Gz»R«*Pr*Dh 

Hb«.664*Kair  C 1 .  l*Dh>*SQR«:Gz*<  1+7.  3*SQR<Pr^Gz>  >/"Pr>*Fr   !  Eq  13.43  KNUTZEN 
TZ  LAMINAR  IN  DEVELOPMENT  REGION 
IF  UGz>.05  THEN  2390 
GOTO  2480 
Hb=Hbf 

flr«gtot«2*Rreg-fle< I  • 
Rb  corw<  I  "'  =  l  CHb*Aregt  ot  ) 

!  NOW  WE  GET  TO  THE  ELEMENTS  ON  THE  BOARD 

IF  Itype<I>=0  THEN  2640         '  IF  NO  ELEMENT  THEN  SET  VERY  HIGH  RESISTANCE 
IF  Itype<IXU  THEN  2430 
D1$*. l*Bi2 
Lpi  n=. l*Lpi n 
Aspi n=.  1 *flsp i  n 

Rpi  r»s_cond=Lpi  n/CKpi  n*flxpi  rn*Npi  n<  I  type1;  I  >  >  > 
Rgap_c o  n  d  =  D i  s^  <  Kai  r *fl*  < I ) ) 
IF  [?yp«<  IX 11  THEN  2540 
Di s=10*Di i 
Lpi n=10*Lpi  n 
Aspi  n=lo*Aipi  n 

Re_cond<  I  )=Rpi  n»__cond*Rgap_cond-/  <Rp  i  ns_cond+Rgap_cond) 
Abl  a=-  =  Ae<  I  'j+2*L*':  I  >*He+fispi  n*Npin<  I  type'!  I )  > 
Hebl as=Hb 

IF  I<6  THEN  Hb<  I  '=Hb 
Rebl as»l/CHeb1as*Rb1 ii  • 

H«st ag«. 57*Kai r*Pr". 4*SQR<Vai rv<.We< I )*Gnu) > 
flstag=2*We< I>*He< I ) 
Restag=l/<Hest ag*flst ag) 
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2620    Re_conw( I  >=Peb1 ii*Rest ag^ <Peb1 as+Restag> 

2630    GOTO    2668 

2640    Re_cony( I )=1E30 

2650  Re_conoM  >  =  1E30 

2660  !  BELOW  HERE  CALCULATE  THE  BOARD  CONDUCTIVE  RESISTANCE 

2670  J=l   '  HERE  IS  BOTTOM  OF  REGION 

2680  IF  I*Nxr>*Nreg  THEN  2720 

2690  IF  NKJ,I>»8  THEN  2748 

2700  Rl  <  J,  I)=Rl_u*r--N1  (  J,  I>*Rb_ver .■•  <R1  yer/NKJ,I)+Rb  yer) 

2710  GOTO  2750 

2720  RKJ,I>-Rb  ver'2E-40 

2730  GOTO  2750 

2740  Rl (J, l)»Rb_ver 

2750  J=2   !  HERE  IS  RIGHT  SIDE  OF  REGION 

2760  IF  I  MOD  N*r  =  0  THEN  2800 

2770  IF  NKJ,I>»0  THEN  2320 

2730  Rl  <J,  I)»RI_hor/NHJ,I  >#Rb_hor-'<Rl_hor^Nl  CJ,  O+Rb  hor) 

2790  GOTO  2330 

2303  Rl  <  J,  n=Pb_hor -2E-40 

2310  GOTO  2830 

2820  Rl (J, I>=Rb_hor 

2830  J  =  3    '  HERE  13  THE  TOP  OF  THE  REGION 

2840  IF  I :=Nxr  THEN  2380 

2350  IF  MKJ,I>»0  THEN  2900 

2860  RKJ,  I>»R1  ver>  NKJ,  I  >*Rb_ver/  (Rl  wer/Nl  CJ,  I)+Rb_uer) 

2870  GOTO  2910 

2830  RKJ,  I)«Rb  v«r/2E-49 

2890  GOTO  2910 

2900  Rl CJ, [)>Rb  vsr 

2910  J=4   '  HERE  IS  THE  LEFT  SIDE  OF  THE  REGION 

2920  IF  (1-1)  MOD  Nxr»8  THEN  2960 

2930  IF  M1<J,I)»8  THEN  2930 

2940  Rl < J, I )»R1  hop- Nl (J, I )*Rb  hor/(RI  hor/Nl CJ, I)+Rb  hor ) 

2950  GOTO  3008 

2960  Rl  (J,  U=Pb_hor-  2E-48 

2970  GOTO  30OO 

2930  Rl (J, I>»Rb  hor 

2990  GOTO  3010 

3000  Rb_eonw< I >=l/<Hb*Ar#gtot > 

!  END  OF  CflLi 

GOSUB  Error 

3060  Calc_t.;i  THIS  SECTION  DETERMINES  SURFACE  TEMP  FROM  JUNCTION  TEMP  RND  Rj_C 
3070  ON  ERROR  GOTO  Calc_t* 
3080  IF  Tem  SOl*l  THEN  RETURN 
3090      FOR  1=1  TO  Nreg 
3100      T><I>=Tj<I)-Pou<  I)*Rj_c<I) 

3110      IF  T*CIXTair<I>*l.  1  THEN  Te<  I>*Tair<I  >*l.  1 
3128      NEXT  I 

3130  RETURN  'END  OF  Cal C 

3140  Cal. : 2* : GOSUB  Error 

3150  Calc2:   !  THIS  SECTION  CONTAINS  THOSE  PARAMETERS  WHICH  CHANGE  WITH  TEMPS 
3160  ON  ERROR  GOTO  Calc2« 

3170  DISR  "WORKING  ON  CHANGING  PARAMETERS" 
3180  FOR  1=1  TO  Nrtg 
3190  IF  ItypcCDOG  THEN  32  30 
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NEXT  I 

3020 

! 

3038 

RETURN 

304O 

i 

3050 

Calc  ti 
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3200  Rtop_r<n  =  lE50 

3210  Rgap~rad*lE58 

3220  GOTO  3250 

3230  Rgap_r.=id=(Epi.?  +  Epsb-Ep*e*Ep*bV.'<4*$i  g*Epse*Epsb*fl«  <  I  >*Te<  I  )A3) 

3240  Rtop  r<  I  ■  =  •  ( l-Epsb)'  •:  Epsb-^r*g>-t-10/i=lr*g+<.  1-Epse  )  '  h*  •.  I  >*Epse  ' > '  <4*Si  g*Te<  I  > 

-3) 

3250    Rtote_b<  I  >=Rgap_r ad*Re_cond<  I  )>"<Rgap_rad+Re_corid<  I )  > 

3260    Rele<7)*Re_conv7l>+Rto7e    b<I)  !  !  !  !    MUST    BE    SUM    FOR    RLGEERflIC    REP.SONS 

3270    NEXT    I 

3280  RETURN  !  END  OF  CALC2 

3290  ! 

3300  Calc_ai re: GOSUB  Error 

3310  Calc~air:>  THIS  SECTION  FINDS  THE  AIR  TEMPS  FOR  EACH  REGION  BASED  ON  POWER 

3320  ON  ERROR  GOTO  Calc_aire 

3330  Cfr*Cpa*Fair*Rho 

3340  Powtot=0 

3350     FOR  1=1  TO  Nreg 

3360     Pout  ot  =Powt  ot  +Pow<  I  -1 

3370     Tai r< I >*Tai r+. 5*Pow< I >*Nyr^Cfr 

3330     IF  (1-1)  MOD  Nxr=0  THEN  3420 

3390         FOR  K  =  I - 1  TO  I-CI-1)  MOD  Nxr  STEP  -1 

3400        Tair<I)=»Tair<I  >+Pow<K>*Nyr/'Cf  r 

3410         NEXT  K 

3420     NEXT  I 

3430  Tout^Tair+Powtot'Cfr 

3440  RETURN  !  END  OF  Cl=iLC_flIR 

3450  ' 

3460  Solv»e:GOSUB  Error 

3470  Solve:  !  THIS  SECTION  SOLVES  THE  PROBLEM 

3480  ON  ERPOP  GOTO  Solves 

3490  DISP  "WORKING  ON  SETTING  UP  THE  MATRIX" 

35O0  Itt=Itt+l 

3510  BEEP 

3520  WRIT  300 

3530  BEEP 

3540  GOSUB  Set  up 

3550  GOSUB  Elu 

3563  Err=0 

3570  Tmax-0 

3530  P»t-n*l000 

3590  IF  Tern  sol* 1  THEN  3730 

3600     FOR  1=1  TO  Nreg 

3610      Pr,ew=(T?',  I  >-B<  I  >  >'Rtote_b<  I  >+<Te<  I  >-Tai  r<  I  >  )/Re_conv(  I  >  +  <Te<  I  )-B<I  >  >/Rt 

op_r< I ) 

3620      IF  Pnew<9  THEN  Pnew=Pow< I )/2 

3630      IF  Pneu<Pmin  THE'.'  Pmin=Pnew 

3640      IF  flBS<<Pow<I)-Pnew)/Pow<l)>>Err  THEN  Err=flBS<  <Pow<  I  )-Pnew)^Pow<  I  )  ) 

3650     Pou<I)*Pnew 

3660      NEXT  I 

3670   Taax»500 

3630   MAT  Tb=B 

3690   IF  Err<Errmax/100  THEN  3340 

370O   GOSUB  Calc_t 

3710   GOSUB  Calc  air 

3720   GOTO  3920 

3730       FOR  1=1  TO  Nreg 

3740      Tnew=CPow<  I  >*Rtote_b<  I  )*Re_conv<  I>+Ba>*Re_conv<  I)+Tai  r<  I  >*Rtote_b<  I > ) 

/Rele<I> 
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3750       IF  ABS<T*<I  >-Tnew>  >Err  THEN  Err=AB3  <  Te  <  I )  -Tnew  ) 

3760       T*<I>=Tnew 

3770       Tj<I)»T*<I)+Pow<I)*Rj  c<I) 

3780      IF  Tj<I)>T»ax  THEN  Tmax=T j < I) 

3790       NEXT  I 

3800  Pmin=0 

3810  MAT  Tb=B 

3820  MAT  A  =  IEP 

3830  IF  Err>Errmax  THEN  392© 

3840  DISP 

3850  FOR  C-l  TO  4 

3860  BEEP 

3870  WAIT  98 

3830  BEEP 

3890  WAIT  15W 

3900  NEXT  C 

3910   RETURN 

3920  G03UB  T«iip_print 

3930  IF  CItt>38)  OR  CErr>200)  OR  <Tout>400>  THEN  GOTO  Bomb 

3940  GOSUB  Calc2 

3950  GOTO  Solw*  !  END  OF  SOLVE 

3960  ! 

3970  Set_up: !  THIS  SECTION  SETS  UP  THE  Nreg  SIMUL  EQUATIONS  IS  MATRIX  FORM 

3980  FOR  1=1  TO  Nreg 

3990  IF  Tem_iol=0  THEN  B<  I >»Te<  I  J^Rtote  b<  I  >+Tai  r<  I  )-'Rb_conv<  I  ) 

4O0O  IF  Te»~sol*l  THEM  B<I)»«:Pow<I)*R«_conv<I)+Tair<I))~R*le<I)+Tair<I)/'Rb_conv< 

I) 

4010  L«dge*Redge»Tedge«Bedge=l 

4020  IF  <I-1>  MOD  Nxr»0  THEN  Ledge*2 

4030  IF  I  MOD  Nxr«0  THEN  Redge=2 

4040  IF  K*Nxr  THEN  T*cjg*  =  2 

405O  IF  I+Nxr  *Nreg  THEN  Bedge«2 

4060  IF  tl-i>  MOD  N  r  =  0  THEM  4080 

4070  A<I,  1-1  *=-l  'RK4,  I  ■ 

4030  IF  I  MOD  Nxr»0  THEN  4100 

4090  A(I,  I  +  n=-l  PI  .'2,  I> 

4100     IF     lONxr    THEN    4  120 

4110    A(I,  I-Nxr  >  — l'Rl  <3,  I) 

4120  IF  [+Nxr>»Nreg  THEN  4140 

4130  FKI,  I+Nxr >=-l  PI  <  1,  I) 

4140  AC  I, I)»<B«dg**l)/Rl CI, I)+<Redge«l >/Rl<2, !>+<Tedge=l >-Pl (3, I >+<L*dge=l >/R1  (4 

,  I  >*l/Rtote  b<  I)  +  l/  Rb  cor.v<  I)-<T«m  ao1  =1  >*Re  conv<  I  :>■'<  R*  1  e  < I>*Rtote_b<  i )  -■ 

4150  NEXT  I 

4  160  RETURN  !  END  OF  SETJJP 

4170  ! 

4130  Elu:    !  THIS  SECTION  PERFORMS  A  LU  DECOMPOSITION  OF  THE  'A 'MATRIX 

4190  DISP  "WORKING  ON  ITERATION  NUMBER  "Jltt 

4200  N»l*Nreg-l 

4210  FOR    K  =  l     TO    Nnil 

4220  .    Kpl=K+l 

4230  FOR  I=Kpl  TO  Nreg 

4240  G  =  -A<  I  ,  K  )--A';K,K> 

4250  A<I,K)=G 

4260  FOR  J  =  Kpl  TO  Nreg 

4270  A<  I,  J)=fl<I,  J)+G*fKK,J> 

4280  NEXT  J 

4290  NEXT  I 

4300      NEXT  K 
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4313  Solder:        !  THIS  SECTION  SOLVES  THE  NEW  MftTRIX  AND  PLACES  THE  ANSWERS 

4320  Npl=Nreg+l    '  INTO  MATRIX  B  TO  BE  PASSED  BACK  TO  SOLVE 

4330      FOR  K=l  TO  Nml 

4340      Kpl=K+l 

4350  FOR  I=Kpl  TO  Nreg 

4360  B<I>»B<I>+H<I,IO*B<IO 

4370  NEXT  I 

4330      NEXT  K 

4390  B<Nreg>*B<Nreg)  fl<Nreg,Nreg) 

4400      FOR  K=2  TO  Nreg 

4410      I=Npl-K 

4420      Ji-I+1 

4430  FOR  J= J 1  TO  Nreg 

4440  BCI>»B< I)-fl<  l,J>*BCJ) 

4450  NEXT  J 

4460      B(I)=B(r"ft>l,I) 

4470      NEXT  K 

4480  RETURN  !  END  OF  SOLVER 

4490  T*«p  print:!  THIS  SECTION  USED  FOR  INTERMEDIATE  OUTPUT 

4500  EXIT  GRAPHICS 

4510  IF  <Temprt»l)  AMD  Cflns=PI)  THEN  PRINTER  IS  0 

4520       PRINT  SPA-' H?  S  "DATA  FOR  "  ;  Pi  ct  1  bl  *J  "     *  "  |  1 1 1  I  "  I  TTERAT  I  ON  " 

4530       PRINT  LIN<2) 

4540  FIXED  2 

4550  PRINT    "PEG    »  Tc  as*     ;  DegC?  Tjunc     (BegC)       Treg    <LS*2£2       Po'A' 

<U)       R  j-c  <\A'  C  '  " 

4560       FOR  1=1  TO  Nreg 

4570       PRINT  I;TfiB<12>,Te<I)-2?3;TRB<2S),Tj<I)-273;TRB<:41?,Tb<I>-273;TflB<53>, 

Pou<  I  )  ;  TAE'  62  i  ,  Rj_C  '  I  "> 

4580       NEXT  I 

4590       PRINT  LIN<2) 

4600       PRINT  "BOARD  THICKNESS=" J 1000*Thi ck  b;  "mm   AND  CONDUCTIVITY  *";Kb;"Wat, 

ts/M-K" 

4610       PRINT 

4620  GOTO  4680    '    PRINT  "  I        RB_CONV       RE  CONV        TE  Hb" 

4630       FOR  1=1  TO  5 

4643       PRINT  USING  4660; I , Rb  conv< I ) , Re  c orw< I > , Te< I )-273 , Hb< I ) 

4650       NEXT  I 

4660    IMAGE  DD,4X,4<4D.4D,4X) 

4670    PRINT 

4683       FIXED  4 

4690    PRINT  "COOLING  AIR  FLOW  OF  ":Fair;"M"?  p~r    SEC     VEL=  " ;  V  ai  r  ;  "  M  ■  Sec  " ,  "  (  "  ; 

Vair*3?'3;  "FTVS>",LIN<1) 

4733    PRINT  "INLET  AIR  TEMP= " ; T ai r-273 ; "  deg  C    OUTLET  AIR  TEMP«" ; Tout-273; "de 

g  C" 

4710  PRINT  LIN<2>, "LARGEST  DIFFERENCE  BETWEEN  ITTEPATIONS  ="; Err ;" *****" 

4720  PRINT  " 

"  ,  L  I  N  ■.  2  • 


4733       PRINTER  IS  16 

4740       STANDARD 

4750       RETURN 

4763  Debug:  Ans**"N"!  TEMP  DEBUG  FOR  RESISTANCES 

4773   RETURN  " 

4738  INPUT  "DO  YOU  WISH  TO  HAVE  A  LIST  OF  ALL  THE  RESISTANCES  PRINTED  (NO  or 

>?", AnsJ 

4793  IF  UPC*<Ans*Ci, 12>-"N"  THEN  RETURN 

4303  PRINTER  IS  3 

4313  FIXED  5 
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4320  PRINT  "    RECg     RE  CONE    RE  CONV    PTOTE  E     RTCiP  g       RB 

CONV" 

4330  FOR  1=1  TO  Nreg 

4840  !  PRINT  TflB<2>  ,  I ,  TfiB<4>  , Re  condC I ) ,  TflB<4>  ;  Re  convC I> , TflB<4> , Rtote_b< I > J TRB< 

5>;Rtop_r<I>;TfiB<4>;Rele<:i> 

4858   PRINT  TflB<2>, I ;Re_cond< I ) : Re  conu( I>;Rtote_b< I> ;Rt op_r< I > ; Rb  conu< I) 

4860  NEXT  I 

4870  PRINTER  IS  16 

4330  STANDARD 

4390  RETURN 

4960  Units:!  THIS  SECTION  CONVERTS  FROM  METER  TO  mm 

4910  Bdl=1000*Bdl !  HEPE  CORRECT  UNITS  FOP  USE  IN  GRAPHICS  AND  TO  REWORK  ANALYSIS 

4920  Bdh*1000*Bdh 

4930  MAT  X»<1000>*X 

4940  MAT  Y=< 1000>*Y 

4950  MAT  Le*<l000)*Le 

4960  MAT  Ue«<1000)*We 

4970  Thick  b*1000*Thick  b 

4930  A1=A1*1E6 

4990  RETURN  !  END  OF  Units 

5000  ' 

5010  Output*!  GOSUE  Error 

5020  Output:'  THIS  SECTION  OUTPUTS  TO  GRAPHICS  ON  A  BLANK  BOARD 

5030  ON  ERROR  GOTO  Outputs 

5040  PLOTTER  IS  "GRAPHICS" 

5050  GRAPHICS 

5060  MSCALE  0,10 

5070  C3IZE  2 

5030  LORG  5 

5090  MOVE  13, 130 

5100  LABEL  " air  flow >  " 

5110  CSIZE  3 

5120  BdhaBdh'Sca 

5138  BdlsBdl^Sca 

5140  IF  Sca=l  THEN  5220 

5150  LORG  6 

5160  LDIR  PI/2 

5170  MOVE  175.70 

5180  IF  ScaM  THEN  LABEL  "THIS  PICTURE  IS  1^"&VAL*< Sca>&"  SIZE" 

5190  IF  Sca<l  THEN  LABEL  "THIS  PICTURE  IS  2X  SIZE" 

52O0  LORG  5 

5210  LDIR  0 

5220  MOVE  100-LEH'  "OUTPUT  DATA  FOP  »&Pi  c t  1  bl  *>/-2,  135 

5230  LABEL  USING  "K"; "OUTPUT  DATA  FOR  "&Pictlbl$ 

5240  MOVE  0,0 

5250  DRAW  0,Bdh 

5260  DPAU  Bd1 ,  Bdh 

5270  DRAW  Bdl ,0 

5230  DRAW  0,0 

5290  LINE  TYPE  3 

5300  Lr  =  Bdl/N-:r  !  LR  =  LENGTH  OF  EACH  REGION 

5310  Hr«Bdh/Nyr  !  HR  =  HEIGHT  OF  EACH  REGION 

5328     FOR  1=1  TO  Nxr-1 

5330     MOVE  I*Lr,0 

5340     DRAW  I#Lr,Bdh 

5350   NEXT  I 

5368  FOR     1=1     TO    Nvr-1 

5370  MOVE    0, I*Hr 
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5380 

3390 

5430 

5410 

5420 

5430 

5440 

5450 

5460 

5470 

5480 

5490 

5590 

5510 

5520 

5530 

5540 

5550 

5560 

5570 

5580 

5590 

5600 

5610 

5620 

5630 

5640 

5650 

5660 

5670 

5630 

5690 

5700 

5710 

5720 

5730 

5740 

5750 

5760 

deg  C 

5770 

5780 

IONS 

5790 

ONS  » 

5800 

5810 

5820 

5830 

CHRSC 

5340 


DRAW  Bdl , I*Hr 
NEXT  I 

LIME  TYPE  1 

LORC  3 

Nr  =  0 

CSIZE  2.4 

FOR  1=1  TO  Nreg 

MOVE  X< t>-.35*Lr,Y<I>+.35*Hr 

LABEL  USING  "K";VAL*<I> 

MOVE  X<I>+.3*Lr,Y< I)+.2S*Hr 

LABEL  USING  "K";Type*<Itype<I>) 

IF  I  type  U. '=0  THEN  5570 

IF  Tj<I>-273<.95*<Ti»ax-273>  THEN  L 

IF  Tj<I)-273>*.95*<Tmax-273)  THEN 

IF  Pow<I>O1.05*Pmin  THEN  LABEL  US 

IF  Pow<I>>1.05*Pmin  THEN  LABEL  USI 

LABEL  USING  5618; Te< I >-273 

GOTO  5590 


ABEL 
LABEL 
ING  5 
NG  56 


J3ING 
USINC 


20; Pot 


5610;  Tj  < 

5600; Tj 

..  <  I  > 

:  I  > 


LABEL  USING  "K"; " " 
LABEL  USING  5610;Tb<I) 
NEXT  I 
IMAGE  "+*",DDD.D,"  C" 
IMAGE  DDD.D. "  C" 
IMAGE  D.DDD, "  M" 
IMAGE  »**",  D.DDD,"  14" 
FIXED  4 
DUMP  GRAPHICS 
PRINTER  IS  0 
PRINT  SPA<20>,  " 
PRINT  SPA<20>, 
PRINT  SPA<28>, 
PRINT  SPA  20  )  . 
PRINT  SPA..  20  i  . 
PRINT  SPA  i.  20/, 
PRINT  SPfi<20>, 
PRINT  LIN'  2  i 
PRINT  SPA':  14  > 
PP  I  NT  US  IN 
';  rout-273; "deg  C" 
STANDARD 
IF  T«m_sol=0  THEN  PRINT  USING  5810; "LARGEST 

#"; Itt-i; "  and  *";  Itt;"  =  " ; Err 

IF  Tem_SOl»l  THEN  PRINT  USING  53  1  O  ;"  LARGEST 

" ; Itt-1 ; "  and  #"; Itt; "  =  ";Err 

IMAGE  1 1A, 2X, . 4D.7A, 2X, DD . 4D , 5A , 2X , 3D . D. 5A 
IMAGE  .'/45A,DD,*A,  DD,  A.DD.4D 

!  MOPE  OF  THE  CIRCUIT  DESCRIPTION  CAN  BE  P 
PRINT  LIN<1), "CIRCUIT  BOARD  DESCRIPTION  IS 
128) 
PRINT  " 


TYPE 

T  j  un  c 

POW 

Tcise 

FLOW  PATE 

5380; "COOL 


A  IP 


■■■'ELOCIT' 
";  Fair- 


Tout 


c"  ;  Tiir-273; 


CHANGE 

IN  POWER 

BETWEEN 

ITTERAT 

CHANGE 

IN  TEMP 

BETWEEN 

ITTEPATI 

. 2X,3D. I 

,  5A 

RESENTEI 

HEPE 

STORED 

UNDER  " ; 

CHR*< 132 

) ; Name*; 

5350 
5360 
5370 
5330 
5890 
5900 


",LINCJ 

PRINTER  IS  16 
EXIT  GRAPHICS 
RETURN 


What_ 
What," 


GOSUB  Error 
THIS  SECTION  PRESENTS  THE  VARIOUS  OPTH 


END  OF  OUTPUT 


'AILABLE  AND  DIRECTS 
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5910  !  PROGRAM  CONTROL  AS  NEEDED 

5920  ON  ERROR  GOTO  What_nowe 

3930  S*="WHAT  NOW  ?" 

3940  GOSUB  Pagehead 

595Q  PRINT  TflB<15),HYou  havt  completed  the  thermal  analysis  of    :" 

5960  PRINT  TflB<3?-LEN<Pict lbl  *)^2>;CHR*<  129)  ?  Pi  Ct  1  bl  *',  CHR*<  123) 

5970  PRINT  TflB<17>, "the  options  available  are  listed  below:" 

5980  PRINT  LIN<1>,TRB<5),  " 1.   CHANGE  COOLING  AIR  PARAMETERS  AND  RE-ANALYZE  " : CHR 

$< 132) ;Name*; CHR*< L23) : " . " 

5990  PRINT  LIN<l),TflB<r5),"2.   SENSITIVITY  ANALYSIS  FOR  CHANGES  IN  COOLING  AIR  PA 

RflMETERS  OF  " ;  CHR*«C  132)  ;  Name*;  CHR*<  123)  ;  "  .  " 

6000  PRINT  LIN<1),TRB<5), "3.   lifiKE  CHANGES  TO  BOARD  DESCRIPTION  IN  FILE  NAME  " ;  C 

HR*< 132); Name* ;CHR*<  128); " . " 

6010  PRINT  LIN<l),TflBC5),"4.   RETRIEVE  A  NEW  BOfiRD  DESCRIPTION  FROM  MASS  STORAGE 

6020  PRINT  LINC 1 ). TAB<5> . "5.   INPUT  A  NEW  BOARD  DESCRIPTION  FROM  THE  KEYBOARD." 

6030  PRINT  LIN<1>,TRB<5>, "6.   TERMINATE  SESSION" 

6040  Map=Bomb=u 

6050  Ansi="l" 

6060  INPuT  "YOUR  CHOICE  FROM  ABOVE  ( 1 , 2, 3, 4, 5, 6)?" , fins* 

6070  IF  Rns*»"BRCK  UP"  THEN  RETURN 

6080  Rna»:NT<VRL<flns*)) 

6090  IF  Cfins>0)  AND  (flns<»6)  THEN  6120 

6100  GOSUB  Err  in 

6110  GOTO  Mhat_nou 

6120  ON  Ana  GOSUB  Redo, Sensi , Change, Get  new, Key  new, Termi nate 

6130  IF  Rns-i  THEN  RETURN 

6140  GOTO  What_now 

6130  Redo:  RETURN  !  THIS  WILL  ALLOW  RESTART  OF  THIS  PROGRAM  WITH  SAME  BOARD 

6160  ' 

6170  Sensie: GOSUB  Error 

6180  Sensit!  THIS  SECTION  PRODUCES  PLOTS  OF  OUTPUT  VS  COOLING  AIR  PARAMETERS 

6190  ON  ERROR  GOTO  Sensie 

6200  S*«"SENSITIVTY  RNflLYSIS" 

6210  GOSUB  Pagehead 

6220  PRINT  "     This  section  allows  you  to    investigate  the  effects  of  variations 

in  the  air  flow  rate.   "; 
6230  IF  Tem_sol=l  THEN  PRINT  "Plots  of  Maximum  Junction  Temperatura  vs  Flow  Rate 

of  the  air  are  pr o d u c e d " 
6240  IF  Tem_sol*S  THEN  PRINT  "Plots  of  Minumum  Power  vs  Flow  rate  of  the  air  are 

produced. " 
6250  PRINT  "     You  specify  the  maximum  flow  rate  p^r    board  (M"3^sec)  and  five  s 
eparate  " 

6260  PRINT  "analyses  are  performed  and  the  results  plotted.   NOTE:  selecting  a  m 
ax i mum  flow" 

6270  PRINT  "»hat  is  evenly  divisable  by  five  (5,10,30)  will  result  in  better  loo 
king  axes. ",LIN<1) 
6280  PRINT  "These  plots  may  be  produced  on  either  the  screen  (with  hard  copy  via 

Key3>  or  " 
6290  PRINT  "on  a  peripheral  plotter  such  as  the  Hp  9S72A.   In  addition  a  printed 

t  abul ar " 
6300  PRINT  "output  of  the  results  at  each  of  the  airflow  rates  may  be  produced." 
6310  PRINT 
6320  flns*="NULL" 

6330  INPUT  "WHAT  IS  THE  UPPER  LIMIT  ON  THE  AIR  FLOW  RATE  FOR  THE  SENSITIVITY  ANA 
LYSIS?", flns* 

6348  IF  flns*="BflCK_UP"  THEN  What_now 
6350  IF  flns*-"NULL"  THEN  6320 
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6360  flnssRBS<VflLCflns*>) 

6370  Xmax=Ans 

6380  Xmin=0 

6390  Xstep*<X»ax-Xmin) -  5 

6400  Ymm  =  0 

6410  Mult=100OO0 

6420  IF  Xmax>».001  THEN  Mult =10000 

6430  IF  Xmax>=.01  THEN  Multal000 

6440  IF  Xmax>=.l  THEN  Mult»100 

6450  AnsS="N" 

6460  INPUT  "DO  YOU  WISH  to  HAVE  TABULAR  RESULTS  FOR  EACH  OF  THE  AIR  FLOWS?<N  or 

Y)",Ans* 

6470  IF  Ans-*="BACK_UP"  THEN  Sensi 

6480  Temprt=0 

6490  IF  UPC*<fins*>=«"Y"  THEN  Temprt=l 

6500  GOSUB  Database 

6510  Map=t 

6520  IF  Tern  sol'l  THEN  6550 

6530  FOR  F*Tr«Xmax  TO  Xstep  STEP  -Xstep 

6540  GOTO  6560 

6550  FOP  Fair*Xstep  TO  Xmax  STEP  Xstep 

6560  GOSUB  CalcQ 

6570  GOSUB  Calcl 

6530  GOSUB  Cak_t 

6590  GOSUB  Calc2 

660O  GOSUB  Ca1c_air 

6610  GOSUB  Solve 

6620  IF  Tem_sol=l  THEN  Y*Tmax-273 

6630  IF  Tem"~so1=0  THEN  Y=Pmin 

6640  IF  Temprt=l  THEN  PRINTER  IS  0 

6650  IF  Temprt,  =  1  THEN  GOSUB  Temp_print 

6660  IF  Map=l  THEN  GOSUB  Plot. 

6670  IF  Rns*="BflCK_UP"  THEN  6450 

6680  Rns*Itt«Map=0 

6690  GRAPHICS 

67O0  PEN  1 

6710  MOVE  Fair.Y 

6720  LORG  5 

6730  LABEL  USING  "K": "*" 

6740  PEN  0 

6750  NEXT  Fair 

6760  EXIT  GRAPHICS 

6770  DUMP  GRAPHICS 

6780  PRINTER  IS  0 

6790  PRINT  LIN< 1 i 

6800  PRINTER  IS  16 

6810  GOSUB  Units 

6820  GOTO  What_now  !  END  OF  SENSI 

6830  ! 

6340  Change: !  THIS  SECTION  LOADS  BOARDS  TO  ALLOW  CHANGES  TO  THE  CIRCUIT  DESCPIPT 

6850  Map=l 

6360  PRINT  PAGE 

6870  DISP  "WORKING  LOADING  BOARDS" 

6830  LOAD  "BOARDS", 1 

6890  Get  new: !  THIS  SECTION  EXPLAINS  THE  TWO  WAYS  TO  GET  A  NEW  BOARD  DESCRIPTION 

6900  S*=" INPUT  OF  NEW  BOARD  DESCRIPTION" 

6910  GOSUB  Page  he  id 

6920  PRINT  "     The  new  board  description  may  be  read  in  from  nasi  storage  in  tw 

O  different" 
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6930  PRINT  "programs.   THERML  'the  program  in  core  now)  will  not  allow  visual 

Checking  or" 

6940  PRINT  "modi fi c at 1 on  of  descriptive  data  for  the  circuit  board.   The  data  is 

read  in" 
6953  PRINT  "at  your  direction  but  only  question*  concerning  the  environment  exte 
rior  to  the" 

6969  PRINT  "circuit  board  are  asked.   The  second  method  of  input  from  mass  stora 
ge  indoles" 

6970  PRINT  "the  program  BOARDS  which  allows  both  data  checking  through  graphic 
s  and  " 

6930  PRINT  "opportunities  to  modify  the  descriptive  data." 

6990  PRINT 

7000  PRINT  LIN<l),TflB<5>,"l.   READ  A  NEW  BOARD  DESCRIPTION  USING  THERML.  (no  c 

orrect  i  ons) " 

7010  PRINT  LIH<l>,TflB<5), "2.   READ  A  NEW  BOARD  DESCRIPTION  USING  BOARDS.  Callo 

us  corrections)" 

7020  PRINT  LIN<2>,"  REMEMBER  THE  THERMELEX  SYSTEM  MUST  BE  IN  THE  DEFAULT  MASS  ST 

ORAGE  DEVICE" 

703O  Ans*="l" 

7040  INPUT  "YOUR  CHOICE  FROM  ABOVE  <  1  ,  2  )  ?  "  ,  AnsS 

7050  IF  Ans*="BACK_UP"  THEN  What_now 

7060  Ans  =  VAL<Ans-J) 

7070  IF  (Ans>0)  AND  (flnsO)  THEN  7100 

708O  GOSUB  Err in 

709O  GOTO  Get  new 

7100  IF  Ans=2  THEN  71 30 

7110  Map=l 

7120  RETURN 

7130  ! 

7140  Map=2 

7150  DISP  "WORKING  LOADING  BOARDS" 

7160  LOAD  "BOARDS",  1 

7170  ! 

7130  Key  new: !  THIS  SECTION  LOADS  BOARDS  WITH  THE  INTENT  TO  INPUT  HEW  BOARD  DESC 

7190  Map=3 

7200  DISP  "WORKING  LOADING  BOARDS" 

7210  LOAD  "BOARDS", 1 

7220  ! 

7230  Terminate:   !  THIS  SECTION  TERMINATES  THE  SESSION 

7240  GCLEAR 

7250  DISP  "WORKING  LOADING  STANDARD  KEY  DEFF I N I T IONS " 

7260  LOAD  KEY  "STDKEY" 

7270  PRINT  LIN<23),3Pfl<13>, "NORMAL  TERMINATION" 

7280  PRINT  LIN<2>,SPfl<15>, "     THANK  YOU 

7290  DISP 

7300  END  !END  OF  TERMINATE 

7310  ! 

7320  Pagehead: !  THIS  ROUTINE  PLACES  THE  PAGE  HEADINGS  FOR  THE  INSTRUCTIONS 

7330  PRINT  PAGE, TflB<34-LEH<S*>  '2>, "**  " ; CHR* (1 32 > ; S* ; CHR* <  1 23 ) ; "  **",LIN<2) 

7340  RETURN 

7350  ' 

7360  Errir,:'  THIS  SECTION  ALERTS  THE  USER  TO  AN  ATTEMPT  TO  INPUT  BAD  DATA 

7370  BEEP 

7330  DISP  "******  INPUT  OUT  OF  RANGE TRY  AGAIN" 

7390  WAIT  1500 

7400  BEEP 

7410  RETURN  !  END  OF  ERR  IN 

7420  i 
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7430  Error: !  THIS  SECTION  IS  THE  ERROR  TRAPPING  ROUTINE  FOR  THE  ENTIRE  PROGRAM 

7449  !  PROGRAM  FLOW  RESUMES  AT  THE  TOP  OF  THE  SECTION  IN  WHICH  THE  ERROR 

7450  !  OCCUPED  AFTER  THE  USER  PRESSES  CONT 
7460  EXIT  GRAPHICS 

7470  PRINTER  IS  16 

7480  PRINT  PAGE 

7490  BEEP 

7500  WAIT  30O 

7510  BEEP 

7520  IF  EPRN=56  THEN  Err_nam« 

7530  PRINT  LIN<20>,SPA<10>, "ERROR  NUMBER" J ERRNJ "HAS  OCUPRED  IN  LINE" ; ERRL; " .   PR 

ESS  CONT  WHEN  READY" 

7540  DISP 

7559  BEEP 

7560  PAUSE 
7570  RETURN 
7530  ! 

7590  Err  name: !  THIS  SECTION  FOR  IMPROPER  FILE  NAME 

7600  PRINTER  IS  16 

7610  PRINT  PAGE 

7620  Msus*=" DEFAULT  MASS  STORAGE" 

7630  FOR  1=2  TO  LEN<Name* • 

7640  IF  Nam«*CI, I3»": "  THEN  7670 

7650  NEXT  I 

7660  GOTO  7700 

7670  Hsus*=Name*C I ] 

7680  CAT  Msus* 

7690  GOTO  7710 

7700  CAT 

7710  PRINT  LIN<1.,  "FiU  Name  ■ ;  CHR*<  132)  ;  NamefC  1 ,  I  -1  3  ;  CHR*<  123)  ;  "  is  NOT  on  " ;  CH 

R*<132);Hsus*;CHR*<128):  "  '.nth  that  spelling " 

7720  PRINT  LINC1), "CHECK  OVER  THE  DIRECTORY  ABOVE  FOR  CORRECT  NAME  OR  SPELLING.. 

7730  DISP  "PRESS  CONT  WHEN  READY" 

7740  PAUSE 

7750  RETURN  !  END  OF  ERROR 

7760  Plote:GOSUB  Error 

777Q    Plot:  i  PLOTTING  ROUTINE  FOP  THE  AXES 

7780  ON  ERROR  GOTO  Plot* 

7790  PLOTTER  IS  "GRAPHICS" 

7800  IF  T«m_so1*0  THEN  7830 

7810  Yitep=10 

7820  FOR  1=0  TO  7 

7830  IF  Tmax-273>50+I*25  THEN  Ystep»15+I*5 

7840  NEXT  I 

7S50  Pltlb1*="Tjunc  vs  Air  Flow" 

7860  Ylbl*»"Junc  Temp  (deg  C>" 

7370  GOTO  7940 

7830  Ystsp=. 1 

7890  FOR  1=1  TO  10 

7900  IF  Prnin>.5*I  THEN  Ystep*.  1*<  1  +  1) 

7910  NEXT  I 

7920  Plt1bl*»"Pmin  ui  Air  Flow" 

7930  YlblJ="Mir  Power/Cornp  CWatts)" 

7940  Y»ax=5*Ystep 

7950  GRAPHICS 

7960  LOCATE  15,120,10,95 

7970  SCALE  Xmi n, Xmax, Y»i n, Ymax 
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7938  AXES  Xstep,Ystep,Xmin,Y»in 

7993  CSIZE  3 

8000  LDIR  0 

8810  LORG  6 

8020  FOR  Xpos=Xmin  TO  Xmax  STEP  Xst«p 

3038  MOVE  Xpos, Ymin-. l#Ystep 

8048  LABEL  USING  "K" ; Xpos*Mul t 

3058  NEXT  Xpos 

8068  MOVE  2.5*Xstep,-.4*Yst*p 

8079  LABEL  USING  '  k  " :  "  A  i  r-F  1  ow  Board  <  "&VF»L*<  l'Mult  >&"  MA3^Sec)" 

8038  LORG  3 

8890  FOR  Ypos*Ymin  TO  Ymax  STEP  Ystep 

8108  MOVE  Xmin-. l*Xstep, Ypos 

8110  LABEL  USING  "K";Ypos 

8128  NEXT  Ypos 

8138  LINE  TYPE  1 

3148  LORG  4 

3158  LDIR  PI--2 

8168  MOVE  -.4*Xst*p,2.5*Ystep 

8178  LAEEL  USING  "K-jYlbl* 

8188  MOVE  Xmin+2.3*Xstep,5. l*Ystep 

8190  CSIZE  4 

8280  LDIR  3 

3218  LABEL  USING  "K";Pictlb1* 

8220  LABEL  USING  "K";Plt Ibl* 

8238  RETURN  !  END  OF  PLOT 

8248  ! 

8258  Bomb:'  THIS  SECTION  DELIVERS  MESSAGE  TO  THE  USER  OF  FAILURE  TO  CONVERGE 

8268  WAIT  2000 

3270  BEEP 

8288  PRINT  PAGE  ,  L  I  N  <  1  8  ;•  ,  "  UNABLE  TO  ACHIEVE  CONVERGENCE  DUE  TO  NUMERICAL  IN3TABIL 

ITIES" 

3298  PRINT  LINC3),"I  SUGGEST  A  CHANGE  IN  EITHER  THE  INSTALLATION  PARAMETERS  OR  :' 

8300  PRINT  LIN' 1', "THE  CIRCUIT  BOARD  PARAMETERS  ...." 

8318  PRINT  LIN<2), "THERMELEX  PREDICTS  TEMPERATURES  MUCH  MUCH  BETTER  THAN  POWER  L 

EVELS" 

8328  PRINT  LINc  1),"TRY  SPECIFYING  THE  COMPONENT  POWER  LEVELS .  "  , L I N <  3 > 

3338  DISP  "Press  CONT  when  ready  to  return  to  option  list" 

8348  BEEP 

8350  PAUSE 

8360  DISP 

8378  Bomb=l 

8388  RETURN  !  END  OF  BOMB 
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